program esfera
        implicit none
 
        integer :: i, nc
        real :: xi, xf, tp, y, x, n, k, dydx
 
        xi=0.0
        xf=5.0
        tp=0.1
        y=0.0
        x=xi
        n=2.0
        k=-5.0
 
        nc=20
 
 
        print*,x,y
 
        do i=1,nc
 
                dydx=9.81-9.81*k*(y**n)
                y=y+tp*dydx
                x=x+tp
 
                Print *, x, y
 
        enddo
 
        end program esfera