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


# Source files
library(MixStable)

# 🎲 Generate synthetic data
X1 <- rstable(1200, 1.2, 0.5, 1, 3, pm=1)

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