10 '•û’öŽ®‚Ì�ª(NEWTON–@)
20 PRINT "f(x(1),x(2))=0 ‚ð�@Line900 ‚É�@“ü‚ê‚ĉº‚³‚¢�B"
25 DIM V(2),A(2),W(2),G(2),F(2),X(2),C(2),D(2),L(2),B(2),Y(2)
30 INPUT "Žn“_1=";V(1)
40 INPUT "”÷�¬‹æŠÔEPS1=";A(1)
50 INPUT "‹æŠÔ•�1=";W(1)
51 INPUT "Žn“_2=";V(2)
52 INPUT "”÷�¬‹æŠÔeps2=";A(2)
53 INPUT "‹æŠÔ•�2=";W(2)
55 K=0:L(1)=0:L(2)=0
56 X(2)=V(2)
57 OPEN "NEWTON.DAT" for create as #2
60 G(1)=V(1):F(1)=V(1):Z=0
65 G(2)=V(2):F(2)=V(2)
70 IF Z=0 THEN GOTO 90
80 G(1)=G(1)-W(1):C(1)=G(1):'GOTO 100
85 G(2)=G(2)-W(2):C(2)=G(2):GOTO 100
90 C(1)=G(1):Z=1
95 C(2)=G(2)
100 GOSUB 700
110 F(1)=F(1)+W(1):C(1)=F(1)
115 F(2)=F(2)+W(2):C(2)=F(2)
120 GOSUB 700
125 K=K+1
130 IF K<50 then GOTO 70
135 CLOSE #2
140 END
700 '
704 L(1)=0
706 X(1)=C(1):'GOSUB 900
707 X(2)=C(2):GOSUB 900
710 Y(1)=B(1):X(1)=A(1)+C(1)
715 Y(2)=B(2):X(2)=A(2)+C(2)
720 GOSUB 900
730 D(1)=C(1):C(1)=D(1)-A(1)*Y(1)/(B(1)-Y(1))
735 D(2)=C(2):C(2)=D(2)-A(2)*Y(2)/(B(2)-Y(2))
740 IF L(1)>30 THEN RETURN:'GOTO 704 
745 :ELSEIF (ABS(D(1)-C(1)))>0.001 then L(1)=L(1)+1:GOTO 706 ENDIF
746 IF (ABS(D(2)-C(2)))>0.001 then L(1)=L(1)+1:GOTO 706 ENDIF
750 PRINT #2,"X(1)=";C(1)
755 PRINT #2,"X(2)=";C(2)
756 L(1)=L(1)+1
760 return
900 'B=((X-2)*X-1)*X+2
901 'B=X+LOG(X):REM ŠÖ�”Ž® f(x)=0
902 'B=X^4+2*X^3-2*X^2+3*X-4:' a0*x^0+a1*x^1+...+an*x^n=0
903 'B=(X(1)-3)^2+(X(2)+1)^2
904 'B(1)=-900*X(1)+0.026*LOG(2)+0.026*LOG((X(2)*8800/33200-0.7)/5900/X(1))
905 'B(2)=X(2)*8800/33200-X(1)*11800-0.7
906 'B=-900*(X(2)*(8400+X(1))/33200/11800-0.7/11800)+0.026*LOG(2*118/59)
907 'B(1)=X(1)+LOG(X(1))
908 'B(2)=X(2)+LOG(X(2))
909 B(1)=-11*10^3*X(2)+0.026*LOG((2/X(2)*(X(1)-0.7)*0.5-0.7)/100/10^3)
910 B(2)=(X(1)-0.7)*0.5-X(2)*100*10^3-0.7
930 RETURN