Search
#台北飲酒記 哥調的是酒，飲的是品味
上星期到台北出差，忙完一輪，約朋友到酒吧坐坐，放鬆下。今次去咗大安區的EAST END酒吧，這裡裝修有格調，藏酒多，靜靜地又可傾到計，酒保OK有型，手勢熟練，試了他們剛攞獎的cocktail「哈利法」，紅茶加水果配製，味道清新，好味。而更好的是，這間吧是Asia Miles夥伴商戶，即是話消費可同時儲分，即時拿出電話開番個app，向侍應show番我的Asia Miles 亞洲萬里通會員號碼，再用渣打Asia Miles信用咭埋單，儲埋多重積分，實行「Pay Once Earn Twice」。
#Asiamiles #kirkyiptravel #asiatop50bars
【#余博講物理】天文基礎數學工具 球面極座標 多重積分
【#余博講物理】天文數學 球面極座標 多重積分
唔係呃你 唔係點你 唔係講笑 詢眾要求 今晚愚人節 11pm 教大家研究天文學嘅基本數學工具 #嗚呀終於講數學 #究竟係講天文定數學
#physics #fysik #astronomy #astronomi #mathematics #matematik #phycat
多重積分 問題. 這一題積分範圍由accos（√3/2）至accos（1/2），可是Φ不是z跟半徑夾角嗎？這樣應該要是acsin吧？ 1. 留言2. 右轉數學版. B22020年6月8日. ... <看更多>
請教一下是否有人能夠解釋一下多重積分變數替換?我搞不懂如何設改變的變數eg. u=.... V=....假設integral 下有R region 要怎樣找出改變後的S region如這題我嘗試 ... ... <看更多>
多重積分 在 [程式] R code 多重積分的問題- 看板Statistics - 批踢踢實業坊 的推薦與評價
[軟體程式類別]:
R code
[程式問題]:
多重積分
[軟體熟悉度]:
高(1年以上)
[問題敘述]:
目前要做一個多重數值積分的結果
內層兩個函數
1. f1(v1,y1) 要把 v1 積掉
2. f2(v2,w,y1) 要把 w 積掉
3. 函數 1 & 2 積完後相乘，接著把 y1 積掉
從頭到尾 w 都是 input 定數
積分結果出現了錯誤訊息
"evaluation of function gave a result of wrong length"
[程式範例]:
n1=n2=5;CPU1=CPU2=1
A1=2/sqrt(2*pi/(9*n1))*gamma((n1-1)/2)*(2/(n1-1))^((n1-1)/2)
A2=2/sqrt(2*pi/(9*n2))*gamma((n2-1)/2)*(2/(n2-1))^((n2-1)/2)
A=A1*A2
F1=function(y1)
{
f1=function(v1)
{
exp(-9*n1/2*(v1*y1-CPU1)^2)*(v1^(n1-1))*exp(-(n1-1)/2*v1^2)
}
return(f1)
}
int1=function(y1) {integrate(F1(y1),lower=0,upper=Inf)$value}
F2=function(w,y1)
{
f2=function(v2)
{
exp(-9*n2*(v2*(w+y1)-CPU2)^2/2)*v2^(n2-1)*exp(-(n2-1)/2*v2^2)
}
return(f2)
}
int2=function(w,y1) {integrate(F2(w,y1),lower=0,upper=Inf)$value}
F3=function(w)
{
f3=function(y1)
{
A*int1(y1)*int2(w,y1)
}
return(f3)
}
int3=function(w) {integrate(F3(w),lower=0,upper=Inf)}
先寫 Global，都是一樣的寫法，int1, int2 函數都沒問題。
但 int3 結果就會出毛病，似乎是資料型態的問題，嘗試了一整天
sapply, mapply, Vectorize 等函數，還是無法克服
希望版上能有大大提供點想法
儘管只是個簡單的三變數 x,y,z 積分也可以，感恩~
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.113.108.134
... <看更多>