use PDL; $n = 50; $vec = zeroes ($n+1); index($vec,0) .= 1;