diff --git a/galpy/potential_src/SCFPotential.py b/galpy/potential_src/SCFPotential.py
index <HASH>..<HASH> 100644
--- a/galpy/potential_src/SCFPotential.py
+++ b/galpy/potential_src/SCFPotential.py
@@ -402,7 +402,7 @@ def compute_coeffs(dens, N, L):
             
             phi_nl = -NN* (2*l + 1.)**.5 * (1. + xi)**l * (1. - xi)**(l + 1.)*_C(xi, L, N)[:,:,nu.newaxis] * Legandre
             
-            return dens(R,z, phi) * phi_nl[nu.newaxis, :,:,:]*nu.array([nu.cos(m*phi), -nu.sin(m*phi)])*dV
+            return dens(R,z, phi) * phi_nl[nu.newaxis, :,:,:]*nu.array([nu.cos(m*phi), nu.sin(m*phi)])*dV
             
                
         Acos = nu.zeros((N,L,L), float)
Removed a negative sign from the integrand in compute_coeffs