MATLAB Mathematics - SERC - Index of

Boundary Value Problems for ODEs
Falkner−Skan equation, positive wall shear, β = 0.5.
1.2
1
0.8
df/dη
0.6
0.4
0.2
0
0 1 2 3 4 5 6
η
5 Use continuation to solve the problem and plot subsequent solutions.
The example then solves the problem for infinity = 4, 5, 6. It uses bvpinit
to extrapolate the solution sol for one value of infinity as an initial guess
for the next value of infinity. For each iteration, the example prints the
computed value of f′′ ( 0)
and superimposes a plot of the solution in the
existing figure.
for Bnew = infinity+1:maxinfinity
end
solinit = bvpinit(sol,[0 Bnew]); % Extend solution to Bnew.
sol = bvp4c(@fsode,@fsbc,solinit);
eta = sol.x;
f = sol.y;
fprintf('Value computed using infinity = %g is %7.5f.\n', ...
Bnew,f(3,1))
plot(eta,f(2,:),eta(end),f(2,end),'o');
drawnow
5-79