I wrote these R codes which to help me compute some statistical properties of the probability distribution I am working on. Please, find the R codes below;

library('numDeriv')

library('Matrix')

beta=1.5

alpha=1.6

lambda=1.5

theta=0.6

#Define the quantile of exponentiated WeiGeo distribution

quantile=function(beta,alpha,lambda,theta,u){

  result