File FGS.FT (FORTRAN source file)

Directory of image this file is from
This file as a plain text file

	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



Feel free to contact me, David Gesswein djg@pdp8online.com with any questions, comments on the web site, or if you have related equipment, documentation, software etc. you are willing to part with.  I am interested in anything PDP-8 related, computers, peripherals used with them, DEC or third party, or documentation. 

PDP-8 Home Page   PDP-8 Site Map   PDP-8 Site Search