策划和For循环ODEFUN BVP4C初始猜测
43岁的观点(30天)
显示旧的评论
% % Thetadot(0)和公关
m = 0;
Prinf = 1000;
为Prinf公关= linspace (0, 1)
solinit = bvpinit(公关,[0 0 0 0 0.05]);
索尔= bvp4c (@odefun @odefun_bc solinit);
Prinf xint = linspace (0, 1);
Sxint =德瓦尔(溶胶,xint);
图(19)
持有在
标题(“HeatFlux(0)和公关的)
包含(“公关”)
ylabel (“热流”)
情节(xint Sxint (5、1));%情节qdot (0)
结束
我想阴谋y(5、1)在不同的公关;然而似乎我最初的猜测是不够的。
接受的答案
Torsten
2023年4月29日18:54
全球男公关
etainf = 20;%寻找临时和收敛速度
% % Thetadot(0)和公关
m = 0;
公关= 0:0.5:50;
为i = 1:元素个数(PR)
公关=公关(i);
solinit = bvpinit (linspace (0 etainf 100), [0 0 0 0 0.05]);
索尔= bvp4c (@odefun @odefun_bc solinit);
qdot0 (i) = sol.y (5、1);
结束
情节(公关、qdot0)
标题(“HeatFlux(0)和公关的)
包含(“公关”)
ylabel (“热流”)
函数yprime = odefun(η,y)
yprime = 0 (5、1);
全球男公关
% Blasius Eqn
yprime (1) = y (2);
yprime (2) = y (3);
yprime (3) = (1/2) * (m + 1) * (1) * y (3) + m * y(2) ^ 2米;
%能源Eqn
yprime (4) = y (5);
yprime(5) =(1/2) *公关* y (1) * (m + 1) * y (5);
yprime = yprime ';
结束
函数res = odefun_bc (ya,)
res =[丫(1);丫(2);yb (2) 1;丫(4);yb (4) 1);
结束
4评论
Torsten
2023年4月29日在7:30
看起来像一个+ b * sqrt (Pr)在我看来。
全球男公关
etainf = 20;%寻找临时和收敛速度
% % Thetadot(0)和公关
m = 0;
公关= 0:0.5:50;
为i = 1:元素个数(PR)
公关=公关(i);
solinit = bvpinit (linspace (0 etainf 100), [0 0 0 0 0.05]);
索尔= bvp4c (@odefun @odefun_bc solinit);
qdot0 (i) = sol.y (5、1);
结束
持有在
情节(公关、qdot0)
=[(元素个数(PR), 1) sqrt(公关。');
b = qdot0。”;
x = \ b;
情节(公关、x (1) + (2) * sqrt (PR))
持有从
函数yprime = odefun(η,y)
yprime = 0 (5、1);
全球男公关
% Blasius Eqn
yprime (1) = y (2);
yprime (2) = y (3);
yprime (3) = (1/2) * (m + 1) * (1) * y (3) + m * y(2) ^ 2米;
%能源Eqn
yprime (4) = y (5);
yprime(5) =(1/2) *公关* y (1) * (m + 1) * y (5);
yprime = yprime ';
结束
函数res = odefun_bc (ya,)
res =[丫(1);丫(2);yb (2) 1;丫(4);yb (4) 1);
结束