function showarg(ivar,iinp)
% make a semi-log plot of the phase of variable ivar responding to input iinp

  freq=logspace(5,9,200);
  for i = 1:200
    gain(i) = equations(ivar,iinp,freq(i));
  end
  semilogx(freq,angle(gain))
