how to plot 3rd order differential equation in Matlab?

Question

Joy Salomi 2020-2-13

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/505266-how-to-plot-3rd-order-differential-equation-in-matlab

评论： Joy Salomi 2020-2-16

Hello,

What is wrong in this coding? I got some error. Please help me in this. I have attached the problem and code below. Thanks in advance.

My equation and bcs are:

f'''+f f'' -1.3333 (f')^2 - f' = 0

f(0)=0 , f'(0)=1 , f'(infinity)=0

function heattransfernew
options=bvpset('stats','off','RelTol',1e-6);
solinit=bvpinit(linspace(0,6,20),[0,6]);
sol=bvp4c(@f, @OdeBC, solinit, options);
plot(sol.x,sol.y(1,:));
end
  function dydx = f(x,y)
dydx(1) = y(2);
dydx(2) = y(3);
dydx(3) = -y(1).*y(3)+1.3333.*(y(2).^2)+y(2);
 end
function res = OdeBC(ya,yb)
res = [ya(1)
    ya(2)-1
    yb(2)
    ];
end