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]