read('s33odr10Eq'):
read('c3fixed'):
readlib(write):

#                                        Ren-Cang Li, June 1, 1996
#                                        na.rcli@na-net.ornl.gov

Digits:=100: DumpDigits:=20:

x0:=vector([0.705265711652884181E-01, 0.944488520327729231, 
	                             0.933093349000206129, 
	    0.306824155844901414, 0.543507631609470465, 
	    0.121449884910852221, 0.363646387537339133, 
	    0.411052830136318548, 0.954768545706830563, 
	    0.224628372697403024E-01, 0.192020631368324729, 
	    0.907948413610440874, 0.107956257679331730, 
	    0.345694008299979427, 0.224639999570846738E-01]);
c3:=1;
err:=c3fixed(x0,c3):
print(`Residual error of determining equations`);
print(err):

c0:=vector(16):
c0[1]:=x0[1]: c0[2]:=x0[2]: c0[3]:=c3:
for i from 4 to 16 do: c0[i]:=x0[i-1]: od:

# Dump out the solution to s33odr10a

appendto(s33odr10a):
writeln():
read('WriteItOut'):
close():
