MATLAB: The following program involving double integration hangs…Tried everything in R2014a

MATLABmultiple integrals

syms x y z;
f1=@(x)(0.6*dirac(x-0.01)+0.4*(x-100));
f2=@(y)(0.6*dirac(y-0.01)+0.4(y-100));
intfn=matlabFunction((1-(1/(1+(2.87*10^(-5)*z/(x*y)))))*z);
innerintegral=@(x,y)(int(intfn,z,0,inf)*-6.3*10^(-3));
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

Lookup integral3() : 
syms x y z;
f1=(0.6*dirac(x-0.01)+0.4*(x-100));
f2=(0.6*dirac(y-0.01)+0.4*(y-100));
intfn=((1-(1/(1+(2.87*10^(-5)*z/(x*y)))))*z);
mainfun=matlabFunction(((1-(1/(1+(2.87*10^(-5)*z/(x*y)))))*z)*...
    0.4*exp(-5*10^(-10)/(x*y))*(0.6*dirac(x-0.01)+...
    0.4*(x-100))*(0.6*dirac(x-0.01)+0.4*(x-100)));
result=integral3(mainfun,0.01,1000,0.01,1000,0,1000) % the upper limits where reduced
Related Question