Monte Carlo 仿真 ，求救

杨学硕

用Monte Carlo方法计算如下积分

我用的方法如下
a=-1;b=1;c=g(-1/2);d=g(1);s_0=(b-a)*(d-c);n=10^4;x=runif(n);y=runif(n);mu_n=sum(y<=f(x));J=mu_n/n;J_0=s_0*J+c*(b-a);
其中又定义了两个函数function [ y ] = f( x )(g(a+(b-a)*y)-c)/(d-c);end

function [ y ] = g( x )
1/(x.^2+x+1);
end



错误是：Output argument "y" (and maybe others) not assigned during call to
"C:\Users\Administrator\Documents\MATLAB\g.m>g".

file:///C:\Users\Administrator\AppData\Roaming\Tencent\Users\1012652646\QQ\WinTemp\RichOle\G2XOW3(XK6W0I8)LQ~H]S82.jpg

jeffej

what too you use???