PTT數位生活區 / Fortran

[問題] IMSL積分與陣列

看板Fortran作者 (好天氣)時間13年前 (2012/01/18 17:28), 編輯推噓5(5推 0噓 8→)
留言13則, 3人參與, 最新討論串1/2 (看更多)
我使用IMSL做積分運算 方程式中有些變數需要用到陣列傳入比較方便 但是不知道要怎麼寫... 目前是將FUNCTION裡面的變數用成一般的陣列寫法 但都出現錯誤 The syntax of this substring is invalid. 我是新手查書都找不太到QQ 想請教板上各位高手要怎麼寫才正確 感謝! -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 218.171.62.93
推
01/18 18:51, , 1F
有例子可以看看嗎？
01/18 18:51, 1F
→
01/18 22:37, , 2F
例如方程式F=a*(t**3)對t積分，a是一維陣列(/5,10,15/)
01/18 22:37, 2F
→
01/18 22:37, , 3F
要求得積分後的3個答案，我是用call QDAGS
01/18 22:37, 3F
推
01/19 08:39, , 4F
是否可以把你的程式碼放上來？
01/19 08:39, 4F
推
01/19 10:54, , 5F
你的意思是F並不是傳回一個real值而是一個real的陣列嗎？
01/19 10:54, 5F
→
01/19 12:49, , 6F
對對對，我就是想要傳回一個陣列的答案!
01/19 12:49, 6F
推
01/19 14:25, , 7F
那你是怎麼寫的？qdags裡F就是一維函數F(X)，你塞一個陣列
01/19 14:25, 7F
→
01/19 14:26, , 8F
給它的話它當然就不會接受
01/19 14:26, 8F
→
01/19 15:59, , 9F
我只是單純的寫F=a(i)*(t**3)，也不知道F要怎麼重新宣告
01/19 15:59, 9F
→
01/19 15:59, , 10F
所以一直跑出array bounds exceeded
01/19 15:59, 10F
→
01/19 15:59, , 11F
看樣子不能用QDAGS寫我要的積分方程式了，謝謝!
01/19 15:59, 11F
推
01/19 23:44, , 12F
你可以自己再寫一個用qdags的副程式以達到你的需求
01/19 23:44, 12F
→
01/20 17:37, , 13F
後來用全域變數寫出來了，謝謝t大熱心指教!
01/20 17:37, 13F
‣ 更多 wusudevi 的文章
文章代碼(AID): #1F5f2r-q (Fortran)
討論串 (同標題文章)
以下文章回應了本文：
完整討論串 (本文為第 1 之 2 篇)：
排序： 最舊先 | 最新先 | 留言數
在新視窗開啟完整討論串 (共2篇)
Fortran 近期熱門文章
更多 近期熱門文章 >>
PTT數位生活區 即時熱門文章
更多 即時熱門文章 >>
‣ 更多 wusudevi 的文章
文章代碼(AID): #1F5f2r-q (Fortran)