# Figs. 1, 2, 3 are not coming while using the following code

2 views (last 30 days)
MINATI on 7 Jun 2019
Commented: Anil Gautam on 23 Jul 2019
Pr=1;Sc=1;
R=0.1;
xa=0;xb=6;
Lv = [-2.5:0.1:0];
for D=[ 0 0.5 1]
for i=1:length(Lv)
L = Lv(i);
fODE = @(x,y) [y(2); y(3); y(2)^2-y(3)*y(1)-1; y(5); -3*Pr*y(1)*y(5)/(3+4*R); y(7); -Sc*y(1)*y(7)];
BCres= @(ya,yb)[ya(1); ya(2)-L-D*ya(3); ya(4)-1; ya(6)-1; yb(2)-1; yb(4);yb(6)];
xint=linspace(xa,xb,101);
solinit1=bvpinit(xint,[0 1 0 1 0 1 0]);
solinit2=bvpinit(xint,[-1 -1 -1 -1 -1 -1 -1]);
sol1=bvp4c(fODE,BCres,solinit1);
sol2=bvp4c(fODE,BCres,solinit2);
sxint1=deval(sol1,xint);
sxint2=deval(sol2,xint);
Sk1(i,1) = sxint1(3,1);
Sk2(i,1) = sxint2(3,1);
end
figure(1)
plot(Lv,Sk1,'-',Lv,Sk2,'--','Linewidth',1.5);
xlabel('\bf \lambda');
ylabel('\bf C_{f}');
hold on
end
%% while using the code, Figs. 1, 2, 3 are not coming as curves, it came as lines.
%% Another help required is how to give a DOT mark at the end of 1st solution (as in PDF)
Anil Gautam on 23 Jul 2019
In given line
solinit1=bvpinit(xint,[0 1 0 1 0 1 0]);
why you take [0 1 0 1 0 1 0] ??
can anyone take [0 0 0 0 0 0 0] ??
or [0 1 1 0 1 1 0] ??
Besically i want to understand the role of A in following line
solinit1=bvpinit(xint, A);

KALYAN ACHARJYA on 7 Jun 2019
Edited: KALYAN ACHARJYA on 7 Jun 2019
Keep figure, only, not figure(1)
MINATI on 17 Jul 2019
can I withdraw the code unsolved