<jupyter_start><jupyter_text>Shuey<jupyter_code>import matplotlib.pyplot as plt
import numpy as np
plt.style.use(['science', 'notebook', 'grid'])
#shuey para angulos menores que 30 grados
angulo=0.15
#gas sand
Vp1=2280
nu1=0.42
rho1=2190

#wet sand
Vp2=2217
nu2=0.42
rho2=2150
A=(Vp2*rho2-Vp1*rho1)/(Vp2*rho2+Vp1*rho1) #=Rpp
print('El valor de A es igual a   = '+str("{:.2e}".format(A)))
Vp=(Vp1+Vp2)/2
print('El valor de vp es igual a   = '+str("{:.2e}".format(Vp)))
deltaVp=Vp2-Vp1
print('El valor de delta vp es igual a   = '+str("{:.2e}".format(deltaVp)))
rho=(rho1+rho2)/2
print('El valor de rho es igual a   = '+str("{:.2e}".format(rho)))
deltarho=rho2-rho1
print('El valor de drho es igual a   = '+str("{:.2e}".format(deltarho)))
D=(deltaVp/Vp)/(deltaVp/Vp+deltarho/rho)
print('El valor de D es igual a   = '+str("{:.2e}".format(D)))
nu=(nu1+nu2)/2
print('El valor de nu es igual a   = '+str("{:.2e}".format(nu)))
delthanu=nu2-nu1
print('El valor de d nu es igual a   = '+str("{:.2e}".format(delthanu)))
B=A*(D-2*(1+D)*((1-2*nu)/(1-nu)))+(delthanu/(1-nu)**2)
print('El valor de B es igual a   = '+str("{:.2e}".format(B)))
Rpp=A+B*(np.sin(angulo)**2)
print('El valor de Rpp es igual a = '+str("{:.2e}".format(Rpp)))
print('El valor de B*sin(angulo) es igual a = '+str("{:.2e}".format(B*(np.sin(angulo)**2))))<jupyter_output>El valor de B*sin(angulo) es igual a = 1.46e-04