MATLAB: The following code is not working….cant reduce it to integral 3

double integralMATLABsymbolicSymbolic Math Toolbox

syms x y z;
f1=@(x)(0.93*dirac(x-0.01)+0.07*dirac(x-100));
f2=@(y)(0.93*dirac(y-0.01)+0.07*dirac(y-100));
intfn=matlabFunction(exp((1-(1/(1+(2.87*10^(-5)*z/(x*y)))))*z*-6.3*10^(-3)));
innerintegral=@(x,y)int(intfn,z,0,inf);
mainfun=@(x,y)innerintegral(x,y)*0.4*exp(-5*10^(-10)/(x*y))*f1(x)*f2(y);%this is the main function to be integrated which utilizes f1, f2 and innerintegral above
resultintermediate=int(mainfun,x,0.01,inf);
result=int(resultintermediate,y,0.01,inf);
double(result)

Best Answer

It‘s just a guess just change all your upper limits to 1 meaning change all inf to 1 you will get inbetween 0 and 1
Related Question