计算 cosx 时出现 Fortran 错误

Question

我被分配了以下问题：

制作一个 Fortran 程序，该程序将能够读取学位 [0-360] 检查有效范围（不是类型），它将能够从以下等式计算和打印 cos(x)，其中 x 是以弧度为单位：

cos(x)=1-x^2/2! + x^4/4!-x^6/6!+x^8/8!-...

作为收敛标准假设10^(-5) 使用两次连续重复之间的绝对误差（我想这意味着做）。

对于! 的计算，应该使用最大可能的整数类型。最后应该在屏幕上打印重复的总数。

所以我的代码是这样的：

    program ex6_pr2
implicit none
!Variables and Constants
integer::i
real*8::fact,fact2 !fact=factorial
real,parameter::pi=3.14159265
double precision::degree,radiants,cosradiants,s,oldcosradiants,difference                         !degree,radiants=angle
print*,'This program reads and calculates an angle`s co-sinus'
print*,'Please input the degrees of the angle'
read*,degree
do while(degree<0 .or. degree>360) !number range
read*,degree
print*,'Error input degree' 
cycle
end do
radiants=(degree*pi/180)
fact=1
fact2=1
s=0
cosradiants=0
!repeat structure
do i=2,200,1
fact=fact*i
fact2=fact2*(i+2)
oldcosradiants=cosradiants
cosradiants=(-(radiants)**i/fact)+(((radiants)**(i+2))/fact2)
difference=cosradiants-oldcosradiants
s=s+cosradiants
if(abs(difference)<1e-5) exit
end do
!Printing results
print*,s+1.

end program

对于 45 度（或 pi/4）等角度，我得到正确的结果，而对于 90 度或 180 等其他角度，我得到错误的结果。

我检查了我认为隐藏错误的阶乘（至少对我而言）。

由于以下错误，我创建了另一个似乎无法运行的代码：FUNCTION 名称，（PROJECT2_EX6~FACT 的结果），在未预期的地方使用，可能缺少“（）”

    program project2_ex6
implicit none
integer(kind=3)::degrees,i,sign
integer::n
double precision::x,err_limit,s_old,s
real,parameter::pi=3.14159265359
print*,'This program calculates the cos(x)'
print*,"Enter the angle's degrees"
read*,degrees
do
if(degrees<0.or.degrees>360) then
  print*,'Degrees must be between 0-360'
  else
   x=pi*degrees/180
    exit
    end if
    end do
    sign=1
    sign=sign*(-1)
    err_limit=1e-5
 n=0
    s=0
    s_old=0
do
  do i=1,n
    end do
  s=(((-1.)**n/(fact(2.*n)))*x**(2.*n))*sign
 s=s+s_old

n=n+1
if(abs(s-s_old)<1e-5) then
  exit
  else
    s_old=s
    cycle
    end if
  end do
  print*,s,i,n
contains
real function fact(i)
double precision::fact
integer::i
if(i>=1) then
  fact=i*fact(i-1)
  else
fact=1
    end if
return
end function
  end program

Answer 1

虽然这是你的功课，但我会在这里帮助你。第一个错误是你需要替换的阶乘

 fact  = 1
 do j = 1,i
    fact = fact*j
 enddo

其次，如果你让你的 do 循环完成这项工作会更容易，所以运行它

do i=4,200,2

并使用

在 do 循环之外预定义 cosradians

cosradiants = 1-radiants**2/2

此外，您还需要考虑可以在循环中使用的不断变化的符号

   sign = sign*(-1)

并在循环之前以sign = 1 开始

在循环中

   cosradiants= cosradiants+sign*radiants**i/fact

如果你已经包含了这些东西，它应该可以工作（至少我的代码可以）