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=1.0
nc=50
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
ICAgICAgICBwcm9ncmFtIGVzZmVyYQogICAgICAgIGltcGxpY2l0IG5vbmUKIAogICAgICAgIGludGVnZXIgOjogaSwgbmMKICAgICAgICByZWFsIDo6IHhpLCB4ZiwgdHAsIHksIHgsIG4sIGssIGR5ZHgKIAogICAgICAgIHhpPTAuMAogICAgICAgIHhmPTUuMAogICAgICAgIHRwPTAuMQogICAgICAgIHk9MC4wCiAgICAgICAgeD14aQogICAgICAgIG49Mi4wCiAgICAgICAgaz0xLjAKIAogICAgICAgIG5jPTUwCiAKIAogICAgICAgIHByaW50Kix4LHkKIAogICAgICAgIGRvIGk9MSxuYwogCiAgICAgICAgICAgICAgICBkeWR4PTkuODEtOS44MSprKih5KipuKQogICAgICAgICAgICAgICAgeT15K3RwKmR5ZHgKICAgICAgICAgICAgICAgIHg9eCt0cAogCiAgICAgICAgICAgICAgICBQcmludCAqLCB4LCB5CiAKICAgICAgICBlbmRkbwogCiAgICAgICAgZW5kIHByb2dyYW0gZXNmZXJh