% Define the functionf = @(x,y) x.^2-15*y.^2+ x.*y.^3;% Use integral2 to evaluate the double integralI = integral2(f,0,4,-2,2);% Display the result rounded to 3 decimal placesfprintf('The solution of the double integral is %.3f', I);