
Simulationsprogramm zur harmonischen Schwingung mit Reibung
************************************************************************
ZUNCHST EINE DER NCHSTEN ZEILEN  UND DANN "WIEDERHOLE BIS....." (siehe unten ) DOPPELT ANKLICKEN !

Anfangsbedingungen
|x|t|D|h|v|k|m|u|f|=|4|0|1|0.02|0|0|1|0|0|: ohne Dmpfung
|x|t|D|h|v|k|m|u|f|=|4|0|1|0.02|0|0.5|1|0|0|
|x|t|D|h|v|k|m|u|f|=|4|0|1|0.02|0|  1 |1|0|0|
|x|t|D|h|v|k|m|u|f|=|4|0|1|0.02|0|1.5|1|0|0|
|x|t|D|h|v|k|m|u|f|=|4|0|1|0.02|0|2|1|0|0|: aperiodischer Grenzfall
|x|t|D|h|v|k|m|u|f|=|4|0|1|0.02|0| 3|1|0|0|

wiederhole bis t>13
a=(-D*x-k*v)/m
b=a+(a-A)/2*f
f=1
x=0.5*b*h^2+v*h+x
v=b*h+v
A=a
t=t+h
_t;x;;1
zurck

w=(D/m-(k/m)^2/4)
p=arctan(2*wrz(w)*m/k)
f=0.5*k/m+wrz(-w)
e=0.5*k/m-wrz(-w)
f(t)=4/sin(p)*exp(-k*t/(2*m))*sin(wrz(w)*t+p);;10:  Gedmpfte Schwingung
f(t)=4*0.5*(D/m/wrz(-w))*exp(-0.5*k/m*t)*(exp(wrz(-w)*t)/e-exp(-wrz(-w)*t)/f);;10: Kriechfall
f(t)=4*exp(-0.5*k/m*t)*(1+0.5*k/m*t);;10 : Aperiodischer Grenzfall
