SUBROUTINE NEWTON(X,Y,GRAD)
C BERECHNET DAS INTERPOLATIONS-POLYNOM MAXIMAL 10.GRADES
C X(GRAD+2) IST DIE INTERPOLATIONSSTELLE - Y(GRAD+2) DAS ERGEBNIS
C FORTRAN II 
	DIMENSION X(12),Y(12),C(10)
	K=1
1004	S=0.
	I=1
	P=1.
1013	IF((K-1)-I)1001,1007,1007
1007	P=P*(X(K+1)-X(I))
	S=S+C(I)*P
	I=I+1
	GOTO 1013
1001	P=P*(X(K+1)-X(K))
	C(K)=(Y(K+1)-Y(I)-S)/P
	IF(K-IGRAD)1003,1006,1003
1003	K=K+1
	GOTO 1004
1006	S=Y(I)
	P=1.
	K=1
1002	P=P*(X(IGRAD+2)-X(K))
	S=S+C(K)*P
	IF(K-IGRAD)1017,1009,1017
1009	Y(IGRAD+2)=S
	RETURN
1017	K=K+1
	GOTO 1002
	END