作者wusudevi (好天氣)
看板Fortran
標題[問題] IMSL積分與陣列
時間Wed Jan 18 17:28:18 2012
我使用IMSL做積分運算
方程式中有些變數需要用到陣列傳入比較方便
但是不知道要怎麼寫...
目前是將FUNCTION裡面的變數用成一般的陣列寫法
但都出現錯誤 The syntax of this substring is invalid.
我是新手查書都找不太到QQ
想請教板上各位高手要怎麼寫才正確
感謝!
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 218.171.62.93
1F：推 terryys:有例子可以看看嗎？ 01/18 18:51
2F：→ wusudevi:例如方程式F=a*(t**3)對t積分，a是一維陣列(/5,10,15/) 01/18 22:37
3F：→ wusudevi:要求得積分後的3個答案，我是用call QDAGS 01/18 22:37
4F：推 laa7352:是否可以把你的程式碼放上來？ 01/19 08:39
5F：推 terryys:你的意思是F並不是傳回一個real值而是一個real的陣列嗎？ 01/19 10:54
6F：→ wusudevi:對對對，我就是想要傳回一個陣列的答案! 01/19 12:49
7F：推 terryys:那你是怎麼寫的？qdags裡F就是一維函數F(X)，你塞一個陣列 01/19 14:25
8F：→ terryys:給它的話它當然就不會接受 01/19 14:26
9F：→ wusudevi:我只是單純的寫F=a(i)*(t**3)，也不知道F要怎麼重新宣告 01/19 15:59
10F：→ wusudevi:所以一直跑出array bounds exceeded 01/19 15:59
11F：→ wusudevi:看樣子不能用QDAGS寫我要的積分方程式了，謝謝! 01/19 15:59
12F：推 terryys:你可以自己再寫一個用qdags的副程式以達到你的需求 01/19 23:44
13F：→ wusudevi:後來用全域變數寫出來了，謝謝t大熱心指教! 01/20 17:37