read('s33odr10Eq'):
read('c8fixed');
readlib(write):

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

Digits:=100; DumpDigits:=20:

# Compute s33odr10b.

print(`Compute s33odr10b`):
read('dc16overdc8');
read('iterndc16');
print(`Initial Guess interval is`);
x0U:=vector([0.122839250325515528, 0.899280987883841476, 
	     1.04810650385086190, 0.875707994968721870, 
	     0.328204094328295282, 0.473238688848100275, 
	     0.157658731844991323, 
	     0.457679325816196358, 0.900326815952593451, 
	     0.809745751866159241E-01, 0.215401861690595914, 
	     0.859866512457313936, 0.140566031107548145, 
	     0.350899642820214464, 0.818535036310124170E-01]);
c80U:=0.278751810751748319; 
x0L:=vector([0.122812754681043204, 0.899263694155489457, 
	     1.04807091277297926, 0.875685303235611223, 
	     0.328258442181786514, 0.473460089871075240, 
	     0.157764264750173883, 
	     0.457698203729908659, 0.900322555520043255, 
	     0.809679729327775072E-01, 0.215443714133805642, 
	     0.859869370554803170, 0.140566921616425006, 
	     0.350962943154270290, 0.818535612591022449E-01]);
c80L:=0.278462697549717475;
xM:=vector(15): 
c8M:=c16Itern(x0L,c80L,x0U,c80U,xM):

c0:=vector(16):
for i from 1 to 7 do c0[i]:=xM[i]: od:
c0[8]:=c8M:
for i from 9 to 16 do c0[i]:=xM[i-1]: od:

# Dump out the solution to s33odr10b

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

#
#  To verify that the computed solution by computing
#  the residual error of determining equations.
#
read('ToCheck');
ResiErr(c0);

#
# Compute s33odr10c.

print(`Compute s33odr10c`):
read('dcxsoverdc8'):
read('iterndcxs');
print(`Initial Guess interval is`);
x0U:=vector([0.123291736829209020, 0.900030788004469895, 
	     1.04926722074548473, 0.876792793531598891, 
	     0.326546386604520344, 0.467632832100765228, 
	     0.154405217617667623, 
	     0.456597760311263368, 0.900783167626624692, 
	     0.808028134989672436e-01, 0.214156202877121421, 
	     0.860063462501298726, 0.140451250566628688, 
	     0.349216608299185549, 0.814879993896481386e-01]);
c80U:=0.287224763121685156;
x0L:=vector([0.123082070524670031, 0.899508023240143895, 
	     1.04851675310963466, 0.876021653795914013, 
	     0.327596778939115707, 0.470932669439457163, 
	     0.156472747296796111, 
	     0.457391313010302725, 0.900423646975197056, 
	     0.809921689719436261E-01, 0.214937946935169710, 
	     0.859877428829806667, 0.140542425353092221, 
	     0.350228336099198512, 0.817979993866679023E-01]);
c80L:=0.281927725439035326;
xM:=vector(15): 
c8M:=cxsItern(x0L,c80L,x0U,c80U,xM):

c0:=vector(16):
for i from 1 to 7 do c0[i]:=xM[i]: od:
c0[8]:=c8M:
for i from 9 to 16 do c0[i]:=xM[i-1]: od:

# Dump out the solution to s33odr10c

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

#
#  To verify that the computed solution by computing
#  the residual error of determining equations.
#
read('ToCheck');
ResiErr(c0);

