Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #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,M,K,i,j,broj;
- gsl_matrix *m,*m2;
- gsl_vector *v;
- printf("Vnesete vrednost za N i M:\t");
- scanf("%d %d",&N,&M);
- m=gsl_matrix_alloc(N,M);
- printf("Vnesete go K");
- scanf("%d",&K);
- v=gsl_vector_alloc(K);
- for(i=0;i<N;i++);
- {
- for(j=0;j<M;j++)
- {
- printf("Vnesete gi broevite vo matricata: \n");
- scanf("%d",&broj);
- gsl_matrix_set(m,i,j,broj);
- }
- }
- printf("Pocetnata matrica e:\n");
- for(i=0;i<N;i++)
- {
- for(j=0;j<M;j++)
- {
- printf("%g",gsl_matrix_get(m,i,j));
- }
- }
- int x,y,k,l;
- for(i=0;i<K;i++)
- {
- printf("Vnesi i,j par za odreduvanje na pozicija: \n");
- scanf("%d%d",&x,&y);
- if(x==N-1||y==M-1)
- {
- gsl_vector_set(v,i,-1);
- }
- else
- {
- gsl_matrix *pom;
- pom=gsl_matrix_alloc(2,2);
- for(k=0;k<2;k++)
- {
- for(l=0;l<2;l++)
- {
- gsl_matrix_set(pom,k,l,gsl_matrix_get(m,x+k,y+l));
- }
- }
- gsl_vector_set(v,i,gsl_matrix_max(pom));
- gsl_matrix_free(pom);
- }
- }
- printf("Krajniot vektor izgleda vaka:\n");
- for(i=0;i<K;i++)
- {
- printf("%g",gsl_vector_get(v,i));
- }
- gsl_matrix_free(m);
- gsl_vector_free(v);
- return 0;
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement