Re: [問題] 誰能夠提供一程式如何寫
※ 引述《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
討論串 (同標題文章)
Fortran 近期熱門文章
PTT數位生活區 即時熱門文章