;Allg. L”sung der DGL der Ged„mpften Schwingung
DSOLVE2_IV(k/m,d/m,0,t,0,s0,0)

S(t):=s0*#e^(-t*SQRT(k^2-4*d*m)/(2*m)-k*t/(2*m))*(SQRT(k^2-4*d*m)-k)/(2*SQRT(~
k^2-4*d*m))+#e^(t*SQRT(k^2-4*d*m)/(2*m)-k*t/(2*m))*(k*s0/(2*SQRT(k^2-4*d*m))+~
s0/2)

;Masse
m:=0.5

;Federh„rte
d:=20

;Anfangsauslenkung
s0:=4

; 
"Bsp.1: Ged„mpfte Schwingung: k^2 < 4Dm"

;Reibungsfaktor
k:=0.3

; 
S(t)=80*SQRT(39910)*#e^(-3*t/10)*SIN(ATAN(SQRT(3991)/3)+SQRT(3991)*t/10)/3991

; 
"Bsp.2: aperiodischer Grenzfall   k^2 = 4Dm"

; 
k:=SQRT(40)-0.001

; 
S(t)=115397/513*#e^(-7114/1125*t)*SIN(2291/20371*t+587/33009)

; 
"Bsp.3: Kriechfall   k^2 > 4Dm"

; 
k:=200

; 
S(t)=#e^(t*(6*SQRT(1110)-200))*(20*SQRT(1110)/333+2)+#e^(-t*(6*SQRT(1110)+200~
))*(2-20*SQRT(1110)/333)