/*
** NL4.E
*/

library nlsys;

proc f(x);
    local f1,f2,f3,f4,f5;
    f1 = sin(x[1])*cos(x[2]) + 4*cos(x[4]);
    f2 = cos(x[4]) - sin(x[5])*x[2]^2;
    f3 = 1/x[5] - x[3]*x[4];
    f4 = sin(x[1]^2) - cos(x[5])*sin(x[2]);
    f5 = x[1]^2 + x[2]*x[3] + x[4]^3 -  x[5];
    retp(f1|f2|f3|f4|f5);
endp;


output file = nl4.out reset;
nlset;
let x0 = 1 2 3 4 5;
_nlchpf = 1;
_nlalgr = 2;
{ x,fvp,jc,tcode }  = nlprt(nlsys(&f,x0));

output off;