10 REM **********   ΠΡΟΓΡΑΜΜΑ Π206   **********

20 REM *                                       *

30 REM *        Μέθοδος Newton-Raphson        *

40 REM *       Για πολυωνυμικές εξισώσεις      *

50 REM *                                       *

60 REM *****************************************

70 CLS

80 EPS=.0000005

90 DIM A(20),B(20),C(20)

100 INPUT "ΔΩΣΕ ΤΟ ΒΑΘΜΌ ΤΟΥ ΠΟΛΥΩΝΎΜΟΥ";N:CLS

110 NN=N

120 ND=0

130 FOR I=0 TO N

140 READ A(I)

150 NEXT I

160 PRINT

170 X=1

180 FOR I=1 TO 100

190 GOSUB 360

200 IF C(N-1)<EPS THEN GOTO 470

210 XN=X-B(N)/C(N-1)

220 IF ABS(XN-X)<EPS THEN GOTO 290

230 X=XN

240 ON ERROR GOTO 260

250 NEXT I

260 IF ND=0 THEN GOTO  450

270 PRINT "Δεν υπάρχει άλλη πραγματική ρίζα"

280 END

290 PRINT "H X(";ND;") είναι  = ";X;

300 PRINT " Επαναλήψεις";I:PRINT

310 ND=ND+1

320 IF ND=NN THEN GOTO 270

330 N=N-1

340 FOR L=0 TO N:A(L)=B(L):NEXT L

350 GOTO 170

360 B(0)=A(0)

370 FOR J=1 TO N

380 B(J)=B(J-1)*X+A(J)

390 NEXT J

400 C(0)=B(0)

410 FOR J=1 TO N-1

420 C(J)=C(J-1)*X+B(J)

430 NEXT J

440 RETURN

450 PRINT "Δεν υπάρχει καμιά πραγματική ρίζα"

460 END

470 IF I=1 THEN GOTO 490

480 IF I>1 THEN GOTO 260

490 X=X+1

500 GOTO 180

510 DATA 1,-2,-3,8,-4