#
#  Input:        c0:16-vector, the computed solution. 
#  Output:       the infinity norm of the 15-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(15,[]);
    for i from 1 to 15 do
        Fv[i]:=subs(c[ 1]=c0[ 1],c[ 2]=c0[ 2],c[ 3]=c0[ 3],
		    c[ 4]=c0[ 4],c[ 5]=c0[ 5],c[ 6]=c0[ 6],
		    c[ 7]=c0[ 7],c[ 8]=c0[ 8],c[ 9]=c0[ 9],
		    c[10]=c0[10],c[11]=c0[11],c[12]=c0[12],
		    c[13]=c0[13],c[14]=c0[14],c[15]=c0[15],
		    c[16]=c0[16],f[i]):
    od: 
    print(`Residual error of determining equations`);
    norm(Fv);
end:
