How to draw fig.1 with this code

MINATI

2019 5 25

0 个回答

8 次查看（30 天）

0 个投票

shrinking.pdf

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

3 个评论
显示 1更早的评论隐藏 1更早的评论

MINATI PATRA 2019-5-26

Actually Fig.1 is not coming properly. I am getting straight lines instead of curves as in Fig.1 Any help

Rik 2019-5-26

http://matlab.wikia.com/wiki/FAQ#Can_you_program_up_the_algorithm_in_this_article_for_me_and_explain_it_to_me.3F

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.

请先登录，再进行评论。

请先登录，再回答此问题。

Follow Question