tps/chap1/tpapprox.py

   1 from math import *
   2
   3
   4 b= 160
   5 c= 1
   6 delta = 160*160 - 4
   7 sqdelta = sqrt(delta)
   8 s_sqd = str(sqdelta)
   9
  10 x1_approx = (-b-float(s_sqd))/2
  11 x1_exact = (-b-sqdelta)/2
  12
  13 x2_approx = (-b+float(s_sqd))/2
  14 x2_exact = (-b+sqdelta)/2
  15
  16 x2p_exact = c/x1_exact
  17 x2p_approx = c/x1_exact
  18
  19 print x1_exact,x1_approx,x2_exact,x2_approx,x2p_exact,x2p_approx
  20
  21