def POTDEF1(z, TheTa2): TheTa = np.sqrt(TheTa2**2 + theta1[l]**2) R = D_d * TheTa MN_Bulge_p = MiyamotoNagaiPotential(amp=M1, a=0, b=b1, normalize=False) Sigma = MN_Bulge_p.dens(R, z) kappa = Sigma / SIGMA_CRIT return (4 / theta2[l]) * TheTa2 * kappa / SIGMA_CRIT**2
def integ(z, TheTa1, TheTa2): TheTa = np.sqrt(TheTa1**2 + TheTa2**2) R = D_d * TheTa MN_Bulge_p = MiyamotoNagaiPotential(amp=MN_M, a=MN_a, b=MN_b, normalize=False) Densidad = MN_Bulge_p.dens(R, z) Kappa = 2 * Densidad return Kappa / (SIGMA_CRIT**2)