ftest.f Listing
|
program test
implicit real*16 (a-z)
a = 3.q0/4.q0
x = 3.q0**(-40.q0)
apx = a+x
write(*,*) 'Parameters: a = 3/4 and x = 3**-40'
write(*,1001) '.8225263339969959081282058400607250240380E-19'
write(*,1000) 'x= ', x
write(*,1000) '(a+x) - a= ', apx-a
write(*,1000) '4*(a+x) - 3*x - 3= ',
1 4*apx-3*x-3
write(*,1000) '0.2*(a+x) + 0.8*x - 0.15= ',
1 .2q0*apx+0.8q0*x-.15q0
write(*,1001) '.1096701778662661210717332307615448891326E-18'
write(*,1000) 'x/(a+x)= ', x/apx
write(*,1000) '1 - a/(a+x)= ', 1-a/apx
write(*,1000) '1/(a/x+1)= ', 1/(a/x+1)
write(*,1000) '(1-(a+0.1*x)/(a+x))/0.9= ',
1 (1-(a+0.1q0*x)/apx)/0.9q0
write(*,1001) '.9999999999999999998903298221337338789283 '
write(*,1000) 'a/(a+x)= ', a/apx
write(*,1000) '1 - x/(a+x)= ', 1-x/apx
write(*,1000) '1/(1+x/a)= ', 1-x/apx
write(*,1000) '(a+x/10)/(a+x) - 0.1*x/(a+x)= ',
1 (a+x/10)/apx - 0.1q0*x/apx
write(*,1001) '.1096701778662661210777470047181541128353E-18'
write(*,1000) ' log((a+x)/a))= ', log(apx/a)
write(*,1000) ' log(a+x)-log(a)= ', log(apx) - log(a)
write(*,1000) ' log(1+x/a)= ', log(1+x/a)
write(*,1000) '-log(1-x/(a+x))= ',-log(1-x/apx)
write(*,*) 'End ftest'
1000 format(a32,g46.40)
1001 format(32x,a46)
end
|