# Python, perhaps an Infinite Loop

I am new to python so I am sure this is just a rocky mistake, When I try to run the program It just keeps running and doesn’t solve anything so I believe that I have made an infinite loop but I am not sure how.

So my objective is to determine a “me” so my max-height “h” becomes 3000, I have defined my functions “h” and “v” When it enters the loop I want it to compute h and v for the first 20 seconds and then checks if my max-height “Apogeu” as got to 3000, if not I will increase my “me” and try it again.

I also had to increase my recursion limit because it was saying that it was being passed, I have tried to put prints inside the loop to check if they appeared but nothing has shown, I have been trying to figure out what was the problem but I couldn’t. If you could point me in the right direction it would be a huge help.

``````import sys
import math
import numpy as np

"""
Input
"""

ISP    = 130    #(s)
t_burn = 3.0    #(s)
g0     = 9.81   #(m/s^2)
mf     = 12.0   #(Kg)
D_R    = 0.12   #(m)
C_D    = 0.4
rho    = 1.225
A      = math.pi*D_R**2/4
me     = 0.1

sys.setrecursionlimit(10**4)

"""
Variables
"""
def v(t):

if t==0.0:
return 0

elif t<=t_burn:
return (ISP*g0*np.log(m0/(m0-me*t))-g0*t-C_D*0.5*A*v(t-0.01)**2*
np.log(m0-me*t)/me)

elif t>t_burn:
return (v(t_burn)-g0*t-C_D*0.5*A*v(t-0.01)**2/mf*t)

def h(t):

if t==0.0:
return 0

elif t<=t_burn:
return (ISP*g0/me*((m0-me*t)*np.log((m0-me*t)/m0)+me*t)-g0*0.5*t**2-
C_D*0.5*A*v(t-0.01)**2/me*((t-m0/me)*np.log(m0-me*t)+m0/me-t))

elif t>t_burn:
return (h(t_burn) + v(t_burn)*t -0.5*g0*t**2-
C_D*0.5*A*v(t-0.01)**2*t**2*0.5/mf)

"""
Resolution
"""

Apogeu = 0
while (Apogeu<3000):

m0   = mf + me
t    = 0.0

while (t<=20):

v(t)
h(t)
t=t+0.1

Apogeu=max(h)
print("Apogeu=",Apogeu,"m")
me = me + 0.1
``````