Re: [問題] Ei?
※ 引述《lwhsyh (KINKI)》之銘言：
: syms x real
: f=((sin(2*x)))/x)*exp(6*x);
: F=int(f,x,-10000,10000)
: F =
: 1/2*i*Ei(1,-60000-20000*i)-1/2*i*Ei(1,-60000+20000*i)
: -1/2*i*Ei(1,60000+20000*i)+1/2*i*Ei(1,60000-20000*i)
: -1/2*i*log(3+i)+1/2*i*log(3-i)+1/2*i*log(-3-i)-1/2*i*log(-3+i)
: 請問各位高手先進
: 我這積分出來的答案表示什麼
: Ei??
: 有無辦法可以表示數字
: 感謝大家幫忙
這邊的Ei是Symbolic Math Toolbox裡面的mfun的東西
使用方法為mfun( 'Ei' , n , z )即可算出數字
其中n必須大於等於0
real(z)必須大於等於0
可是你的答案裡面有-60000-20000*i這種東西
real(z)是小於0的
這樣子跑下去mfun( 'Ei' , 1,-60000-20000*i )的答案是NaN - Infi
而就算real(z)是大於0的拿去跑mfun( 'Ei' , 1,60000-2000*i )答案也是0而已
而這個積分我用quad和quadl去處理也是跑不出來（答案是NaN）
使用高斯基分法也是跑出來NaN
（不過我沒有幫你做變數變換轉換區間，只有稍為算一下而已）
你這個可能真的沒辦法積分
Symbolic Toolbox常常會跑出來mfun的東西
要查詢各mfun細節的話可以help mfunlist
--
Deserves death! I daresay he does. Many that live deserve death. And some die
that deserve life. Can you give that to them? Then be not too eager to deal out
death in the name of justice, fearing for your own safty. Even the wise cannot
see all ends.
Gandalf to Frodo
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.120.25.235
推
09/14 14:00, , 1F
09/14 14:00, 1F