#!/usr/bin/env python2.5
from matplotlib.pyplot import *
from numpy import sin, exp,  absolute, pi, arange
from numpy.random import normal

def f(t):
    s1 = sin(2*pi*t)
    e1 = exp(-t)
    return absolute((s1*e1))+.05


t = arange(0.0, 5.0, 0.1)
s = f(t)
nse = normal(0.0, 0.3, t.shape) * s

plot(t, s+nse, 'b^')
vlines(t, [0], s)
xlabel('time (s)')
title('Comparison of model with data')
show()

