MATLAB: How to find solution for Blasius Equation

function blasius
xa = 0; xb = 7;
solinit = bvpinit(linspace(xa,xb,1000),[0 0 0]);
sol = bvp4c(@blasius_equation,@bc,solinit);
xint = linspace(xa,xb,1000);
Sxint = deval(sol,xint);
plot(xint,Sxint([2 2],:));
function res = bc(ya,yb)
res = [ ya(1); ya(2); yb(2)-1];
function dydx = blasius_equation(x,y)
dydx = [y(2); y(3); (-y(1)*y(3))/2];
I did not get the error. But if i change the line "plot(xint,Sxint([2 2],:));" to "plot(xint,Sxint([1 2],:));", i get two lines in the graph. What does this line in the program means?

Best Answer

plot(xint,Sxint([2 2],:));
plots y(2) over [xa xb].
I don't know why you or the author of the code does not simply use
plot(xint,Sxint(2,:));
Best wishes
Torsten.
Related Question