#
#  Input:        c0:6-vector, the computed solution. 
#  Output:       the infinity norm of the 8-vector whose components
#                are the value of determining functions evaluated at 
#                the input c0.
#
#                              Ren-Cang Li, June 1, 1996
#                              na.rcli@na-net.ornl.gov

ResiErr:=proc(c0:vector)
    local Fv, tmp, i, j;
    Fv:=vector(8,[]);
    for i from 1 to 8 do
        Fv[i]:=subs(cx[1]=c0[1],cx[2]=c0[2],cx[3]=c0[3],
                    cy[1]=c0[4],cy[2]=c0[5],cy[3]=c0[6], f[i]):
    od:
    print(`Residual error of determining equations`);
    print(norm(Fv));
end:
