gsl_matrix2k17/2k18Ikolok termin3

#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()
{
    gsl_matrix *m;
    gsl_vector *v1;
    int i,j,K,L,N,Z;
    printf("Vnesete gi dimenziite na matricata KxL:\t");
    scanf("%d%d",&K,&L);
    m=gsl_matrix_calloc(K,L);
    printf("Vnesete ja dimenzijata na vektorot N:");
    scanf("%d",&N);
    v1=gsl_vector_alloc(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++)
    {
        double u=gsl_rng_uniform(r);
        int izbor=(int)(u*1000)%145+25;
        gsl_vector_set(v1,i,izbor);

    }
    gsl_rng_free(r);
    printf("Vektorot na pocetok izgleda vaka: \n");
    for(i=0;i<N;i++)
    {
        printf("%g",gsl_vector_get(v1,i));
    }
    int k,br=1;
    printf("Vnesete go Z: ");
    scanf("%d",&Z);
    for(i=0;i<K;i++)
    {
        for(j=0;j<L;j++)
        {
            gsl_vector *pom;
            pom=gsl_vector_alloc(Z);
            int brojac=0;
         for(k=(br*Z)-Z;k<br*Z;k++)
         {

            gsl_vector_set(pom,brojach,gsl_vector_get(v1,k));
            brojac++;
         }
         gsl_matrix_set(m,i,j,gsl_vector_max(pom));
         br++;
         gsl_vector_free(pom);

        }
    }
    printf("Na kraj matricata izgleda vaka: \n");
    for(i=0;i<K;i++)
    {
        for(j=0;j<L;j++)
        {
            printf("%g",gsl_matrix_get(m,i,j));
        }
    }
    printf("\n");
    gsl_matrix_free(m);
    gsl_vector_free(v1);
    return 0;
}