def func_i_alpha(r, l, lam): ia_p = [] for i in range(len(parm.TIME)): ia = parm.get_i_alpha_at(i) w = theta = parm.get_theta_at(i) ia_p.append(r / l * ia + parm.P * lam / l * w * math.sin(theta) + 1 / l * parm.U_ALPHA) ia_p = np.asarray(ia_p) return ia_p * SCALER
def func_w(j, lam): w = [] for i in range(len(parm.TIME)): ia = parm.get_i_alpha_at(i) ib = parm.get_i_beta_at(i) theta = parm.get_theta_at(i) ci = parm.get_ci_at(i) w.append(-(lam / j * (-ia * math.sin(parm.P) * theta + ib * math.cos(parm.P) * theta) - parm.F / j * theta - ci / j)) w = np.asarray(w) return w * SCALER