最近遇到重積分的問題
我想要積分的函數如下
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)
但是還是沒有改進..
不知道大家有沒有什麼好的方法 還是我的寫法可以改進
謝謝