Please someone correct the error
Show older comments
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
Answers (1)
Walter Roberson
on 25 Jul 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.
Categories
Find more on Ordinary Differential Equations in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!