pure function zzeros_2(n1,n2) result(A) integer, intent(in) :: n1,n2 complex(8) :: A(n1,n2) A = zero end function zzeros_2