In [1]:
import imp
import modular_theory as mt
import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
plt.style.use('fivethirtyeight')
%matplotlib inline
In [5]:
imp.reload(mt)
Out[5]:
In [6]:
procedure = "S1 60T O.1 s1"
simulation = mt.Simulation(procedure)
In [7]:
simulation.run()
In [8]:
np.set_printoptions(formatter={'float': '{: 0.3f}'.format})
simulation.array_notation[:,:,0]
Out[8]:
In [9]:
from numpy import linspace
time = linspace(0, 2*7200, 2*int(7200 / 0.01))
plt.plot(time, simulation.reference_memory)
plt.xlim((0, 70))
Out[9]:
In [10]:
plt.plot(time, simulation.strength_memory)
plt.xlim((0, 70))
Out[10]:
In [11]:
plt.plot(time, simulation.packet_initiation_rate)
plt.xlim((0, 70))
Out[11]:
In [12]:
simulation.data
Out[12]:
In [ ]:
In [ ]: