MATLAB Mathematics - SERC - Index of

5 Differential Equations
3 Solve on the initial interval. The example obtains the solution for
infinity = 3. It then prints the computed value of f′′ ( 0)
for comparison
with the value reported by Cebeci and Keller [2]:
sol = bvp4c(@fsode,@fsbc,solinit);
eta = sol.x;
f = sol.y;
fprintf('\n');
fprintf('Cebeci & Keller report that f''''(0) = 0.92768.\n')
fprintf('Value computed using infinity = %g is %7.5f.\n', ...
infinity,f(3,1))
The example prints
Cebeci & Keller report that f''(0) = 0.92768.
Value computed using infinity = 3 is 0.92915.
4 Setup the figure and plot the initial solution.
figure
plot(eta,f(2,:),eta(end),f(2,end),'o');
axis([0 maxinfinity 0 1.4]);
title('Falkner-Skan equation, positive wall shear, \beta = 0.5.')
xlabel('\eta')
ylabel('df/d\eta')
hold on
drawnow
shg
5-78