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

Question

MINATI 2019-6-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/465936-figs-1-2-3-are-not-coming-while-using-the-following-code

评论： Anil Gautam 2019-7-23

sol.pdf

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
%%Please help with the following questions:
%% 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)