matlab 整理公式,已知一段连续的数据，如何用matlab实现公式里参数的优化...

func=@(a,f)(((((a(1)+a(2)*(10^-2/a(3))^2*(((10^-2/a(3))^2-(10^-2/f)^2)/((((10^-2/a(3))^2-(10^-2/a(3))^2))^2+(a(4)^2*(10^-2/a(3))^2*(10^-2/f)^2))))^2+(a(2)*(10^-2/a(3))^2*(a(4)*((10^-2/f))*(10^-2/a(3))/((((10^-2/a(3))^2-(10^-2/a(3))^2))^2 + a(4)^2*(10^-2/a(3))^2*(10^-2/f)^2 )))^2)^0.5+(a(1)+a(2)*(10^-2/a(3))^2*(((10^-2/a(3))^2-(10^-2/f)^2)/((((10^-2/a(3))^2-(10^-2/a(3))^2))^2+(a(4)^2*(10^-2/a(3))^2*(10^-2/f)^2))))/2)^0.5-1)^2+(((a(1)+a(2)*(10^-2/a(3))^2*(((10^-2/a(3))^2-(10^-2/f)^2)/((((10^-2/a(3))^2-(10^-2/a(3))^2))^2+(a(4)^2*(10^-2/a(3))^2*(10^-2/f)^2))))^2+(a(2)*(10^-2/a(3))^2*(a(4)*((10^-2/f))*(10^-2/a(3))/((((10^-2/a(3))^2-(10^-2/a(3))^2))^2 + a(4)^2*(10^-2/a(3))^2*(10^-2/f)^2 )))^2)^0.5-(a(1)+a(2)*(10^-2/a(3))^2*(((10^-2/a(3))^2-(10^-2/f)^2)/((((10^-2/a(3))^2-(10^-2/a(3))^2))^2+(a(4)^2*(10^-2/a(3))^2*(10^-2/f)^2))))/2))/(((((a(1)+a(2)*(10^-2/a(3))^2*(((10^-2/a(3))^2-(10^-2/f)^2)/((((10^-2/a(3))^2-(10^-2/a(3))^2))^2+(a(4)^2*(10^-2/a(3))^2*(10^-2/f)^2))))^2+(a(2)*(10^-2/a(3))^2*(a(4)*((10^-2/f))*(10^-2/a(3))/((((10^-2/a(3))^2-(10^-2/a(3))^2))^2 + a(4)^2*(10^-2/a(3))^2*(10^-2/f)^2 )))^2)^0.5+(a(1)+a(2)*(10^-2/a(3))^2*(((10^-2/a(3))^2-(10^-2/f)^2)/((((10^-2/a(3))^2-(10^-2/a(3))^2))^2+(a(4)^2*(10^-2/a(3))^2*(10^-2/f)^2))))/2)^0.5+1)^2+(((a(1)+a(2)*(10^-2/a(3))^2*(((10^-2/a(3))^2-(10^-2/f)^2)/((((10^-2/a(3))^2-(10^-2/a(3))^2))^2+(a(4)^2*(10^-2/a(3))^2*(10^-2/f)^2))))^2+(a(2)*(10^-2/a(3))^2*(a(4)*((10^-2/f))*(10^-2/a(3))/((((10^-2/a(3))^2-(10^-2/a(3))^2))^2 + a(4)^2*(10^-2/a(3))^2*(10^-2/f)^2 )))^2)^0.5-(a(1)+a(2)*(10^-2/a(3))^2*(((10^-2/a(3))^2-(10^-2/f)^2)/((((10^-2/a(3))^2-(10^-2/a(3))^2))^2+(a(4)^2*(10^-2/a(3))^2*(10^-2/f)^2))))/2));

[a,r,j]=nlinfit(f,y,func,[1.5^2 1*10^-4 9 0.005]);

当我运行这个程序的时候，出现了下面的错误