texto: metodos numericos para ecuaciones diferenciales ordinarias

=0.3;c=0.1;m=1;g=9.8;k=1600;n=0.1;u=0.6;N=m*g;v0=0.5;
if j==1,
end
y1=x2;
if j==2,
end
y1=(-k/m)*x1-(n/m)*x2+(u*N)*(1-b*(v0-x2)+c*(v0-x2)^3)/m;
Se realiza tres simulaciones para v0 = 0.5m/s y v0 = 1m/s tenemos las figuras (3.2) y
(3.3) donde el sistema es inestable, para v0 = 1.5m/s el sistema es estable según la figura
(3.4).
posicion (m)
0.05
0.04
0.03
0.02
0.01
0
−0.01
−0.02
−0.03
velocidad v0=0.5 m/s
−0.04
0 1 2 3 4
tiempo (s)
5 6 7 8
Figura 3.2: Caso v0 = 0.5m/s, gráfica de desplazamientos versus tiempo
32