#!/usr/bin/python
# 5.4

import scipy as sp
from pylab import *

w = arange(0, pi, 0.01)
G = 1-exp(-1j*w)
subplot(2,1,1)
loglog(w,abs(G))
ylabel('Magnitude')

subplot(2,1,2)
semilogx(w, 180*angle(G)/pi)
xlabel('w (rad/s)')
ylabel('Phase')

show()
#axis tight, label('',18,'w (rad/s)','Phase',18)