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)); endfunction 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
Best Answer