```/* nag_summary_stats_freq (g01adc) Example Program.
*
* Copyright 2014 Numerical Algorithms Group.
*
* Mark 7, 2001.
*/

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

int main(void)
{

/* Scalars */
double   xsd, xskew, xkurt, xmean;
Integer  exit_status = 0, i, j, k, kmin1, n, nprob;

NagError fail;

/* Arrays */
double   *x = 0;
Integer  *ifreq = 0;

INIT_FAIL(fail);

/* Skip heading in data file */
scanf("%*[^\n] ");
scanf("%ld%*[^\n] ", &nprob);
for (j = 1; j <= nprob; ++j)
{
scanf("%ld%*[^\n] ", &kmin1);
k = kmin1 + 1;

/* Allocate memory */
if (!(x = NAG_ALLOC(k, double)) ||
!(ifreq = NAG_ALLOC(k, Integer)))
{
printf("Allocation failure\n");
exit_status = -1;
goto END;
}

for (i = 1; i <= kmin1; ++i)
scanf("%lf%ld", &x[i - 1], &ifreq[i - 1]);
scanf("%lf%*[^\n] ", &x[k - 1]);

printf("\nProblem %4ld\n", j);
printf("Number of classes %4ld\n", kmin1);

* Mean, variance, skewness, kurtosis, etc., one variable,
* from frequency table
*/
nag_summary_stats_freq(k, x, ifreq, &xmean, &xsd, &xskew, &xkurt, &n,
&fail);

if (fail.code == NE_NOERROR)
{
printf("Successful call of "
printf("          Class          Frequency\n\n");
for (i = 1; i <= kmin1; ++i)
printf("%10.2f%10.2f%12ld\n", x[i-1], x[i],
ifreq[i-1]);

printf("\n Mean %16.4f\n", xmean);
printf(" Std devn%13.4f\n", xsd);
printf(" Skewness%13.4f\n", xskew);
printf(" Kurtosis%13.4f\n", xkurt);
printf(" Number of cases%8ld\n", n);
}
else
{