作者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