dblquad and double integral problem

i want to use dblquad command to solve my problem but it is a symbolic Algebra matrix which is very big. so i use a simple 2by2 matrix to disp it .i wanto to expand 16by16 matrix .can anyone help me.

i have a file which is my 16by16 matrix.thanks,i want to double integral it(the varaible is ks and yt and region is -1<=ks<=1、-1<=yt<=1)

example:

a={@(x,y)(2.*x+y.^3),@(x,y)(x+2*y.^3);@(x,y)(x+3*y.^3),@(x,y)(x+4*y.^3)};

for i=1:4 matrix(i)=dblquad(a{i},-1,1,-1,1); end reshape(matrix,2,2)