作者tss1030 (josh)
看板Fortran
标题Re: [问题] 谁能够提供一程式如何写
时间Fri May 1 14:51:51 2009
※ 引述《tss1030 (josh)》之铭言：
: 最近老师出了一个题目
: 要用fortran写cosx=x的数值分析
: 小弟fortran算是初学者
: 烦请高手能够提供我程式或指导如何编写
: 另外fortran有办法能够算到50位吗?
不好意思再请问一下
我下面这样写是跟我题目是一样意思吗
program main
real*8 tcos, x
write (*,*) 'Please enter x:'
read (*,*) x
write (*,*) 'x =', x
write (*,*) 'cos(x) by calculation =', tcos(x)
write (*,*) 'cos(x) by function call =', cos(x)
end program main
real*8 function tcos(x)
integer*2 n, n2
real*8 x, fac, term
tcos = 1d0
fac=1d0
n = 0
1 n = n+1
n2 = n*2
fac = fac * (n2-1) * n2
term = x**n2 / fac
tcos = tcos + (-1)**n*term
if (term .ge. 1d-10) go to 1
return
end function tcos
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 59.115.160.219