library(MixStable)
x <- seq(0.1, 2, length.out = 5)
alpha <- 1.5; beta <- 0.5; delta <- 1.0; omega <- 1.0

res_g     <- integrate_sine_weighted(x, alpha, beta, delta, omega)
res_h_tld <- integrate_sine_log_weighted(x, alpha, beta, delta, omega)
res_h     <- integrate_cosine_log_weighted(x, alpha, beta, delta, omega)

cat("Result sine_weighted_exp_ralpha:\n")
print(res_g)
cat("Result sine_log_weighted_exp_ralpha:\n")
print(res_h_tld)
cat("Result cosine_log_weighted_exp_ralpha:\n")
print(res_h)