FUNCTION cbrt(x) ! ! Return the cubic root for positive or negative x IMPLICIT NONE REAL :: x,cbrt cbrt=sign(1.,x)*(abs(x)**(1./3.)) END FUNCTION cbrt