作者Semisphere (所指千歌音之处)
看板Fortran
标题Re: [问题] 可以判断整数吗?
时间Sun Dec 26 17:07:11 2010
program Err
implicit none
integer(4) :: i,j,k,iNum
real(4) :: x,rTemp,dx
!real(8) :: x,rTemp,dx
iNum = 25
x = 1.0
do i=1,iNum
rTemp = sqrt(x)
x = (tan(atan(exp(log(rTemp*rTemp))))) + 1.0
dx = abs(x - float(iNum + 1))
write(*,*) dx
end do
! real(4) , Final dx = 3.8146973E-06
! real(8) , Final dx = 1.776356839400250E-014
stop
end
: 以下的 C语言程式设计，
: double 就是 fortran 的 real*8,
: 经过 25次的 复杂运算，产生的
: 绝对误差如下
: 谨 提供参考
: #if 0
: dx= 5.684342e-014
: 请按任意键继续 . . .
: #endif
: // ----------------------------------------------
: #include <stdio.h>
: #include <math.h>
: #include <process.h>
: // ----------------------------------------------
: template <class T>
: T sqr(T x)
: {
: return(x*x);
: }// end of sqr()
: // ----------------------------------------------
: void main()
: {
: int no, i;
: double x, dx;
: no= 25;
: x= 1.0;
: for (i=1;i<=no;i++) {
: x= (tan(atan(exp(log(sqr(sqrt(x))))))) + 1.0;
: }
: dx= fabs(x - (no + 1));
: printf("\n dx= %.6le\n", dx);
: system("pause");
: }// end of main()
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.115.60.193