// Programa para hallar numericamente la posicion en funcion del tiempo
// para un surfeador de nubes
// Secuencia Fuerzas de rozamiento - Actividad 3

m=75;
u=0.2;

// Condiciones iniciales
// En este caso tomaremos un sistema de referencia con origen en el piso 
// y una posicion inicial de 3000 m
// y con sentido negativo hacia abajo

// posicion
x(1)=3000;
// velocidad
v(1)=0;
// aceleracion, F/m = - g - (u/m) * v
a(1)= - 10 - (u/m) * v(1)* abs(v(1));

tiempo(1)=0;

for i=1:10000
  tiempo(i+1)= tiempo(i) + (20/10000);
  x(i+1)= x(i) + v(i) * (20/10000);
  v(i+1)= v(i) + a(i) * (20/10000);
  a(i+1)= - (10) - (u/m) * v(i+1)*abs((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")