[發問] 二重積分的運算問題

[thereocr (眼睛閉上就可以睡著)]

 2π 2π             a1*a2*cos(α-β)
∫　∫  _____________________________________dαdβ
 0　 0  (a1^2+a2^2+d^2 - 2*a1*a2*cos(α-β))

以上是題目!!

我試了好久用Matlab記算這個方程式都沒算出來

不知道問題出在哪裡...

希望高人可以指點我ㄧ下

--------------以下程式-------

----------------設定function----

function out=integrn(X,Y)
D=0.8;
a1=55;
a2=30;
up=(a1*a2*cos(X-Y));
down=sqrt(a1^2+a2^2+D^2-2*a1*a2*(cos(X-Y)));
out=up/down;
end
------------------------------------

---------運算--------------

Xmin=0;
Xmax=2*pi;
Ymin=0;
Ymax=2*pi;
ans=dbquad(@integrn,Xmin,Xmax,Ymin,Ymax)

--------------------------------------


--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 150.69.51.201
※ 編輯: thereocr        來自: 150.69.51.201        (12/18 16:59)
※ 編輯: thereocr        來自: 150.69.51.201        (12/19 09:59)
推 webbber:把integrn寫在ablquad(@(x,y)函數,Xmin,Xmax,Ymin,Ymax)
→ webbber:雖然這樣會變得很長但是不會有錯 函數裡除了XY不能有其他
→ webbber:變數 ex: a1,a2,D
clear
Xmin=0;
Xmax=2*pi;
Ymin=pi;
Ymax=3*pi;
M12=dblquad(@(X,Y)
(55*30*cos(X-Y))/sqrt(55^2+30^2+0.8^2-2*55*30*cos(X-Y)),Xmin,Xmax,Ymin,Ymax);

這是我的程式,還是會出現錯誤!!
請問我應該怎麼改呢?
如果方便的話W大可以幫我試試看嗎??
※ 編輯: thereocr        來自: 150.69.51.201        (12/19 16:10)
→ webbber:./(sqrt(55^2+30^2+0.8^2-2*55*30*cos(X-Y)))
→ webbber:除號前面要有點 sqrt()外面要有()包起來
W大謝謝 我弄出來了XDDD 真的超級感謝你的
※ 編輯: thereocr        來自: 150.69.51.201        (12/20 11:25)