[ Foro de Python ]

que hice mal?

G    =   6.67e-11 # Constante de gravitacion universal
rho  =   2000 # Densidad
prof =   3 # Profundidad
Radio =  2
x     = np.linspace (-10, 10, 101) # Este comando produce un arreglo
                                  # de 101 números del -10 al 10.

def g_cilindro(x, z, rho, R):
   g_cil = ((2*(np.pi)*G*rho*(4)/prof))
   return g_cil

gz = g_cilindro(x, prof, rho, Radio)
fig, ax = plt.subplots(1,1, figsize=(12,5))
ax.plot(x,gz, linewidth = 3)
ax.set_xlabel('Distancia [m]', fontsize=16)
ax.set_ylabel(r'$g_{z}[mGal]$', fontsize=16)
ax.grid()