File: ddho.ode

package info (click to toggle)
plotutils 2.4.1-15
  • links: PTS
  • area: main
  • in suites: etch, etch-m68k
  • size: 11,072 kB
  • ctags: 6,952
  • sloc: ansic: 76,305; cpp: 12,402; sh: 8,475; yacc: 2,604; makefile: 894; lex: 144
file content (32 lines) | stat: -rw-r--r-- 644 bytes parent folder | download | duplicates (14) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
 
# You may run this example by doing:
#
#    ode < ddho.ode | graph -T X -C
# 
# or alternatively, to get a real-time plot,
#
#    ode < ddho.ode | graph -T X -C -x 0 25 -y -0.5 1 0.5
# 
# This example simulates a dam-driven harmonic oscillator (DDHO) with
# damping, the equation for which is
# 
# y'' = -k/m * y - R/m * y' + cos(w*t)
# 
# If R^2 > 4km, motion is overdamped
# If R^2 = 4km, motion is critically damped
# If R^2 < 4km, motion is damped
# 
# With the choice of parameters below, motion is damped but not
# overdamped.

y'	= vy
vy'	= -k/m * y - R/m * vy + cos(w*t)

y	= 1
vy	= 0
k	= 1
m	= 1
R	= 0.5
w	= 2*PI

print t, y
step 0,25