压电仿真，有一点点小问题

代码说明：

%主函数 function  text_greens%通过广义格林函数求有效介电常数clear;%弹性刚度常数矩阵c1=[2.0961.205 1.046 0     0     0;...   1.205 2.096 1.046 0     0    0;...   1.046 1.046 2.106 0     0    0;...   0     0     0    0.423 0     0;...   0     0     0    0     0.423 0;...   0     0     0    0     0   0.4455];  % Elasticity Matrix of ZnO (10^11N/m^2)  c2=[11.5310.863  0.863   0      0     0;...