求y=a1.*exp(-sqrt((a2-lab)^2/a5 (a3-cap)^2/a6 (a4-xin)^2/a7))