if __name__=="__main__":
import numpy as np
import scipy as sp
import matplotlib.pyplot as plt
theSolver=Integrator('daisy.ini')
timeVals,yVals,errorList=theSolver.timeloop5Err()
whiteDaisies=[frac[0] for frac in yVals]
thefig=plt.figure(1)
thefig.clf()
theAx=thefig.add_subplot(111)
points=theAx.plot(timeVals,whiteDaisies,'b+')
points[0].set_markersize(12)
theLines=theAx.plot(timeVals,yVals)
theLines[0].set_marker('+')
theLines[1].set_linestyle('--')
theLines[1].set_color('k')
theAx.set_title('lab 5 interactive 6')
theAx.set_xlabel('time')
theAx.set_ylabel('fractional coverage')
theAx.legend(theLines,('white daisies','black daisies'),loc='best')
timeVals,yVals,errorList=theSolver.timeloop5fixed()
whiteDaisies=[frac[0] for frac in yVals]
plt.show()
