pure function zones_1(n) result(A) integer, intent(in) :: n complex(8) :: A(n) A = one end function zones_1