from visual import *
from visual.graph import *
l=5
k=89
k=k*pi/180
g=9.81
dt=0.01
vk=0
w=500
display(title='pendulum', x=0, y=0, width=w, height=w)
kulx=l*sin(k)
kuly=l*cos(k)
suf=box(pos=vector(0, 0.5/2, 0), size=(1, 0.5, 1), color=color.red)
sznur=cylinder(pos=(0,0,0), axis=(kulx, -kuly,0), radius=0.02)
kula=sphere(pos=(kulx,-kuly,0), axis=(0,0,1), radius=0.15, color=color.green)
h=w/2
gdisplay(title='chart', x=w+1, y=0, width=w, height=h)
st=gcurve(color=color.blue)
wt=gcurve(color=color.yellow)
t=0
while True:
rate(100)
vk=vk-dt*g/l*sin(k)
k=k+vk
kulx=l*sin(k)
kuly=l*cos(k)
kula.pos=vector(kulx,-kuly,0)
sznur.axis=vector(kulx, -kuly,0)
s=l*k
t +=dt
st.plot(pos=(t,s))
wt.plot(pos=(t,vk))