REAL FUNCTION CHGSIN(A)
	CHGSIN=-A
	RETURN
	END
	REAL FUNCTION FALOG(A)
	INTEGER PRT,DIV,LN,SQR
	COMMON /JOBLOG/ LOG
	COMMON /MATHER/ PRT,DIV,LN,SQR
	IF (A .LE. 0) GOTO 10
	FALOG=ALOG(A)
	RETURN
10	LN=LN+1
	IF (PRT .EQ. 0.0) WRITE(LOG,100) A
100	FORMAT(' *******LOG OF NON-POSITIVE NUMBER*******',G10.4)
	FALOG=0.0
	RETURN
	END
	REAL FUNCTION FAND(A,B)
	FAND=B
	IF (B .NE. 0.0) FAND=A
	RETURN
	END
	REAL FUNCTION FEQUAL(A,B)
	FEQUAL=0.0
	IF (A .EQ. B) FEQUAL=1.0
	RETURN
	END
	REAL FUNCTION FAREP(A,B)
	A=B
	FAREP=B
	RETURN
	END
	REAL FUNCTION FGREAT(A,B)
	FGREAT=0.0
	IF (A .GT. B) FGREAT=1.0
	RETURN
	END
	REAL FUNCTION FGREQ(A,B)
	FGREQ=0.0
	IF (A .GE. B) FGREQ=1.0
	RETURN
	END
	REAL FUNCTION FLESEQ(A,B)
	FLESEQ=0.0
	IF (A .LE. B) FLESEQ=1.0
	RETURN
	END
	REAL FUNCTION FLESS(A,B)
	FLESS=0.0
	IF (A .LE. B) FLESS=1.0
	RETURN
	END
	REAL FUNCTION FNOTEQ(A,B)
	FNOTEQ=0.0
	IF (A .NE. B) FNOTEQ=1.0
	RETURN
	END
	REAL FUNCTION FOR(A,B)
	FOR=A
	IF (A .EQ. 0) FOR=B
	RETURN
	END
	REAL FUNCTION FSQRT(A)
	INTEGER PRT,DIV,LN,SQR
	COMMON /JOBLOG/ LOG
	COMMON /MATHER/ PRT,DIV,LN,SQR
	IF (A .LT. 0) GOTO 10
	FSQRT=SQRT(A)
	RETURN
10	SQR=SQR+1
	IF (PRT .EQ. 0) WRITE(LOG,100)A
100	FORMAT(' *******SQUARE ROOT OF NEGATIVE NUMBER*******',G10.4)
	FSQRT=SQRT(-A)
	RETURN
	END
	REAL FUNCTION FXTOY(X,Y)
	FXTOY=X**Y
	RETURN
	END
	REAL FUNCTION FNOT(A)
	IF (A .EQ. 0.0) GOTO 10
	FNOT=0.0
	RETURN
10	FNOT=1.0
	RETURN
	END
	REAL FUNCTION PDIV(A,B)
	INTEGER PRT,DIV,LN,SQR
	COMMON /JOBLOG/ LOG
	COMMON /MATHER/ PRT,DIV,LN,SQR
	IF (B .EQ. 0.0) GOTO 10
	PDIV=A/B
	RETURN
10	DIV=DIV+1
	IF (PRT .EQ. 0) WRITE(LOG,100)A,B
100	FORMAT(' *******REAL DIVISION BY ZERO*******',G10.4,'/',G10.4)
	PDIV=0.0
	RETURN
	END
	REAL FUNCTION PADD(A,B)
	PADD=A+B
	RETURN
	END
	REAL FUNCTION PSUB(A,B)
	PSUB=A-B
	RETURN
	END
	REAL FUNCTION PMUL(A,B)
	PMUL=A*B
	RETURN
	END
	REAL FUNCTION CONF01(X)
	REAL CONF(20)
	DATA CONF/63.657,9.925,5.841,4.604,4.032,3.707,3.499,
	1 3.355,3.250,3.169,3.106,3.055,3.012,2.977,2.947,2.921
	2,2.898,2.878,2.861,2.845/
	I=X-0.5
	IF (I .GE. 1) GOTO 10
	CONF01=1E38
	RETURN
10	IF (I .GT. 20) GOTO 20
	CONF01=CONF(I)
	RETURN
20	CONF01=2.576
	RETURN
	END
	REAL FUNCTION CONF05(X)
	REAL CONF(20)
	DATA CONF/12.706,4.303,3.182,2.776,2.571,2.447,2.365,
	1 2.306,2.262,2.228,2.201,2.179,2.160,2.145,2.131,2.120,
	2 2.110,2.101,2.093,2.086/
	I=X-0.5
	IF (I .GE. 1) GOTO 10
	CONF05=1E38
	RETURN
10	IF (I .GT. 20) GOTO 20
	CONF05=CONF(I)
	RETURN
20	CONF05=1.96
	RETURN
	END