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
ICAgICAgICBwcm9ncmFtIGVzZmVyYQogICAgICAgIGltcGxpY2l0IG5vbmUKIAogICAgICAgIGludGVnZXIgOjogaSwgbmMKICAgICAgICByZWFsIDo6IHhpLCB4ZiwgdHAsIHksIHgsIG4sIGssIGR5ZHgKIAogICAgICAgIHhpPTAuMAogICAgICAgIHhmPTUuMAogICAgICAgIHRwPTAuMQogICAgICAgIHk9MC4wCiAgICAgICAgeD14aQogICAgICAgIG49Mi4wCiAgICAgICAgaz0tNS4wCiAKICAgICAgICBuYz0yMAogCiAKICAgICAgICBwcmludCoseCx5CiAKICAgICAgICBkbyBpPTEsbmMKIAogICAgICAgICAgICAgICAgZHlkeD05LjgxLTkuODEqayooeSoqbikKICAgICAgICAgICAgICAgIHk9eSt0cCpkeWR4CiAgICAgICAgICAgICAgICB4PXgrdHAKIAogICAgICAgICAgICAgICAgUHJpbnQgKiwgeCwgeQogCiAgICAgICAgZW5kZG8KIAogICAgICAgIGVuZCBwcm9ncmFtIGVzZmVyYQ==