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

1 visualización (últimos 30 días)
MINATI
MINATI el 7 de Jun. de 2019
Comentada: Anil Gautam el 23 de Jul. de 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
%%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)
  3 comentarios
MINATI
MINATI el 16 de Jul. de 2019
ok Anil
you can change according to BC but graph should match previous study also.
Anil Gautam
Anil Gautam el 23 de Jul. de 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);
Thank you in advance

Iniciar sesión para comentar.

Respuestas (1)

KALYAN ACHARJYA
KALYAN ACHARJYA el 7 de Jun. de 2019
Editada: KALYAN ACHARJYA el 7 de Jun. de 2019
Keep figure, only, not figure(1)
  3 comentarios
MINATI
MINATI el 9 de Jun. de 2019
Any help will be greatly appreciated
MINATI
MINATI el 17 de Jul. de 2019
can I withdraw the code unsolved

Iniciar sesión para comentar.

Categorías

Más información sobre Time Series Objects en Help Center y File Exchange.

Etiquetas

Community Treasure Hunt

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

Start Hunting!

Translated by