C-FOCAL,1969

01.01 E
01.05 S PI=3.14159
01.10 A "(A+B*I)X^2+(C+D*I)X+E+G*I=0",! ,?A B C D E G ?,!!
01.20 I (A) 1.3,1.25,1.3
01.25 I (B) 1.3,7.1,1.3
01.30 S DR=C^2-D^2-4*(A*E-B*G)
01.40 S DI=2*C*D-4*(A*G+B*E)
01.50 I (DR) 1.6,1.95
01.60 S TH=FATN(FABS[DI/DR])
01.70 I (DR) 1.9
01.80 I (DI) 1.81; G 2.1
01.81 S TH=2*PI-TH
01.90 I (DI) 1.93
01.91 S TH=PI-TH; G 2.1
01.93 S TH=PI+TH; G 2.1
01.95 I (DI) 1.97
01.96 S TH=PI/2; G 2.1
01.97 S TH=3*PI/2

02.10 S Z=FSQT(FSQT[DR^2+DI^2])
02.20 S DR=Z*FCOS(TH/2)
02.30 S DI=Z*FSIN(TH/2)
02.40 S R(1)=DR-C; S I(1)=DI-D
02.50 S R(2)=-DR-C; S I(2)=-DI-D
02.60 F N=1,2; D 4
02.70 F N=1,2; T "W",%1,N,"=",%4.03,RF(N),"+",IF(N)," *I",!
02.80 T !!!; G 1.01

04.10 S VF(N)=[2*A*R(N)+2*B*I(N)]/(4*A^2+4*B^2)
04.20 S IF(N)=(+2*A*I(N)-2*B*R(N))/[4*A^2+4*B^2]

07.10 T "THAT'S NO QUADRATIC!!!",!!!
07.20 G 1.01
*