| 查看: 407 | 回复: 1 | ||
ilovexiaomu金虫 (小有名气)
|
[求助]
二重积分问题
|
|
各位虫友，遇到一个比较复杂含参变量的二重积分，想进行数值编程，但是苦于对mathematica代码的不熟悉，现求助各位尝试修改一下代码，能计算二重积分的mathematica代码 新搜素了相关的数值积分程序，但是运行结果出现了不一样的地方，望各位帮忙解答一下： Clear["Global`*"]; f[x_] := x^2 + y; n = 20; a = 0; b = 1; ya = 2; yb = 3; ssx1 = Sum[(1 + (-1)^i)*f[a + i*(b - a)/n], {i, 1, n - 1}];(*ss1=2y2+2y4+\[Ellipsis]+2yn-2*) ssx2 = Sum[(1 - (-1)^i)*f[a + i*(b - a)/n], {i, 1, n - 1}]; sxy[y_] := N[(f[a] + f + ssx1 + 2 ssx2)*(b - a)/3/n, 2] // FullSimplify; ssy1 = Sum[(1 + (-1)^i)*sxy[ya + i*(yb - ya)/n], {i, 1, n - 1}]; ssy2 = Sum[(1 - (-1)^i)*sxy[ya + i*(yb - ya)/n], {i, 1, n - 1}]; s = N[(sxy[ya] + sxy[yb] + ssy1 + 2 ssy2)*(yb - ya)/3/n, 2] // FullSimplify; Print["s=", s] Integrate[sxy[y], {y, 2, 3}] 输出结果为什么会出现不一样的情况呢？ s=0.33+1.0 y Out[199]= 2.8 还望各位帮忙解答 |
回复此楼» 猜你喜欢
综述论文作为代表作会不会影响评审专家的印象分？
已经有12人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有4人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有29人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有9人回复
-
是这周出结果还是下周出结果？
已经有3人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有9人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有4人回复
-
帮忙看看fileCode
已经有13人回复
-
为什么网上很多人说本周 12号出结果
已经有7人回复
-
读博的好处
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求解一道二重积分，顺便请教些问题。
已经有17人回复
- 三角函数的定积分
已经有16人回复
- 请问一下matlab这个是什么问题
已经有10人回复
- 帮忙求解一个含有绝对值符号的二重积分问题 或者告诉怎么求解
已经有8人回复
- 有关磁等价问题
已经有3人回复
- 谁知道今年数二考试大纲，空间解析几何与向量代数考不考
已经有6人回复
- 求助，15考研，纠结于数学二的高数部分应该看哪几章，
已经有14人回复
- 关于积分区间可加性定理的证明
已经有8人回复
- maple二重积分积不出来，积分结果该怎么化简，急急急！！!
已经有9人回复
- 考研福利！ 2014考研数学必看 张宇各种答疑解惑！
已经有5人回复
- 这个函数积分是多少
已经有15人回复
- 请教MATLAB中下面的提示是什么意思？
已经有4人回复
- 求助一个LMI问题
已经有10人回复
- 二重积分极坐标变换下积分限问题
已经有4人回复
- 华师大版数学分析(第三版)第二型曲面积分 部分内容的混涩
已经有7人回复
- 关于数学符号的问题？请教大虾
已经有16人回复
- matlab中系数是矩阵的积分
已经有5人回复
- matlab 积分计算 算不出结果，求帮助！！！
已经有10人回复
- mathematica请教计算结果问题
已经有4人回复
- 求助定积分 ∫e^(-x^2)dx
已经有8人回复
- 这个复杂式子如何计算出求导结果
已经有5人回复
- 【求助】数学积分
已经有29人回复
- 【交流】大家帮忙看下这道题
已经有7人回复
- 【求助】一个二重积分问题【已解决】
已经有17人回复
ilovexiaomu
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1125.1
- 红花: 1
- 帖子: 234
- 在线: 68.5小时
- 虫号: 443234
- 注册: 2007-10-27
- 专业: 凝聚态物性I:结构、力学和
2楼2015-05-19 11:37:48
100