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

조회 수: 1(최근 30일)
MINATI 7 Jun 2019
댓글: Anil Gautam 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)
##### 댓글 수: 3표시숨기기 이전 댓글 수: 2
Anil Gautam 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);

댓글을 달려면 로그인하십시오.

### 답변(1개)

KALYAN ACHARJYA 7 Jun 2019
편집: KALYAN ACHARJYA 7 Jun 2019
Keep figure, only, not figure(1)
##### 댓글 수: 3표시숨기기 이전 댓글 수: 2
MINATI 17 Jul 2019
can I withdraw the code unsolved

댓글을 달려면 로그인하십시오.

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by