[啦啦啦啦11221]

这是我编的 看看哪里不对呢 公式在附件。有偿求指导
for j=1:n
T=exp(-a1.*(p*x0-x).^2);
E=T;
E=T+E;
end
for i=1:m
Y=exp(-a1.*(o*x0-x).^2);
H=Y;
H=Y+H;
end
[x,z]=meshgrid(x,z);
[E,H]=meshgrid(E,H);
B3=1/2*Q*(E+H);