In [ ]:

    

%matplotlib inline

import numpy as np
import matplotlib.pyplot as plt
import pdb

dt = 0.5
xk_1 = 0
vk_1 = 0
#ak_1 = 0 #gamma term - acceleration component

a = 0.1 # alpha
b = 0.1 # beta
# g = 0.001 # gamma

iteration = 0

fig = plt.figure()
#fig.show()
fig.canvas.show
plt.axis([0, 200, 0, 1.5])

time = np.linspace(0,200)

for i in time
    xm[i] = np.random.rand(1)

    xk[i] = xk_1 + (vk_1 * dt)
    vk[i] = vk_1 # + (ak_1 * dt)
    #ak = (vk - vk_1) /dt

    rk = xm - xk

    xk[i] += a*rk
    vk[i] += (b*rk)/dt
    #ak += (2*g*rk)/np.power(dt,2)

    xk_1 = xk[i]
    vk_1 = vk[i]
    #ak_1 = ak

    #print("xk = " ,xk, "vk = " ,vk, "ak = ",ak)
    #iteration += 1
    #plt.axis([0, 20+iteration, 0, 1.5])  
plt.plot(time,xm,'k+')
plt.plot(time,xk,'b*')
#plt.pause(0.05)
#plt.draw()


    

In [ ]:

    

%matplotlib inline

import numpy as np
import matplotlib.pyplot as plt
import pdb

abc = np.linspace(0,20)
plt.plot(abc,sin(abc))
show()


    

In [ ]: