Usa matlab para dibujar ecuaciones diferenciales ordinarias

Primero, cree una función personalizada para ecuaciones diferenciales

function dy=odefun(t, y, a);

dy(1)=y(2 ) ;

dy(2)=(1/a)*sqrt(1 y(1)^2);

dy=dy(:);

En segundo lugar, use la función ode45() para encontrar su solución numérica y luego use la función plot() para dibujar el diagrama de relación x——y(x), x——dy/dx

a=1;

y0=;

[t, y]=ode45(@(t,y)odefun(t,y,a),[0 20],y0 )

plot(t, y(:, 1), '*', t, y(:, 2), '.-'), cuadrícula activada, espera

legend('y( x)', 'dy/dx', 'Ubicación', 'noroeste') )

xlabel('x'), ylabel('y(x), dy/dx' )

Ejecutar resultados