// Programa para hallar numericamente la posicion en funcion del tiempo
// para un cuerpo unido al extremo de un resorte considerando el rozamiento con el aire
// Secuencia Fuerzas de rozamiento - Actividad 2

// constante el�stica del resorte
k=7.0;
m=2.3;
b=0.35;

// Condiciones iniciales
// posicion
x(1)=0.3;
// velocidad
v(1)=0;
// aceleracion, F/m = - k x/m - (b/m) v
a(1)= - (k/m) * x(1) - (b/m) * v(1);

tiempo(1)=0;

for i=1:100000
  tiempo(i+1)= tiempo(i) + (6/10000);
  x(i+1)= x(i) + v(i) * (6/10000);
  v(i+1)= v(i) + a(i) * (6/10000);
  a(i+1)= - (k/m) * x(i+1) - (b/m) * v(i+1);
end

plot(tiempo, x)
xtitle("posicion en funcion del tiempo")
scf(1)
plot(tiempo, v)
xtitle("velocidad en funcion del tiempo")
scf(2)
plot(tiempo, a)
xtitle("aceleraci�n en funcion del tiempo")