高斯烟团模型matlab实现 求大神指点

function [C Cs] = puff_model(x,y)
syms Q bx by0 bz0 d u H z Zo temp1 temp2 t;
Q = input('请输入污染量(kg)Q = ');
u = input('请输入计算风速(m/s)u = ');
Zo = input('请输入地面粗糙长度(m)Zo = ');
H = input('请输入有效源高(m)H = ');
z = input('请输入预测点有效高度(m)z = ');
d = input('请输入计算精度(m)d = ');
t = input('请输入扩散时间(s)t = ');
[x,y]=meshgrid(0:d:4000,-4000:d:4000);%定义解空间和计算精度.
by0=0.08*x.*(1+0.0001*x).^(-1/2);%计算y轴向的基本扩散参数
bz0=0.06*x.*(1+0.0015*x).^(-1/2);%计算z轴向的基本扩散参数
%高斯扩散浓度模拟计算%
temp1 = (2 * pi)^(2/3)*by0.*bz0;
temp2 = -((x-u.*t).^2)./2 -(y.*y)./(by0.*by0)-...
((z-H).^2)./(2*(bz0.*bz0))-((z+H).^2)./(2*(bz0.*bz0));
C = Q./temp1 .* exp(temp2);
Cs = input('输入所有求解浓度(mg/m^3):');%所输人数以”[”和”]”作为开头和结束。
contour(x,y,C,Cs);%以x,y,C变量分别作为xyz的轴,绘制扩散浓度分布图
shading interp;
colorbar;
grid ;
xlabel('x轴距离(m)');
ylabel('y轴距离(m)');
title('气体下风向浓度分布图');