diff --git a/BLAS/SRC/caxpy.f b/BLAS/SRC/caxpy.f
index 6503a80c1..6879a14bc 100644
--- a/BLAS/SRC/caxpy.f
+++ b/BLAS/SRC/caxpy.f
@@ -103,13 +103,16 @@
 *
 *     .. Local Scalars ..
       INTEGER I,IX,IY
+      COMPLEX CDUM
 *     ..
-*     .. External Functions ..
-      REAL SCABS1
-      EXTERNAL SCABS1
+*     .. Statement Functions ..
+      REAL CABS1
+*     ..
+*     .. Statement Function definitions ..
+      CABS1(CDUM) = ABS(REAL(CDUM)) + ABS(AIMAG(CDUM))
 *     ..
       IF (N.LE.0) RETURN
-      IF (SCABS1(CA).EQ.0.0E+0) RETURN
+      IF (CABS1(CA).EQ.0.0E+0) RETURN
       IF (INCX.EQ.1 .AND. INCY.EQ.1) THEN
 *
 *        code for both increments equal to 1
diff --git a/BLAS/SRC/zaxpy.f b/BLAS/SRC/zaxpy.f
index f9081e957..c46465554 100644
--- a/BLAS/SRC/zaxpy.f
+++ b/BLAS/SRC/zaxpy.f
@@ -103,13 +103,16 @@
 *
 *     .. Local Scalars ..
       INTEGER I,IX,IY
+      COMPLEX*16 ZDUM
 *     ..
-*     .. External Functions ..
-      DOUBLE PRECISION DCABS1
-      EXTERNAL DCABS1
+*     .. Statement Functions ..
+      DOUBLE PRECISION CABS1
+*     ..
+*     .. Statement Function definitions ..
+      CABS1(ZDUM) = ABS(DBLE(ZDUM)) + ABS(DIMAG(ZDUM))
 *     ..
       IF (N.LE.0) RETURN
-      IF (DCABS1(ZA).EQ.0.0d0) RETURN
+      IF (CABS1(ZA).EQ.0.0D+0) RETURN
       IF (INCX.EQ.1 .AND. INCY.EQ.1) THEN
 *
 *        code for both increments equal to 1