Hi, I solved correctly a system of ODE second order with BCs, with the instruction bvp4c.
However, I have some questions about it.
My code is as follows:
xmesh = linspace(1e-5,1.5*1e-3,20);
solinit = bvpinit(xmesh,1e-3*[0,0,0,0,1,1,1,1]);
sol = bvp4c(@bvpfcn,@bcs,solinit);
1) When I plot the solution, I should choose in x interval sol.x or xmesh? What is their difference?
2) I don't understand how Matlab uses the second vector [0,0,0,0,1,1,1,1] in bvpinit.
It should represent the initial guess of solution (I have 8 variables), but it seems that its values don't change final solution. Is it right?
I've already read Matlab guide about bvp4c, but I have not found these answers.