323.sagews

#http://www1.uprh.edu/rbaretti 
#http://www1.uprh.edu/rbaretti/Methodsoftheoreticalphysics.htm 
# www.geocities.com/serienumerica4
# www.geocities.com/serienumerica
f0,f1,f2=var('f0,f1,f2')
solve([.5*f0+.135*f1 +.916*f2==1 ,.5*f0+.135*f1+.00916*f2==.368, .5*f0+.0498*f1+.00124*f2==.111 ],f0,f1,f2)

# A(i,j) 0.500E+00    0.368E+00    0.677E-01

# 0.500E+00    0.135E+00    0.916E-02
 
# 0.500E+00    0.498E-01    0.124E-02

#  bi  0.100E+01    0.250E+00    0.111E+00

R=RealField(20)
f0,f1,f2=var('f0,f1,f2')
solve([.5*f0+.368*f1 +.0677*f2==1 ,.5*f0+.135*f1+.00916*f2==.250, .5*f0+.0498*f1+.00124*f2==.111 ],f0,f1,f2)

[[f0 == 10014991/78556200, f1 == 549265/785562, f2 == 3939125/392781]]

R=RealField(20)
R(3939125/392781)

10.029

R(287945/97554)

2.9516

R(15800/22671)

0.69693