File QUAGL.FT (FORTRAN source file)

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

	SUBROUTINE QUAGL(A,B,C,X1,X2,K)
C UNTERPROGRAMM ZUR LOESUNG QUADRATISCHER GLEICHUNGEN
C PARAMETER : A,B,C  KOEFFIZIENTEN DER NORMALFORM
C	      X1,X2  REELLE LOESUNGEN ODER
C	             REALTEIL UND IMAGINAERTEIL KOMPLEXER LOESUNGEN
C	          K  KENNZIFFER MIT FOLGENDER BEDEUTUNG :
C 		      K = 0  KEINE LOESUNG, DA A = 0
C 		      K = 1  ZWEI REELLE LOESUNGEN
C 		      K = 2  ZWEI KOMPLEXXE LOESUNGEN
C FORTRAN II
	IF(A)1011,1012,1011
1011	GL=.5/A
	DISK=B**2-4.*A*C
	IF(DISK)1013,1014,1014
1014	K=1
	WURZ=SQRT(DISK)
	X1=(-B+WURZ)*GL
	X2=(-B-WURZ)*GL
	RETURN
1013	K=2
	WURZ=SQRT(-DISK)
	X1=-B*GL
	X2=WURZ*GL
	RETURN
1012	K=0
	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