Advertisement
Guest User

Untitled

a guest
Feb 20th, 2018
71
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Python 0.58 KB | None | 0 0
  1. import matplotlib.pyplot as plt
  2. import numpy as np
  3. # plot x as fn of t
  4. k = 5.0
  5. x0 = 10.0
  6. x = x0
  7. m = 2.0
  8. time = 6.0
  9. displacementlist = []
  10. velocitylist = []
  11. timelist = np.arange(0, 6.0, 0.01)
  12. v = 0
  13. dt = 0.001
  14. t = 0
  15. for i in range(len(timelist)):
  16.     acc = (x * -k) / m
  17.     v = v + (acc * timelist[i])
  18.     x = x + (v * timelist[i])
  19.     displacementlist.append(x)
  20.     velocitylist.append(v)
  21. print(timelist)
  22. print(velocitylist)
  23. print(displacementlist)
  24.  
  25. plt.plot(timelist, velocitylist, color = 'r')
  26. plt.plot(timelist, displacementlist, color = 'b')
  27. plt.axhline(0)
  28. plt.show()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement