function cg = sf_cg(J)
//SF_CG([J1 J2 J3; M1 M2 M3])
//
//   Clebsch-Gordan
//
//   GCG, 16-Mar-2003

fac    = (-1)^(J(1,1)-J(1,2)+J(2,3))*sqrt(2*J(1,3)+1)
J(2,3) = -J(2,3)
cg     = fac * sf_threej(J)

endfunction
