24小时热门版块排行榜    

查看: 3332  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

封晨洁

金虫 (初入文坛)

[求助] 求三重积分的详细程序

看过徐士良的fortran常用算法程序集,我的被积函数是一个矢量函数f1(x,y,z)i+f2(x,y,z)j+f3(x,y,z)k,具体表达式知道且表达式复杂(有e指数,分数,根号等)其次在输出结果的时候,要输出三重积分后分别的满足f1,f2,f3的数据,我想得到数据在tecplot里做图像。求大虾能帮忙,指导程序该怎么写
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

封晨洁

金虫 (初入文坛)

引用回帖:
4楼: Originally posted by wangww2011 at 2014-03-20 18:19:15
积分区域从无限大并不是个问题 你可以用高斯 http://en.wikipedia.org/wiki/Gaussian_quadrature 但其实这个也不是很重要 你就用普通的梯形法都行的 就是取一个很大的边界就行了 如1e10之类的

不好意思 你有困难 ...

帮帮我看看吧。http://d.kankan3d.com/file/data/bcs/2014/0320/w367h3049570_1395319100_586.jpg。。。http://d.kankan3d.com/file/data/bcs/2014/0320/w236h3049570_1395319115_122.jpg
求三重积分的详细程序
Z_[~3M6WJL2XIND8RRAGB]3.jpg


求三重积分的详细程序-1
R6_4NH1B_0_ERJEPPM2]{(1.jpg

5楼2014-03-20 20:38:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

wangww2011

木虫 (著名写手)

数值积分还是很成熟的 如果是单个的 上这里最简单
https://www.wolframalpha.com/input/?i=int_0^1+e^x%2Fsqrt%28x%29+dx

如果是根据不同的边界算出不同的积分值 并且最好导出到一定格式的数据 那就侧重最后文件的格式就行 中间的可以自己写个函数 或者调用 一些软件自带的函数 如用matlab来写 最为方便 当然用Fortran 也行的
2楼2014-03-20 08:23:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

封晨洁

金虫 (初入文坛)

引用回帖:
2楼: Originally posted by wangww2011 at 2014-03-20 08:23:46
数值积分还是很成熟的 如果是单个的 上这里最简单
https://www.wolframalpha.com/input/?i=int_0^1+e^x%2Fsqrt%28x%29+dx

如果是根据不同的边界算出不同的积分值 并且最好导出到一定格式的数据 那就侧重最后文 ...

我的积分区域是从有界到无界(这里用高斯-拉盖尔?)。被积函数是矢量,function下(可以编写多个函数表达式吗?)有三个函数表达,输出结果也是分别满足三个函数的输出值f1,f2,f3.在fortran里编写,我有困难,求帮助
3楼2014-03-20 11:30:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangww2011

木虫 (著名写手)

引用回帖:
3楼: Originally posted by 封晨洁 at 2014-03-20 11:30:46
我的积分区域是从有界到无界(这里用高斯-拉盖尔?)。被积函数是矢量,function下(可以编写多个函数表达式吗?)有三个函数表达,输出结果也是分别满足三个函数的输出值f1,f2,f3.在fortran里编写,我有困难,求帮 ...

积分区域从无限大并不是个问题 你可以用高斯 http://en.wikipedia.org/wiki/Gaussian_quadrature 但其实这个也不是很重要 你就用普通的梯形法都行的 就是取一个很大的边界就行了 如1e10之类的

不好意思 你有困难 但也没有说清楚具体的困难在哪里?个人觉得你这种问题 程序最复杂的地方是 如何把数据保存成tecplot可以读取的文件了 因为积分的话 如果你用其他程序 直接调用就行了 如http://docs.scipy.org/doc/scipy/reference/tutorial/integrate.html 所以积分本身你根本不用花时间的
4楼2014-03-20 18:19:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:52 by Z9YWQ5EAO3qp
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:45 by Z9YWQ5EAO3qp
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:33 by Z9YWQ5EAO3qp
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:14 by Z9YWQ5EAO3qp
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 13:35 by C3vABopYbZJS
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +10 zju2000 2026-09-26 13/650 2026-09-30 10:26 by yexuqing
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:37 by etmYJ6d2rquH
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 IDs3scOF0tjC 2026-09-28 7/350 2026-09-29 14:07 by etmYJ6d2rquH
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 09:59 by JzYBbHIXWSrW
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
信息提示
请填处理意见