Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- import pylab
- import random
- import ps8_help
- if __name__ == '__main__':
- random.seed(0)
- x, y = ps8_help.simulationWithDrug(100, 1000, 0.1, 0.05,
- {'guttagonol': False}, 0.005, 100, timing=True)
- pylab.plot(x, 'r', label='Total')
- pylab.plot(y, label='Resistant')
- pylab.title('Resistence and single medication from midpoint')
- pylab.xlabel('Timestep')
- pylab.ylabel('Virus population size')
- pylab.legend()
- pylab.show()
Advertisement
Add Comment
Please, Sign In to add comment