#include <gsl/gsl_vector.h>
#include <gsl/gsl_cblas.h>
#include <math.h>

int main(int argc, char **argv) {
	gsl_vector *vec = gsl_vector_alloc(10);
	double norm = 0.0;
	int i;
	
	for(i=0; i<10; i++) {
		gsl_vector_set(vec, i, (double)i);
	}
	
	norm = gsl_blas_dnrm2(vec);
	printf("norm1 : %f\n", norm);
	
	gsl_blas_ddot(vec, vec, &norm);
	norm = sqrt(norm);
	printf("norm2 : %f\n", norm);
	
	return(0);
}