test_that("multiplication works", {
  expect_equal(2 * 2, 4)
})

test_that("sim_students returns binary matrix of correct dimensions", {
  alpha <- sim_students(N = 10, K = 4)
  expect_true(is.matrix(alpha))
  expect_equal(dim(alpha), c(10, 4))
  expect_true(all(alpha %in% c(0, 1)))
})