ODE Solver in matlab

HI

Answer 1

[darova]

Here is a start. But don't know how to write 3d and 4th initial guess

f = @(t,x) [x(2); -C2/C1];
fb = @(p0,pt) [(1-gamma)*pt(1)-(xt+1+lambda)*pt(2)
    (1-gamma)*p0(1)-(x0+1)*p0(2)
    % dont know how to write 3d
    % dont know how to write 4th];

solinit = bvpinit(0:0.1:5, [1 1 1 1]);      % t span and initial guess
sol = bvp4c(@f, @fv, solinit);              % solve equation