# Clean environment and load libraries
rm(list = ls())  # Clear environment to avoid conflicts


# Source files
library(MixStable)

n <- if (identical(Sys.getenv("NOT_CRAN"), "true")) 1200 else 100
X1 <- rstable(n, 1.2, 0.5, 1, 3, pm = 1)


# 🔍 Estimate parameters using MLE
estimated_params <- Max_vrai(X1)
print("📊 MLE Results:")
print(estimated_params)