```/* nag_univar_robust_1var_trimmed (g07ddc) Example Program.
*
* Copyright 2019 Numerical Algorithms Group.
*
* Mark 27.0, 2019.
*
*/

#include <nag.h>
#include <stdio.h>

#define NMAX 1000
int main(void)
{

/* Local variables */
Integer exit_status = 0, i, k, n;
NagError fail;
double alpha, propn, *sx = 0, tmean, tvar, wmean, wvar, *x = 0;

INIT_FAIL(fail);

printf("nag_univar_robust_1var_trimmed (g07ddc) Example Program Results\n\n");
/* Skip heading in data file */
scanf("%*[^\n] ");
scanf("%" NAG_IFMT " ", &n);
if (n >= 2) {
if (!(x = NAG_ALLOC(NMAX, double)) || !(sx = NAG_ALLOC(NMAX, double)))
{
printf("Allocation failure\n");
exit_status = -1;
goto END;
}
}
else {
printf("Invalid n.\n");
exit_status = 1;
return exit_status;
}
for (i = 1; i <= n; ++i)
scanf("%lf ", &x[i - 1]);
scanf("%lf ", &alpha);

/* nag_univar_robust_1var_trimmed (g07ddc).
* Trimmed and winsorized mean of a sample with estimates of
* the variances of the two means
*/
nag_univar_robust_1var_trimmed(n, x, alpha, &tmean, &wmean, &tvar, &wvar, &k, sx,
&fail);
if (fail.code != NE_NOERROR) {
printf("Error from nag_univar_robust_1var_trimmed (g07ddc).\n%s\n",
fail.message);
exit_status = 1;
goto END;
}

propn = (double) k / n;
propn = 100.0 - propn * 200.0;
printf("Statistics from middle %6.2f%% of data\n\n", propn);
printf("               Trimmed-mean = %11.4f\n", tmean);
printf("   Variance of Trimmed-mean = %11.4f\n\n", tvar);
printf("            Winsorized-mean = %11.4f\n", wmean);
printf("Variance of Winsorized-mean = %11.4f\n", wvar);
END:
NAG_FREE(x);
NAG_FREE(sx);
return exit_status;
}
```