Please someone correct the error
Afficher commentaires plus anciens
maxellfluid
function maxellfluid
beta=1.5;
gamma=0.25;
Pr=3.5;
lam=0.2;
E=0.5;
sol=bvpinit(linspace(0,5,50),[0 0 0 0 0 0 0]);
sol1=bvp4c(@bvpexam2,@bcexam2,sol);
x1=sol1.x;
y1=sol1.y;
figure (1)
plot(x1,y1(2,:),'b'); %plot of prime
hold on
figure(2)
plot(x1,y1(4,:),'--'); %plot of g
hold on
figure(3)
plot(x1,y1(6,:),':'); %plot of theta
hold on
function res=bcexam2(y0,yinf)
res=[y0(1)-0;y0(2)-1;y0(4)-0;y0(6)-1;yinf(2)-0;yinf(4)-0,yinf(6)-0];
end
function ysol=bvpexam2(~,y)
yy1=(1/(1-beta*y(1)^2))*(-y(1)*y(3)+y(2)^2-2*lam*(y(4)-beta*y(1)*y(5))-2*beta*y(1)*y(2)*y(3));
yy2=(1/(1-beta*y(1)^2))*(-y(1)*y(5)+y(2)*y(4)+2*lam*(y(2)+beta*(y(2)^2-y(1)*y(3)+y(4)))-2*beta*y(1)*y(2)*y(5));
yy3=(1/(1+E*y(6)-Pr*gamma*y(1)^2))*(-E*y(7)^2-Pr*y(1)*y(7)+Pr*gamma*y(1)*y(2)*y(7));
ysol=[y(2);y(3);yy1;y(5);yy2;y(7);yy3];
end
end
Réponses (1)
Walter Roberson
le 25 Juil 2022
res=[y0(1)-0;y0(2)-1;y0(4)-0;y0(6)-1;yinf(2)-0;yinf(4)-0,yinf(6)-0];
^
You used a comma in the place marked, instead of a semi-colon.
Catégories
En savoir plus sur Ordinary Differential Equations dans Centre d'aide et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!