program esfera
implicit none
integer :: i
real :: xi , xf, tp, y, x, n, k, nc, dydx
xi= 0.0
xf= 5.0
tp= 0.1
y= 0.0
x= xi
n= 2.0
k= 5.0
nc= ( xf- xi) / tp
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
ICAgICAgICBwcm9ncmFtIGVzZmVyYQogICAgICAgIGltcGxpY2l0IG5vbmUKIAogICAgICAgIGludGVnZXIgOjogaQogICAgICAgIHJlYWwgOjogeGksIHhmLCB0cCwgeSwgeCwgbiwgaywgbmMsIGR5ZHgKIAogICAgICAgIHhpPTAuMAogICAgICAgIHhmPTUuMAogICAgICAgIHRwPTAuMQogICAgICAgIHk9MC4wCiAgICAgICAgeD14aQogICAgICAgIG49Mi4wCiAgICAgICAgaz01LjAKIAogICAgICAgIG5jPSh4Zi14aSkvdHAKIAogCiAgICAgICAgcHJpbnQqLHgseQogCiAgICAgICAgZG8gaT0xLCBuYwogCiAgICAgICAgICAgICAgICBkeWR4PTkuODEtOS44MSprKih5KipuKQogICAgICAgICAgICAgICAgeT15K3RwKmR4ZHkKICAgICAgICAgICAgICAgIHg9eCt0cAogCiAgICAgICAgICAgICAgICBQcmludCAqLCB4LCB5CiAKICAgICAgICBlbmRkbwogCiAgICAgICAgZW5kIHByb2dyYW0gZXNmZXJh