def res_edo(c, b):
cond_ini = [-512, 824]
#phi_ini = -np.pi/2.
p = Perfil(cond_ini, c, b)
N= 100 #3*np.int(1e4)
dphi= np.pi/ N
phi=np.linspace(0 , np.pi/ 2., N)
x=np.zeros(N)
z=np.zeros(N)
x[0]=cond_ini[0]
z[0]=cond_ini[1]
for n in range(1,N):
p.avanza_rk4(dphi)
x[n]=p.r_actual[0]
z[n]=p.r_actual[1]
return [x, z]
def res_edo(c, b):
"""
Calcula edo numericamente dados c, b
utiliza la clase perfil
"""
cond_ini = [0.0001, 0.0001]
# phi_ini = -np.pi/2.
p = Perfil(cond_ini, c, b)
N = 100 # 3*np.int(1e4)
dphi = np.pi / N
phi = np.linspace(0.0001, np.pi / 2.0, N)
x = np.zeros(N)
z = np.zeros(N)
x[0] = cond_ini[0]
z[0] = cond_ini[1]
for n in range(1, N):
p.avanza_rk4(dphi)
x[n] = p.r_actual[0]
z[n] = p.r_actual[1]
return [x, z]
