#!/usr/bin/python # 4.1 import scipy as sp from pylab import * a = 0.9 n = sp.arange(-10,21) y = a**n*(n>=0) stem(n,y,markerfmt='.') xlabel('Time(n)') ylabel(r'$0.9^n1(n)$') show()