syms beta
epsilon_r=2.1;
f=(7e12);
w=(2*3.1416*f);
surface_conductivity = (0.00012164049 - 0.0053497497i);
c=(3e8);
epsilon_naught=(8.854e-12);
solve((1/(sqrt((beta).^2 - (w/c).^2))) + ((epsilon_r)/(sqrt((beta).^2 - ((w/c).^2)*epsilon_r))) + i*(surface_conductivity/(epsilon_naught*w)),beta)
Do verify?