ejercicio 3

ejercicio 3 mov del sist masa-resorte
function ejercicio3
[t1, xx1]=ode45(@muelle, [0,40], [1,0],[ ], 5);
plot(t1,xx1(:,1))

hold on

[t2, xx2]=ode45(@muelle, [0,40], [1,0],[ ], 40);
plot(t2,xx2(:,1))

[t3, xx3]=ode45(@muelle, [0,40], [1,0],[ ], 200);
plot(t3,xx3(:,1))
xlabel('tiempo');
ylabel('desplazamiento');
title('desplaz en funcion del tiempo para cada coef de amortiguamiento');
legend('c=5','c=40','c=200');
end

% funcion muelle
function x=muelle(t,x,c)
x=[x(2); (-c*x(2)-20*x(1))/20];
end