Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- line = input("","s");
- he=strsplit(line);
- n=length(he);
- mat = ones (n, n);
- for i = 1:n
- mat(1,i) = str2num(he{1,i});
- endfor
- for i = 2:n
- for i2 = 1:n
- mat(i,i2) = scanf("%d","C");
- endfor
- endfor
- printf("%.6f\n",det(mat))
- ev=sort(eig(mat));
- if(length(ev)==0)
- exit()
- endif
- rev=[ev(1)];
- for i = 2:length(ev)
- if(ev(i)-ev(i-1)>0.001)
- rev(end+1)=ev(i);
- endif
- endfor
- for i = length(rev):-1:1
- printf("%.6f\n",rev(i));
- temp=mat;
- for i2 = 1:n
- temp(i2,i2)-=rev(i);
- endfor
- ans=rref(null(temp)');
- for i2 = 1:rows(ans)
- for i3 = 1:columns(ans)
- printf("%.6f ",ans(i2,i3));
- endfor
- printf("\n");
- endfor
- endfor
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement