system of coupled ordinary differential equations

1 次查看（过去 30 天）
Priya M
Priya M 2021-1-23
评论： Priya M 2021-1-25
i need matlab code by using bvp4c for the above equations.....
(1+1/beta)f'''(eta)-2*(f'+g')f'+(f+g)f''-M^2f'=0
(1+1/beta)g'''(eta)-2(f'+g')g'+(f+g)g''-M^2g'=0
teta''+pr[(f+g)theta'-A(f'+g')theta+S*theta]=0
phi''+Sc[(f+g)phi'-B(f'+g')phi-r*phi]=0
boundary conditions
f(0)=0 ,f'(0)=1 ,g(0)=0, g'(0)=alpha ,theta(0)=1+r1theta'(0) ,phi(0)=1+r2phi'(0)
f'(infinity)=0 g'(infinity)=0,theta(infinity)=0,phi(infinity)=0
  2 个评论
Priya M
Priya M 2021-1-25
xlow=0;
xhigh=6;
solyinit=bvpinit(linspace(xlow,xhigh,2),[0,1.01,11,1,0,1,0,12,11.05,1.005]);
M=0.5;
beta=0.02;
B=0.5;
Pr=0.7;
A=0.5;
Sc=1.0;
S=0.1;
r=0.3;
r1=0.1;
r2=0.1;
dydx=@(x,y)[y(2);
y(3);
(((M^2*y(2))-2*y(2)^2-2*(y(2)*y(5)-y(1)*y(3)-y(4)*y(3))/(1+1/beta)));
y(5);
y(6);
(((M^2*y(5))-2*(y(2)*y(5))-2*y(5)^2-y(1)*y(6)-y(4)*y(6))/(1+1/beta));
y(8);
((1/Pr)*(A*(y(2)*y(7)-y(1)*y(8)-y(4)*y(8)+A*(y(5)*y(7))-S*y(7))));
y(10);
((1/Sc)*(B*(y(2)*y(9))+B*(y(5)*y(9))-y(1)*y(10)-y(4)*y(10)+r*y(7)))];
res=@(ya,yb)[ya(1); ya(2)-1; ya(4); ya(5)-alpha; ya(7)-r1*ya(8); ya(9)-r2*ya(10); yb(2); yb(5); yb(7); yb(9)];
fsol=bvp4c(dydx,res,solyinit);
X=fsol.x;
Y=fsol.y;
plot(X,Y)

请先登录，再进行评论。

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

回答（0 个）

类别

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

标签

Community Treasure Hunt

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

Start Hunting!

Translated by