FUNCTION ERF(X) C C BÛ×ÈCËEHÈE ÔÓHKÖÈÈ ERF(X) C AËÃOPÈTM C X < 1 - PßÄ TEÉËOPA C X < 4 - AÏÏPOKCÈMAÖÈß ÏOËÈHOMOM C X > 4 - ERF= 1. C REAL C(5) DATA P/0.3275911/ DATA C/0.254829592,-0.284496736,1.421413741,-1.453152027 *,1.061405429/ C AX=X IF (AX.LT.0) AX=-X ERF=1. IF (AX.LT.1) GO TO 40 IF (AX.LT.4) GO TO 41 C 43 IF (X.LT.0) ERF=-ERF RETURN C C PßÄ TEÉËOPA 40 A=AX I=0 K=1 ERF=AX 42 K=K+2 I=I+1 A=-A*AX*AX/I A1=A/K ERF=ERF+A1 IF (ABS(A1).GE.1E-5) GOTO 42 ERF=ERF*1.12838 GO TO 43 C C AÏÏPOKCÈMAÖÈß 41 ERF=0. A=1./(1.+P*AX) A1=1. DO 44 I=1,5 A1=A*A1 44 ERF=ERF+C(I)*A1 ERF=1.-ERF*EXP(-AX*AX) GOTO 43 END