test_that("test-MOGIW works", {
  x <- rMOGIW(10,1,3,2,3)
  fx<- dMOGIW(x,1,3,2,3)
  fx2 <- dMOGIW(x,alpha=1,beta=3,lambda=2,theta=3)
  expect_equal(fx, fx2)

})