File: POLY2.FT of Tape: Various/ETH/eth11-2

(Source file text) 

---

C	POLY2.FTN - SAMPLE POLYNOMIAL OPERATIONS USING:
C  PADD		- ADD TWO POLYNOMIALS
C  PCLD		- LINEAR SYNTHETIC DIVISION
C  PGCD		- GCD OF TWO POLYNOMIALS
C	CALLS PNORM, PDIV
C  PVSUB	- COMPOSITION OF POLYNOMIALS
C	CALLS PMPY,PADDM,PCLA
C  PSUB		- SUBTRACTION OF POLYNOMIALS
C  PILD		- EVALUATE POLYNOMIAL AND ITS FIRST DERIVATIVE
C	CALLS PQSD
C
	DIMENSION P1(4),P2(3),P3(3),P6(4),W1(5),W2(5)
	DATA P1/-72.,-9.,2.,1./,P2/26.,-2.,6./,P3/6.,-8.,8./
C
C	OUTPUT CHANNEL = IOUT
	IOUT=2
	WRITE(IOUT,100)
100	FORMAT(//' ADDING THE FOLLOWING POLYNOMIALS:'/)
	WRITE(IOUT,101) P1(4),P1(3),P1(2),P1(1)
101	FORMAT(10X,'(',F5.1,')X**3 +(',F5.1,')X**2 +(',
	1F5.1,')X +(',F5.1,')'/)
	WRITE(IOUT,102) P2(3),P2(2),P2(1)
102	FORMAT(23X,'(',F5.1,')X**2 +(',F5.1,')X +(',F5.1,')'/)
	CALL PADD(P1,IP1,P1,4,P2,3)
	WRITE(IOUT,103)
103	FORMAT(' SUM IS:'/)
	WRITE(IOUT,101) P1(4),P1(3),P1(2),P1(1)
	WRITE(IOUT,104)
104	FORMAT(//' SUBTRACT THE FOLLOWING POLYNOMIAL:'/)
	WRITE(IOUT,102) P3(3),P3(2),P3(1)
	CALL PSUB(P1,IP1,P1,4,P3,3)
	WRITE(IOUT,105)
105	FORMAT(/' RESULT IS:'/)
	WRITE(IOUT,101) P1(4),P1(3),P1(2),P1(1)
C
C	COMPOSITION
	P2(1)=2.
	P2(2)=1.
	WRITE(IOUT,106)
106	FORMAT(/' NOW PRECOMPOSE WITH X+2:'/)
	CALL PVSUB(P6,IP6,P1,4,P2,2,W1,W2)
	WRITE(IOUT,105)
	WRITE(IOUT,101) P6(4),P6(3),P6(2),P6(1)
	WRITE(IOUT,107)
107	FORMAT(/' TRANSLATING INTO FACTORS OF X-2:'/)
	CALL PCLD(P6,4,2.)
	WRITE(IOUT,105)
	WRITE(IOUT,101) P6(4),P6(3),P6(2),P6(1)
	P1(1)=45.
	P1(2)=12.
	P1(3)=1.
	WRITE(IOUT,108)
108	FORMAT(//' G.C.D. OF ABOVE POLYNOMIAL AND X**2+12X+45 IS:'/)
	IP1=3
	CALL PGCD(P6,4,P1,IP1,W1,1.E-2,IER)
	WRITE(IOUT,105)
	WRITE(IOUT,102) P1(3),P1(2),P1(1)
C
C	EVALUATE POLYNOMIAL AND DERIVATIVE
	CALL PILD(PV,DV,-6.,P1,3)
	WRITE(IOUT,109) PV,DV
109	FORMAT(/' THE FUNCTION AT X=-6 EQUALS'F8.3/
	1' ITS DERIVATIVE AN X=-6 EQUALS'F8.3)
	STOP
	END

---