        function ss=uzero(x)
%
        ss=4.*(x^2-1)*exp(x^2);
