最近遇到重积分的问题
我想要积分的函数如下
1 x
S {[S f(y)dy]g(x)}dx
0 0
我的写法是 :
fun = @(x,y) f(y)g(x)
ymax = @(x) x
ans = quad2d(fun,0,1,0,ymax)
试了一些简单的例子 用手验算没问题
但是函数复杂化之后就好像怪怪的了
不知道是不是有更好的写法
对了 g(x) 里面有 1/(x-1) 也有可能在边界出了问题 导致整体误差
所以我把 范围改成 quad2d(fun,0,1-1e-15,0,ymax)
但是还是没有改进..
不知道大家有没有什么好的方法 还是我的写法可以改进
谢谢