| 4 | 1/1 | 返回列表 |
| 查看: 572 | 回復(fù): 3 | ||
luohubin木蟲(chóng) (正式寫(xiě)手)
|
[求助]
關(guān)于一個(gè)積分算法的問(wèn)題(有代碼)
|
|
請(qǐng)教大家一個(gè)散點(diǎn)積分的問(wèn)題。我將fortran代碼貼在下面。 說(shuō)明一下參數(shù):f是一維矩陣,存放n個(gè)散點(diǎn)。散點(diǎn)間等距,為d。fi是積分值。 問(wèn)題是,我看不出這個(gè)積分是基于什么方法的,在數(shù)值分析書(shū)上沒(méi)找到對(duì)應(yīng)的。希望大家給我點(diǎn)建議。謝謝! SUBROUTINE gensim(f,d,n,fi) IMPLICIT NONE REAL(KIND=8), DIMENSION(*) :: f REAL(KIND=8) :: d, fi, ff, corr INTEGER :: n, n1, j, n2 C IF(n.GE.3) THEN IF(MOD(n,2).EQ.0) THEN n1=n-1 ff=-f(n-2)+8.d0*f(n-1)+5.d0*f(n) corr=ff*d/12.d0 ELSE n1=n corr=0.d0 ENDIF n2=n1-1 ff=0.d0 DO 20 j=2,n2,2 20 ff=ff+(f(j-1)+4.d0*f(j)+f(j+1)) fi=d*ff/3.d0+corr ELSEIF(n.EQ.2) THEN ff=5.d0*f(1)+8.d0*f(2)-f(3) fi=d*ff/12.d0 ELSEIF(n.LE.1) THEN fi=0.d0 ENDIF C RETURN END |

木蟲(chóng) (正式寫(xiě)手)

木蟲(chóng) (正式寫(xiě)手)

木蟲(chóng) (正式寫(xiě)手)

| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +9 | Ymlll 2026-03-18 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專(zhuān)業(yè)(081702)315分求調(diào)劑 +11 | yangfz 2026-03-17 | 11/550 |
|
|
[考研] 一志愿西安交通大學(xué)材料工程專(zhuān)業(yè) 282分求調(diào)劑 +5 | 楓橋ZL 2026-03-18 | 7/350 |
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
[考研] 化學(xué)求調(diào)劑 +3 | 臨澤境llllll 2026-03-17 | 4/200 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] 一志愿天大材料與化工(085600)總分338 +5 | 蔡大美女 2026-03-13 | 5/250 |
|
|
[教師之家] 焦慮 +9 | 水冰月月野兔 2026-03-13 | 13/650 |
|
|
[考研] 材料工程專(zhuān)碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
[考博] 環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
[考研] 265求調(diào)劑 +3 | 梁梁校校 2026-03-17 | 3/150 |
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|