wangtao1211
发表于 2014-5-31 14:33
谢谢啦谢谢啦谢谢啦谢谢啦谢谢啦
ice_and_fire
发表于 2014-7-25 22:06
clc up=@(x)sqrt(1-x.^2/2);%y积分上限 low=@(x)-sqrt(1-x.^2/2)';%y积分下限 fun=@(y,x)exp(-x.^2/2).*sin(x.^2+y);%积分函数 numint=quad2dggen(fun,low,up,-.5,1,eps)%eps为计算精度，比如换成1e-10，这里的eps就当前浮点数的精度 numint = 0.4119
ice_and_fire
发表于 2014-7-25 22:06
clc
up=@(x)sqrt(1-x.^2/2);%y积分上限
low=@(x)-sqrt(1-x.^2/2)';%y积分下限
fun=@(y,x)exp(-x.^2/2).*sin(x.^2+y);%积分函数
numint=quad2dggen(fun,low,up,-.5,1,eps)%eps为计算精度，比如换成1e-10，这里的eps就当前浮点数的精度
numint =
0.4119
離岸風
发表于 2014-7-28 20:19
先谢谢楼主了~
1157829863
发表于 2014-8-28 13:01
是不是真的好啊
流雨星月
发表于 2014-8-29 14:02
感谢楼主分享，顶一下
LYJA
发表于 2014-9-25 21:59
好，楼主辛苦了，给力！！！！！
fwing
发表于 2015-2-9 12:30
DDDDDDDDDDDDDDDDDDDDDDDDDD
曦洋
发表于 2015-9-8 19:31
赞，这个不错的，大家可以下下来看看的
曦洋
发表于 2015-9-8 19:31
赞，大赞特赞的啊，很不错，很好的