"Tao_Real","Tao_Imaginary""Tao_Real","Tao_Imaginary"// function for reflection coefficient.
function[tao]=reflection_coefficient(Zl,Zo)
  tao=(Zl-Zo)/(Zl+Zo);
endfunction