// function to find SWR,
function[SWR]=VSWR(tao)
  SWR=(1+tao)/(1-tao)
endfunction