Added initial z0 estimate, switched squared parameters to exp()