MATLAB: How to find f”(0) where f=y(1),f’=y(2),f”=y(3)

MATLABres=[ya(1) ya(2)-1 ya(4) ya(5)-p ya(7)-1 yb(2) yb(5) yb(7)]

dydx=[y(2); y(3); 2*a1*y(2)*(y(2)+y(5))-a1*y(3)*(y(1)+y(4));
y(5); y(6); 2*a1*y(5)*(y(2)+y(5))-a1*y(6)*(y(1)+y(4));
y(8); A*Pr*a2*y(7)*(y(2)+y(5))-Pr*a2*y(8)*(y(1)+y(4))];
res=[ya(1); ya(2)-1; ya(4); ya(5)-p; ya(7)-1; yb(2); yb(5); yb(7)];

Best Answer

f0 = deval(sol,0);
f0(3)

Best wishes

Torsten.

Related Question