[求助]非线性方程计算的问题

[qingtian]

h=[0
100
200
300
400
500
600
700
800
900
1000];
x0=[0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.0
1.1]
p=zeros(length(h),1);
for i=1:length(h)
y=h(i);
f=inline('log(x^2+(bi)*(Tav)^2)-2*(ai)*y/(Tav)-c1','x','y');
f1=subs(f);
p(i)=fzero(inline(f1),x0);
end
plot(h,p)
求教高手，我想求出每个h对应的压力值，这个程序运行老有错误，inline部分，请告诉到底怎么改吧，谢谢了:tongue: