Integral calling another integral

Question

Orongo 2018-11-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/427909-integral-calling-another-integral

评论： Star Strider 2018-11-4

Hi I get an error when calculating an integral. The error is caused by my "nested" integral. Here is what I have written so far and running it show you an error after the calculation of f_lx. How can I have this running?

age=61;  
f_AD=@(t) f_lx(t,age);
AD(age,1)=integral(f_AD,age,106);
function res=f_lx(t,age)
 param_1955 = [4.67255690389772/1000, 0.0192034319814117/100000, 1.47616811690684/10];
 mu_1955=@(x) f_mu(x,param_1955); % my error appears here
 res = exp(-integral(mu_1955,0,age)); 
end
function res=f_mu(x,param)
 a=param(1); b=param(2); c=param(3);
 res = zeros(size(x));
 ind = x>100;
 res(ind) = a+b*exp(c*100)+(x(ind)-100)*0.001;
 res(~ind) =a+b*exp(c*x(~ind));
end