clear all
clc
syms al bt
for i=1:9
    for j=1:9
        K(i,j)=2*j*al+i*bt;
    end
end
I=zeros(9,9);
al=1:1/10:2;
bt=2/2/10:4;
for p=1:9
    for q=1:9
        C=K(p,q);
        I(p,q)=trapz(al,trapz(bt,C));
    end
end