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

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