double precision y
       real x
       write(*,*) "x?"
       read(*,*) x

       y=3.31132*(X**0)-0.0241461*(x**1)+0.000134561*(x**2) 
     .-(4.00458e-07)*(x**3)+(6.63974e-10)*(x**4)-(5.79031e-13)*x**5 
     .+(2.07175e-16)*x**6 


      write(*,*) 'y=',y
      stop
      end



c       y=(3.31*(x**0))-(0.0241*(x**1))+(0.000135*(x**2))
c     .-((4e-7)*(x**3))+((6.64e-10)*(x**4))+((5.79e-13)*(x**5))
c     .+((2.07e-16)*(x**6))
      
cY = 3.31*X^0  -  0.0241*X^1  +  0.000135*X^2 
c     -  4e-07*X^3  +  6.64e-10*X^4  -  5.79e-13*X^5 
c     +  2.07e-16*X^6