function z=txline(x,a2,b)
z=sqrt(a2 + b*x.^2);