gsl_matrix 2k16/2k17 Ikolok termin1

#include<stdio.h>
#include<stdlib.h>
#include<gsl/gsl_matrix.h>
#include<gsl/gsl_vector.h>
#include<gsl/gsl_rng.h>
#include<time.h>

int main()
{
    int N,i,j,k,l;
    gsl_matrix *m;
    gsl_vector *v1;
    printf("Vnesi ja dolzinata na vektorot N:\t");
    scanf("%d",&N);
    v1=gsl_vector_alloc(N);
    for(k=0;k<N;k++)
    {
        printf("Vnesete gi elementite na vektorot: \n");
        int broj;
        scanf("%d",&broj);
        gsl_vector_set(v1,k,broj);
    }
    m=gsl_matrix_calloc(N,N);
    gsl_rng *r;
    gsl_rng_env_setup();
    r=gsl_rng_alloc(gsl_rng_default());
    gsl_rng_set(r,time(0));
    for(i=0;i<N;i++)
    {
        for(j=0;j<N;j++)
        {
            for(k=0;k<N;k++)
            int br=0;
            if(i==j)
                {
                    gsl_matrix_set(m,i,j,gsl_vector_get(v1,k+br));
                }
            for(l=0;l<N;l++)
                {
                    if(i==j && l>i && l>j)
                        {
                            double u=gsl_rng_uniform(r);
                            int broj=(int)(u*1000)+(gsl_vector_get(v1,k));
                            gsl_matrix_set(m,i,l,broj);
                        }

                }
                            if(i>j)
                                {
                                    break;
                                }
        }
    }
gsl_rng_free(r);
    printf("Novata matrica e: \n");
      for(i=0;i<N;i++)
    {
        for(j=0;j<N;j++)
        {
            printf("%g",gsl_matrix_get(m,i,j));
        }
    }
    gsl_matrix_free(m);
    gsl_matrix_free(v1);
    return 0;
}