matlab 解方程

bb52013147y

哪位matlab大神帮忙编个程解个方程
C_T (x,y,z)=f∙C(x,y,z)+(1-f)∙C(x,y,z_a)
C(x,y,z)=Q/U∙C(y,x)∙C(z,x)
f=0.5∙(1+θ)
θ= (∫_0^H▒C(x,y,z)dz)/(∫_0^∞▒C(x,y,z)dz)
z_a=z-z_i

式中，C_T (x,y,z)——总浓度
      C(x,y,z)——不考虑地形时的浓度
      C(x,y,z_a)——沿地形抬升的烟羽浓度
       θ——烟羽质量与总烟羽质量的比值
Q——源的泄放速率
U——有效风速值
C(y,x)和C(z,x)——分别表示水平方向和垂直方向浓度分布的概率密度函数
f——权函数
      z_a和z_i——分别表示该点地形的有效高度以及高度值
求解当∑_i 〖(C-C_B-C_T)〗^2 最小时候的x0，y0，z0；