[問題]有關積分問題想請教大家一下

[reddragon00 (00)]

請教一下各位大大
剛接觸matlab不久,忽然要寫出這樣的function對我而言有些吃力

對一個連乘的標準常態(theta*x)*Gammapdf做0.001到T的積分

我寫的 function 如下,但是一直跑不出結果,請各位大大能為我指點一下
萬分感謝


function [Pd Ps]=defaul_survival_prob(T,nu,sigma,theta,)
N=length(theta);
    for j=1:N
    if j==1

       phi(j)=@(x)normcdf((theta(j).*x)./(sigma(j).*sqrt(x)));

       else j>1;

       phi(j)=@(x)normcdf((theta(j).*x)./(sigma(j).*sqrt(x)))...
          *normcdf((Li(j-1)-m(j-1)*T-VG_theta(j-1).*x)./(sigma(j-1).*sqrt(x)));
    end
    end
Fzd=@(phi)phi(N).*gampdf(T/nu,1/nu);
Pd=quad(Fzd,0.001,T);

--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 220.132.65.140
推 sunev:Nonscalar arrays of function handles are not allowed;
→ sunev: use cell arrays instead.
→ reddragon00:可以請教一下如何改善這個問題嗎?