Numerical solution for convolution 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,

Answer 1

[Torsten]

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.