Vamos con ello.
A será la aceleración y v la velocidad, ambas dependen del tiempo. Y la aceleración es la derivada de la velocidad respecto del tiempo.
G la aceleración gravitatoria que 9,8 m/s^2
M es la masa
K la constante de proporcionalidad de la fuerza de resistencia
La suma de fuerzas en el eje Y será la masa por la aceleración en dicho ejem
F = ma = m·dv/dt = mg - kv^2
m·dv/dt = mg - kv^2
[m/(mg - kv^2)] dv = dt
Es una ecuación de variables separables aunque no inmediata.
$[m/(mg-kv^2)]dv =
La separaremos en dos fracciones sencillas de la forma
a/[sqrt(mg) - vsqrt(k)]; b/[sqrt(mg) + vsqrt(k)] que sumadas den la original.
Sumamos ambas fracciones con el algoritmo de la suma de fracciones. El denominador quedará igual, luego también debe hacerlo el numerador y cumplirse:
a[sqrt(mg) + vsqrt(k)] + b[sqrt(mg) - vsqrt(k)] = m
El termino en v no existe, luego
a·v·sqrt(k) - b·v·sqrt(k) = 0 ==>
a=b
a·sqrt(mg)+a·sqrt(mg) = m
a = m/[2sqrt(mg)] = (1/2)sqrt(m/g)
Y las integrales de las fracciones son estas dos:
-(1/2)sqrt(m/g)ln[sqrt(mg) - vsqrt(k)] / sqrt(k)
(1/2)sqrt(m/g)ln[sqrt(mg) + vsqrt(k)] / sqrt(k)
Sumadas y aplicadas todas las propiedades de los logaritmos tendremos
ln {[sqrt(mg) + vsqrt(k)] / [sqrt(mg) - vsqrt(k)]}^[(1/2)(sqrt(m/gk))]
Esto es el lado izquierdo, en el derecho habremos integrado y obtenido t+C.
ln ({[sqrt(mg) + vsqrt(k)] / [sqrt(mg) - vsqrt(k)]} ^ [(1/2)(sqrt(m/gk))]) = t + C
Se simplifica un poco si multiplicamos y dídimos por sqrt(k) en la base y por sqrt(m) en el exponente
ln ({[sqrt(mgk)+kv]/[sqrt(mgk)-kv]}^{m/[2sqrt(mgk)]}) = t + C
Ahora tenemos que imponer que para t=0 se tiene v=0
ln(1) = C
C = 0
Con eso solución particular es:
ln ({[sqrt(mgk)+kv]/[sqrt(mgk)-kv]}^{m/[2sqrt(mgk)]}) = t
Yo creo que incluso podemos despejar v
{[sqrt(mgk)+kv]/[sqrt(mgk)-kv]}^{m/[2sqrt(mgk)]} = e^t
[sqrt(mgk)+kv]/[sqrt(mgk)-kv] = e^[2t·sqrt(mgk)/m]
sqrt(mgk)+kv = [sqrt(mgk)-kv] e^[2t·sqrt(mgk)/m]
kv{1+e^[2t·sqrt(mgk)/m]} = sqrt(mgk) {-1 + e^[2t·sqrt(mgk)/m]}
v(t) = (1/k){-1+e^[2t·sqrt(mgk)/m]} / {1+e^[2t·sqrt(mgk)/m]}
Para simplificar al máximo, llamemos d = 2sqrt(mgk)/m
v(t) =(1/k) (-1+e^dt)/(1+e^dt)
el lim t-->+oo v(t) es (1/k) porque el 1 y -1 se hacen despreciables ante la magnitud de e^(+oo) quedando (1/k) ·e^(+oo) / e^(+oo) = 1/k
Luego la velocidad terminal es 1/k. Donde k era la constante de la fuerza de resistencia, es decir, la constante k que hacia que Fr = kv^2.
Y eso es todo, espero que te sirva y lo hayas entendido. Yo creo que está bien hecho, pero un poco lioso sí a sido.