' Wave simulation
' made by Svjatoslav Agejenko
' in 2003.12
' H-Page: svjatoslav.eu
' E-Mail: svjatoslav@svjatoslav.eu
 
SCREEN 13

DIM SHARED yh(1 TO 300)
DIM SHARED yhs(1 TO 300)


FOR x = 1 TO 300
yhs(x) = 0
yh(x) = 50
NEXT x



1

CLS
FOR x = 1 TO 300
'  LINE (x, 0)-(x, 200 - yh(x)), 0
'  LINE (x, 200 - yh(x))-(x, 200), 15
  PSET (x, 200 - yh(x)), 31
NEXT x


FOR x = 10 TO 290
  hk = (yh(x - 1) + yh(x + 1) + yh(x + 2) + yh(x - 2)) / 4
  yhs(x) = yhs(x) + (hk - yh(x)) / 5
  yhs(x) = yhs(x) / 1.01
NEXT x

FOR x = 10 TO 290
  yh(x) = yh(x) + yhs(x)
'  yh(x - 1) = yh(x - 1) + yhs(x) / 2
'  yh(x + 1) = yh(x + 1) + yhs(x) / 2
'SOUND 0, .05
NEXT x

'FOR x = 10 TO 290
'  yh(x) = (yh(x) * 100 + yh(x + 1) + yh(x - 1)) / 102
'NEXT x

IF RND * 100 < 2 THEN
  p = RND * 200
  s = RND * 10 + 2
  FOR x = 0 TO 3.14 STEP 3.14 / s
     yh(p) = yh(p) + SIN(x) * s * 3
     p = p + 1
  NEXT x
END IF

IF INKEY$ <> "" THEN SYSTEM
GOTO 1