File: candelator.ode

package info (click to toggle)
xppaut 8.0%2Bdfsg-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 7,332 kB
  • sloc: ansic: 74,690; makefile: 127; sh: 92
file content (16 lines) | stat: -rw-r--r-- 241 bytes parent folder | download 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
 
!g=980
par k=.1
init m1=60,m2=60
r1=m1/4
r2=m2/4
d=(-r1*m1+r2*m2)/(m1*r1^2+m2*r2^2)
w'=.75*g*sin(th)*d-k*w
th'=w
m1'=scale*f(th1)
m2'=scale*f(th2)
par scale=1e-4
f(t)=if(t<-1.308)then(.556+.8307*t)else(if(t<.6982
@ meth=cvode
aux drive=d
d