#
#  Input:        c0:7-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(cx10, c0:vector)
    local Fv, tmp, i, j;
    Fv:=vector(8,[]);
    for i from 1 to 8 do
        Fv[i]:=subs(cx[1]=cx10,  cx[2]=c0[1], cx[3]=c0[2], cx[4]=c0[3],
                    cy[1]=c0[4], cy[2]=c0[5], cy[3]=c0[6], cy[4]=c0[7],
                    f[i]);
    od:
    print(`Residual error of determining equations`);
    print(norm(Fv));
end:
