#include <math.h>

#include <stdio.h>

#include <gsl/gsl_matrix.h>

#include <gsl/gsl_blas.h>

int

main (void)

{

  size_t i,j;

  gsl_matrix *m = gsl_matrix_alloc (10, 10);

  for (i = 0; i < 10; i++)

    for (j = 0; j < 10; j++)

      gsl_matrix_set (m, i, j, sin (i) + cos (j));

  for (j = 0; j < 10; j++)

    {

      gsl_vector_view column = gsl_matrix_column (m, j);

      double d;

      d = gsl_blas_dnrm2 (&column.vector);

      printf ("matrix column %zu, norm = %g\n", j, d);

    }

  gsl_matrix_free (m);

  return 0;

}