计算的sin（x）使用Fortran语言

 integer n 
    real term , sum , deg 
    write(*,*) 'Enter Degree' 
    read(*,*) deg 
    deg = deg * 3.14 /180 
    n = 3 
    term = deg 
    sum = 0 
2  if (abs(term) .gt. 0.000001) then !<<<<<<<<<<< THIS CONDITION 
     goto 1 
    else 
     goto 3 
    endif 
1  sum = sum + term 
    write(*,*) 'Your', n - 2, ' Term is ' , term 
    term = term *((deg ** 2)/ (n *(n - 1))) * (-1) 
    n = n + 2 
    goto 2 
3  write(*,*) ' YOur final sum ' , sum 
    pause 
    end

我发现这个程序计算的sin（x）很明显的罪（x）由用户进入了我的价值没有得到的条件整点（abs（term）.gt。0.000001）这是否意味着计算机不能比这更精确。纠正我，如果我错了计算的sin（x）使用Fortran语言

来源

2013-04-10 Rohit Sthapit

该程序使用默认的实际变量。他们通常可以精确到约。 6位数字。你可以使用所谓的双精度，这可以允许更多。下面你看到15位数的例子。

integer,parameter :: dp = selected_real_kind(p=15,r=200) 
real(dp) :: term , sum , deg 

deg = deg * 3.14_dp /180

等等...

参见：

http://gcc.gnu.org/onlinedocs/gfortran/SELECTED_005fREAL_005fKIND.html

http://gcc.gnu.org/onlinedocs/gfortran/ISO_005fFORTRAN_005fENV.html（尤其是real64）

在旧的程序，你也可以看到

double precision x

其是过时的，或

real*8 x

这是非标准的。

来源

2013-04-10 08:41:26

条件if (abs(term) .gt. 0.000001)是测试该术语非零的一种方法。使用整数，您只需使用if (term .ne. 0)，但对于实数，它可能不会在内部表示为相同的零。 if (abs(term) .gt. 0.000001)过滤实数精度内非零的数字。

来源

2013-04-26 06:41:58 SethMMorton

计算的sin（x）使用Fortran语言

回答

相关问题