Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- clc
- clear all
- format short e
- M=magic(60);
- I=eye(60);
- C=M+I;
- A=C(:,1:10);
- b=sum(A,2);
- [Q,R]=qr(A);
- xqr=R\(Q'*b);
- x=ones(10,1);
- norm(xqr-x)/norm(xqr)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement