MATEXE:  PROC(A,B,C);

/* PROCEDURE DOES AN ELEMENT BY ELEMENT MULTIPLICATION ON MATRIX
   ARGUMENTS TWO AND THREE RETURNING MATRIX ARGUMENT ONE. 

   AUTHOR:  AL JAWORSKI
   DATE:    26 APRIL 1982 */

DCL (A(*,*),B(*,*),C(*,*)) FLOAT BIN(53);
DCL (I,J,AL1,AH1,AL2,AH2) FIXED BINARY(31);

AL1=LBOUND(A,1);
AH1=HBOUND(A,1);
AL2=LBOUND(A,2);
AH2=HBOUND(A,2);
IF AL1^=LBOUND(B,1) |
   AL1^=LBOUND(C,1) |
   AH1^=HBOUND(B,1) |
   AH1^=HBOUND(C,1) |
   AL2^=LBOUND(B,2) |
   AL2^=LBOUND(C,2) |
   AH2^=HBOUND(B,2) |
   AH2^=HBOUND(C,2) THEN DO;
  PUT EDIT('Error--Incompatible matrix dimensions for multiplication')
          (SKIP,A);
  STOP;
END;

DO I=AL1 TO AH1;
  DO J=AL2 TO AH2;
    A(I,J)=B(I,J)*C(I,J);
  END;
END;

END MATEXE;