import pylab
import math

from sistemi import *

t = 0.0
delta_t = 1e-3 

tempi = [ ]
output = [ ]
u = 1

sistema = BigSystem()

while t <= 8:
    output.append(sistema.evaluate(u, delta_t))
    tempi.append(t)
    t = t + delta_t


pylab.figure(1)
pylab.plot(tempi, output, 'r-+', label='output')
pylab.xlabel('time')
pylab.legend()

pylab.show()