clc s=1.118//m h=1.00//m a=sqrt(1+(s/h)^2) disp(a,"a is= ")