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
