lang/dynamo/1.8.7/src/dyn.steifedgl


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15

NOTE STIFF EQUATIONS  SIEHE: SIMULATION AUGUST 1980, SEITE 38
L    Y1.K=Y1.J+DT*(-21*Y1.J+19*Y2.J-20*Y3.J)
L    Y2.K=Y2.J+DT*(+19*Y1.J-21*Y2.J+20*Y3.J)
L    Y3.K=Y3.J+DT*(+40*Y1.J-40*Y2.J-40*Y3.J)
N    Y1=1
N    Y2=0
N    Y3=-1
NOTE KONSTANTEN MUESSEN GEEIGNET GEWAEHLT WERDEN:  DT  SEHR KLEIN
C LENGTH=20
C DT=.01
C PRTPER=1
C PLTPER=1
PRINT Y1,Y2,Y3
PLOT  Y1,Y2,Y3