sieving

still need solving the matrix