FUNCTION cbrt(x) IMPLICIT NONE REAL x,cbrt cbrt=sign(1.,x)*(abs(x)**(1./3.)) RETURN END