program Caida_Libre
    implicit none
    double precision, dimension(33) :: v, t
    double precision :: m, g, c, t_final, dt
    integer :: i

    print*, "Fisica computacional (Clase 2:00-3:00)"
    print*, "Alumna: Nalleli Iridian Avila Garcia 1941603"
    print*, "Este programa calcula las velocidades de caída libre de un objeto con v0=0m/s, c=12.5kg/s, en un tiempo de 0 a 64 s, con dt=2 seg."

    ! Parámetros
    m = 68.1  ! masa del objeto en kg
    g = 9.81  ! aceleración debido a la gravedad en m/s^2
    c = 12.5  ! coeficiente de arrastre
    t_final = 64.0  ! tiempo final en segundos
    dt = 2.0  ! tamaño del paso de tiempo en segundos

    ! Inicialización de las variables
    t = [(i-1)*dt, i=1,33]  ! vector de tiempo
    v = 0.0  ! vector de velocidad

    ! Método de Euler
    do i = 2, 33
        v(i) = v(i-1) + (g - c/m * v(i-1)) * dt
    end do

    ! Imprimir cada velocidad
    do i = 1, 33
        print '(A,F6.2,A,F6.2,A)', "En el tiempo t=", t(i), "s, la velocidad es ", v(i), "m/s."
    end do
end program Caida_Libre