Solving transcendental equation for a complex variable using newtzero file.

I am trying to solve the complex transcendental equation for the variable (ZsTM). When I tried using your newtzero file, I got an error message about the arguments of the command. It would be so nice of you if you can help me in this issue.

My code is :

     function [output] = EBG2(epsilon_r , freq , pe ,gap ,h,rvia)
    speed = 3e8 ; 
    e0 = 8.85e-12 ; 
    u0 = 4*pi*1e-7 ;
    z0 = 377;
    lambda = speed./ freq 
    k0 = 2*pi./lambda 
    kr=k0*sqrt(epsilon_r)
    epsilon_eff = (epsilon_r +1)/2
    k_eff = k0.*sqrt(epsilon_eff)
    w=2*pi*freq;
    kp=(1/(pe*sqrt((1/2*pi)*log((pe^2)/(4*rvia*(pe-rvia))))))
    zeff=z0/sqrt(epsilon_eff)
    alpha = (k_eff*pe/pi)* log(1/(sin(((pi*gap/(2*pe))))))
    zg=(-1*j*(zeff./(2*alpha)))
  % The transcendental equation for ZzTM

    ZsTM=(j*w*u0*tan(kr*h*((1-((1/epsilon_r)*(1-(ZsTM/z0)^2)))^0.5))*(kr^2*((ZsTM/z0)^2)-kp^2)*(1-((1/epsilon_r)*(1-(ZsTM/z0)^2))))/(kr*((1-((1/epsilon_r)*(1-(ZsTM/z0)^2)))^0.5)*(kr^2-kp^2))
    output=ZsTM

And in the main program I pass the parameters

epsilon_r=6.15; freq=60e9; pe=1.25e-3; gap=0.3e-3; h=0.254e-3; rvia=0.08e-3;