FUNCTION FGS
C AMPLITUDEN-WERTE DER FREIEN GEDAEMPFTEN SCHWINGUNG
C LOESUNG DER DIFFERENTIALGLEICHUNG :
C	DY**2/DT**2 + A*DY/DT + B*Y = 0
C FUER DIE ANFANGSWERTE YT0=Y(T=0) UND VT0=DY/DT(T=0)
C DIE PARAMETER T,A,B,YT0 UND VT0 STEHEN IN COMMON
C FORTRAN II
	COMMON T,A,B,YT0,VT0
	H=0.5*A
	D=A*A-4.*B
	IF(D)1000,1010,1020
1000	D=SQRT(-D)
	FGS=EXP(-H*T)*(YT0*COS(D*T)+(VT0+H*YT0)*SIN(D*T)/D)
	RETURN
1010	FGS=EXP(-H*T)*(YT0+(H*YT0+VTO)*T)
	RETURN
1020	D=SQRT(D)
	C1=((D+H)*YT0+VT0)/(2.*D)
	C2=((D-H)*YT0-VT0)/(2.*D)
	FGS=C1*EXP((D-H)*T)+C2*EXP(-(D+H)*T)
	RETURN
	END