谢谢,辛苦啦~
不错不错~~~拿来借鉴下
非常感谢楼主的福利!
非常感谢楼主的福利!
这都行 笑了
楼主好银啊~~~
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
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
{:3_41:}{:3_42:}{:3_48:}
{:3_59:}{:3_59:}{:3_59:}{:3_59:}{:3_59:}{:3_59:}