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

ResiErr:=proc(xyzw0:vector)
    local Fv, tmp, i, j;
    Fv:=vector(3,[]);
    for i from 1 to 3 do
        Fv[i]:=subs(x=xyzw0[1],y=xyzw0[2],z=xyzw0[3],w=xyzw0[4],f[i]):
    od: 
    print(`Residual error of determining equations`);
    print(norm(Fv));
end:
