MATLAB: Numerical solution for convolution integral

convolutiondouble integral

Hi, I need to solve numerically the attached double integral, where x1 and x2 are variables, a is constant. x1 comprises the boundary of the first integration.
Can someone suggest how to do this?
Best Regards,

Best Answer

a=...;
L=...;
fun=@(x,y)(0.5*x.^3+x.^2).*exp(-x/a).*(0.5*y.^3+y.^2);
upper_y=@(x)x;
q=integral2(fun,0,L,0,upper_y);
Best wishes
Torsten.
Related Question