求二重变上限积分

图中后面那个是N2，更正N2(15，4.5)
求QD的值！！
请给出程序或命令的步骤！
谢谢！！

首先，你把两个文件的后缀名改为.m;其次，你把两个文件都放在当前工作目录下；

最后，运行main.m就OK了

楼主，根据经验，因为正态分布的累计分布在数学上是积不出来，所以不能求出精确解，只能搞个数值解了

约1023.3941，明早给你发程序

mklyzl 发表于 2014-3-25 09:46
首先，你把两个文件的后缀名改为.m;其次，你把两个文件都放在当前工作目录下；

最后，运行main.m就OK了
...

谢谢你，可不可以帮我解释一下，是否你把方程化简了。特别是1023.393这样的数字是怎么来的？

${{Q}_{D}}\approx \text{1023}\text{.394086}$

Mathematica代码如下（简单的二分法）：

1. ff[x_] :=
   
2. SetPrecision[10.4, 10] (CDF[NormalDistribution[1000., Sqrt[300]], x] -
   
3.      CDF[NormalDistribution[1000., Sqrt[300]], 0]) +
   
4.   6.6 NIntegrate[
   
5.     PDF[NormalDistribution[15., Sqrt[4.5]], y] PDF[
   
6.       NormalDistribution[1000., Sqrt[300]], z], {z, 0, x}, {y, 0,
   
7.      x - z}] -
   
8.   14.;(*定义函数*)
   
9. a = 1000.;(*起始左端点*)
   
10. b = 2000.;(*起始右端点*)
    
11. ee = (b - a)/2;
    
12. (*误差*)While[ee > 10^(-7), c = (a + b)/2;
    
13. sa = Sign[ff[a]];
    
14. sb = Sign[ff[b]];
    
15. sc = Sign[ff[c]];
    
16. If[sc == 0, Break, If[sa sc < 0, b = c, a = c]];
    
17. ee = b - a;]
    
18. c(*结果*)

复制代码

mklyzl 发表于 2014-3-25 09:46
首先，你把两个文件的后缀名改为.m;其次，你把两个文件都放在当前工作目录下；

最后，运行main.m就OK了
...

请把解释稍微写一下，放在一个小文件里，设置100论坛币。

客初 发表于 2014-3-25 10:28
请把解释稍微写一下，放在一个小文件里，设置100论坛币。

响应客初的要求，上传个附件
Mathematica程序及详细说明（复制－粘贴－运行）.txt (848 Bytes, 需要: 100 个论坛币)

2014-3-25 10:54:06 上传

需要: 100 个论坛币  [购买]

客初 发表于 2014-3-25 10:28
请把解释稍微写一下，放在一个小文件里，设置100论坛币。

m文件里面的那个quad是数值积分的命令，myfun1是被积函数。
由于我是做统计模型的，非线性方程求解不熟，所以就只能是这种弱爆的方法了。
另外，楼主可采用楼下的用mathmeticas写的程序。
没有完美解决楼主问题，很抱歉！

凡星有梦 发表于 2014-3-25 10:54
响应客初的要求，上传个附件

感谢！已下载。