```/* nag_regsn_quant_linear_iid (g02qfc) Example Program.
*
* NAGPRODCODE Version.
*
* Copyright 2016 Numerical Algorithms Group.
*
* Mark 26, 2016.
*/
/* Pre-processor includes */
#include <stdio.h>
#include <nag.h>
#include <nag_stdlib.h>
#include <nagg02.h>

#define B(i,j) b[(j) * m + i]
#define BU(i,j) bu[(j) * m + i]
#define BL(i,j) bl[(j) * m + i]
#define X(i,j) x[(i) * m + j]

int main(void)
{
/* Integer scalar and array declarations */
Integer i, j, l, m, n, ntau;
Integer *info = 0;
Integer exit_status = 0;

/* NAG structures */
NagError fail;

/* Double scalar and array declarations */
double df;
double *b = 0, *bl = 0, *bu = 0, *tau = 0, *x = 0, *y = 0;

/* Initialize the error structure */
INIT_FAIL(fail);

printf("nag_regsn_quant_linear_iid (g02qfc) Example Program Results\n\n");

/* Skip heading in data file */
scanf("%*[^\n] ");

/* Read in the problem size */
scanf("%" NAG_IFMT "%" NAG_IFMT "%" NAG_IFMT "%*[^\n] ", &n, &m, &ntau);

/* Allocate memory for input arrays */
if (!(y = NAG_ALLOC(n, double)) ||
!(tau = NAG_ALLOC(ntau, double)) || !(x = NAG_ALLOC(n * m, double)))
{
printf("Allocation failure\n");
exit_status = -1;
goto END;
}

/* Read in the data */
for (i = 0; i < n; i++) {
for (j = 0; j < m; j++)
scanf("%lf", &X(i, j));
scanf("%lf", &y[i]);
}
scanf("%*[^\n] ");

/* Read in the quantiles required */
for (l = 0; l < ntau; l++) {
scanf("%lf", &tau[l]);
}
scanf("%*[^\n] ");

/* Allocate memory for output arrays */
if (!(b = NAG_ALLOC(m * ntau, double)) ||
!(info = NAG_ALLOC(ntau, Integer)) ||
!(bl = NAG_ALLOC(m * ntau, double)) ||
!(bu = NAG_ALLOC(m * ntau, double)))
{
printf("Allocation failure\n");
exit_status = -1;
goto END;
}

/* nag_regsn_quant_linear_iid (g02qfc).Quantile linear regression, simple
interface, independent, identically distributed (IID) errors */
nag_regsn_quant_linear_iid(n, m, x, y, ntau, tau, &df, b, bl, bu, info,
&fail);
if (fail.code != NE_NOERROR) {
printf("Error from nag_regsn_quant_linear_iid (g02qfc).\n%s\n",
fail.message);
if (fail.code == NW_POTENTIAL_PROBLEM) {
for (i = 0; i < ntau; i++)
printf("%" NAG_IFMT " ", info[i]);
printf("\n");
}
else {
exit_status = 1;
goto END;
}
}

/* Display the parameter estimates */
for (l = 0; l < ntau; l++) {
printf(" Quantile: %6.3f\n\n", tau[l]);
printf("         Lower   Parameter   Upper\n");
printf("         Limit   Estimate    Limit\n");
for (j = 0; j < m; j++) {
printf(" %3" NAG_IFMT "   %7.3f   %7.3f   %7.3f\n", j + 1, BL(j, l),
B(j, l), BU(j, l));
}
printf("\n\n");
}

END:

NAG_FREE(info);
NAG_FREE(b);
NAG_FREE(bl);
NAG_FREE(bu);
NAG_FREE(tau);
NAG_FREE(x);
NAG_FREE(y);

return (exit_status);
}
```