#include <stdio.h>
#include <gsl/gsl_rng.h>
#include <gsl/gsl_randist.h>
#include <gsl/gsl_statistics.h>
#include <gsl/gsl_sort.h>
#include <math.h>

int main(){
    gsl_rng * r;
    r = gsl_rng_alloc (gsl_rng_mt19937);
    gsl_rng_default_seed=2344;
    int i, n;
    double mean, variance, median, skew, kurt,q1,q3;
    double alpha, beta, mu, sigma;
    double data[1000000];
    n=1000000;
    for(i=0;i<n;i++){
      data[i]=gsl_ran_landau(r);
      /*      printf("%lf \n",data[i]);*/
    }
    gsl_sort(data,1,n);
    mean=gsl_stats_mean(data,1,n);
    variance=gsl_stats_variance(data,1,n);
    median=gsl_stats_quantile_from_sorted_data(data,1,n,0.5);
    q1=gsl_stats_quantile_from_sorted_data(data,1,n,0.25);
    q3=gsl_stats_quantile_from_sorted_data(data,1,n,0.75);
    printf("Mean = %lf\n",mean);
    printf("SD = %lf\n",sqrt(variance));
    printf("Q1 = %lf\n",q1);
    printf("Q2 = %lf\n",median);
    printf("Q3 = %lf\n",q3);
    
    return(0);
}