How to find solution for Blasius Equation?

2 次查看（过去 30 天）
Aruna P
Aruna P 2016-3-13
评论： Geoff Hayes 2017-4-22
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?

请先登录，再回答此问题。

采纳的回答

Torsten
Torsten 2016-3-14
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.
  8 个评论
显示 6更早的评论
Aruna P
Aruna P 2016-3-15
i tried using "vpasolve" but i am getting error
Torsten
Torsten 2016-3-15
I don't understand what you are trying to do.
To print the results to a file, use fprintf after the command
plot(xint,Sxint([2 2],:));
Best wishes
Torsten.

请先登录，再进行评论。

更多回答（1 个）

Maria Hamdani
Maria Hamdani 2017-4-22
what will i do if i have to print table of f,f'and f'' for different values of eta
  1 个评论
Geoff Hayes
Geoff Hayes 2017-4-22
Maria - please post your above question as a question rather than as an answer.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Mathematics 的更多信息

标签

尚未输入任何标签。

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by