program ex1 implicit double precision (a-h,o-z) x = 10864.d0 y = 18817.d0 res = 9.d0*x*x*x*x - y*y*y*y + 2.d0*y*y write(*,*)'P(10864,18817) = ',res,' (exact value : 1)' x = 1.d0/3.d0 y = 2.d0/3.d0 res = 9.d0*x*x*x*x - y*y*y*y + 2.d0*y*y write(*,*)'P(1/3,2/3) = ',res end