// Programa para hallar numericamente la posicion en funcion del tiempo
// para una esfera de 2 g que cae en un aceite con un factor de rozamiento de 390 g/s
// Secuencia Fuerzas de rozamiento - Actividad 2

m=0.002;
b=0.39;

// Condiciones iniciales
// considerando un sistema de referencia con origen en la posicion desde la que se suelta la esfera
// y con sentido negativo hacia abajo
// posicion
x(1)=0;
// velocidad
v(1)=0;
// aceleracion, F/m = - g - (b/m) * v
a(1)= - 10 - (b/m) * v(1);

tiempo(1)=0;

for i=1:10000
  tiempo(i+1)= tiempo(i) + (4/1000000);
  x(i+1)= x(i) + v(i) * (4/1000000);
  v(i+1)= v(i) + a(i) * (4/1000000);
  a(i+1)= - 10 - (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("aceleracion en funcion del tiempo")