MATLAB: I want to know, how to solve the Falkner Skan equation f´´´+ff´´+​b(1-(f´)sq​uare)=0 2D plot

2d plot

Best Answer

syms b f(t) t Y
Df = diff(f);
D2f = diff(Df);
D3f = diff(D2f);
Eq = D3f + f*D2f + b*(1 - Df^2) == 0;
[VF,Sbs] = odeToVectorField(Eq);
odefcn = matlabFunction(VF, 'Vars',{t,Y,b});
tv = [0, 5];
ics = [0, 0, 0];
b = -42; % Choose A Value
[T,Y] = ode45(@(t,Y)odefcn(t,Y,b), tv, ics);
figure
plot(T,Y)
grid
Related Question