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)