previous technique didnt work for the present code

1 visualización (últimos 30 días)
MINATI
MINATI el 24 de Feb. de 2019
Comentada: Walter Roberson el 27 de Feb. de 2019
function main
so=1;
c=-1.25;
% a=1;
Pr=1;
%n=2;
n=input('n=') %%n=-0.3,-0.5
x=[1 1];
x1=fsolve(@solver,x);
function F=solver(x)
[t,u]=ode45(@equation,[0,20],[so c 1 x]);
s=length(t);
F=[u(s,2)-1; u(s,4)];
end
figure(1)
plot(t,u(:,4),'b-');
hold on
function dy=equation(t,y)
dy=zeros(5,1);
dy(1)=y(2);
dy(2)=y(3);
dy(3)=y(2)^2-y(1)*y(3)-1;
dy(4)=y(5);
dy(5)=Pr*(n*y(2)*y(4)-y(1)*y(5));
end
end
ERROR is
code is not running
  21 comentarios
Jan
Jan el 27 de Feb. de 2019
@MINATI: This is a link to asking a new question, but not to your new question.

Iniciar sesión para comentar.

Respuestas (0)

Categorías

Más información sobre Get Started with MuPAD en Help Center y File Exchange.

Community Treasure Hunt

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

Start Hunting!

Translated by