Zeichensatz UTF-8+Z1013()+CTRL()-ohne Umlaute(äöüß)
10 CLS
20 INPUT "Grad des Polynoms:";N
30 PRINT:M=N:DIM A(N),B(N),C(N)
40 INPUT "Koeff.hoechster Potenz:";A0
50 FOR K=1 TO N
60 PRINT:PRINT "A(";K;")=";
70 INPUT A(K)
80 A(K)=A(K)/A0
90 NEXT K:PRINT:PRINT
100 C(1)=1
110 IF N=1 THEN PRINT:PRINT "X( 0 )=";-A(1):END
120 IF N=2 THEN P=A(1):Q=A(2):GOSUB 400:END
130 P=1:Q=0
140 B(1)=A(1)-P
150 B(2)=A(2)-Q-P*B(1)
160 FOR K=3 TO N
170 B(K)=A(K)-Q*B(K-2)-P*B(K-1)
180 NEXT K
190 C(2)=B(1)-P
200 FOR K=3 TO N
210 C(K)=B(K-1)-Q*C(K-2)-P*C(K-1)
220 NEXT K
230 W=C(N-1)*C(N-1)-C(N)*C(N-2)
240 U=(B(N-1)*C(N-1)-B(N)*C(N-2))/W
250 V=(B(N)*C(N-1)-B(N-1)*C(N))/W
260 IF ABS(U)+ABS(V)<1E-5 THEN GOSUB 400:GOTO 280
270 P=P+U:Q=Q+V:GOTO 140
280 N=N-2
290 FOR K=1 TO N
300 A(K)=B(K)
310 NEXT K
320 GOTO 100
400 D=P*P/4-Q:PRINT
410 D0=SQR(ABS(D))
420 PRINT "X(";M-N+1;")=";
430 IF D<0 THEN PRINT -P/2;"+J*";D0:PRINT:PRINT "X(";M-N+2;")=";-P/2;"-J*";D0
440 IF D>=0 THEN PRINT (D0-P/2):PRINT:PRINT "X(";M-N+2;")=";(-D0-P/2)
450 RETURN