How to draw fig.1 with this code

1 view (last 30 days)
MINATI on 25 May 2019
Commented: Rik on 26 May 2019
Pr=1;Sc=1;
R=0.1;
xa=0;xb=6;
Lv = [-2.5:0.1:0];
for D=[ 1 1.5 2]
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

Rik on 26 May 2019
What is your question? Are you successfully implementing the paper, but is your problem with the plotting?
Note that you didn't use hold on, which means the axis is cleared every iteration.
MINATI PATRA on 26 May 2019
Actually Fig.1 is not coming properly. I am getting straight lines instead of curves as in Fig.1 Any help
Rik on 26 May 2019
Try to make your question more precise. I'm not going to read 4-5 pages of a pdf. Try to distil the problem to its essence. Have a read here and here. It will greatly improve your chances of getting an answer.