矩阵N是一个包含多项式函数的2x2矩阵：
N=[(1-x).*(1-y), 0 ;
0 , (1-x).*(1-y) ];
在x区间定义y和[0:1]以制作方形元素
我还有一个运算符矩阵D
D= [d/dx, 0
0, d/dy];
我想评估：
的积分int(D*B)dA
所以我的问题是，我无法将x和y定义为：
x= 0:dx:1
y= 0:dy:1
因为那时我会将矩阵中length(0:dx:1)的向量存储为元素，并且由于尺寸不合适，矩阵乘法将无法解决。
那么将这些多项式作为函数元素存储然后用矩阵微分算子D