program esfera
        implicit none

	integer :: i
	real :: xi, xf, tp, y, x, n, k, nc, dydx

	xi=0.
	xf=5.
	tp=.1
	y=0.
	x=xi
        n=2.
	k=5.

	nc=(xf-xi)/tp  !intervalo


        print*,x,y

        do i=1,nc

        	dydx=9.81-9.81*k*(y**n)
		y=y+tp*dxdy
		x=x+tp

        	Print *, x, y

        enddo

	end program esfera