PTT數位生活區 / Fortran

Re: [心得] fortran 做辛普森積分

看板Fortran作者 (交換關聯)時間16年前 (2008/09/11 11:43), 編輯推噓0(0推 0噓 4→)
留言4則, 2人參與, 最新討論串2/4 (看更多)
※ 引述《sjgau (sjgau)》之銘言： : ! do the simpson integrated : sum= 0.0 : ! for i=1 to (ct-2) step 2 do ... : i= 1 : do while (i .LE. (ct-2)) : sum= sum + y(i) + 4*y(i+1) + y(i+2) : i= i + 2 : end do 這地方有較好一點的寫法: real*8 ll ll = 4.0 sum = 2.0*y(1) i = 2 do i = 2, ct-1 sum = ll*y(i) + sum ll = 6.0 - ll enddo sum = 2.0*y(ct) + sum sum=sum/3.0 ==== 全加起來再除3.0 和2/3 4/3 2/3...差再那，我真的不知道。 反正老code都這樣寫，就學起來吧。 -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 220.132.13.207
→
09/11 13:29, , 1F
全部加完再除應該要比每次都除來的好
09/11 13:29, 1F
→
09/11 13:31, , 2F
因為除法而造成數值準確的喪失只做了一次 這是我個人的想法
09/11 13:31, 2F
→
09/11 13:31, , 3F
打錯 該改為 sum/6.0
09/11 13:31, 3F
→
09/11 13:34, , 4F
沒錯 是 /3.0
09/11 13:34, 4F
‣ 更多 DrStein 的文章
文章代碼(AID): #18o9FhDz (Fortran)
討論串 (同標題文章)
本文引述了以下文章的的內容：
以下文章回應了本文：
完整討論串 (本文為第 2 之 4 篇)：
排序： 最舊先 | 最新先 | 留言數
在新視窗開啟完整討論串 (共4篇)
Fortran 近期熱門文章
更多 近期熱門文章 >>
PTT數位生活區 即時熱門文章
更多 即時熱門文章 >>
‣ 更多 DrStein 的文章
文章代碼(AID): #18o9FhDz (Fortran)