Programación en MATLAB para Runge-Kutta para resolver ecuaciones diferenciales de segundo orden

MATLAB puede resolver x''+0.7x'+0.8x'|x'|+25.6x-25.6x?=0 ecuaciones diferenciales de segundo orden siguiendo los siguientes pasos:

1. Cree una función personalizada func()

función

f

=

func(t,x)

%x''+0.7x'+0.8x'|x'|+25.6x-25.6x?=0

f(1)=x(2);

f(2)=25.6*x(1)^3-25.6*x(1)-0.8*x(2)*abs(x(2))-0.7*x(2);

f=f(:);

2. Establecer la función del algoritmo Runge-Kutta runge_kutta()

Formato de llamada: [t,x]

=

runge_kutta(@(t,x)func(t,x),x0,h,a,b);

3. Datos, dibuja las gráficas de x(t), x′(t).

trama(x(:,1),x(:,2))