Advertisement
yojimbos_law

dumb linear system (n=3,k=6)

Mar 13th, 2015
307
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.81 KB | None | 0 0
  1. r=var('r')
  2. M=Matrix(SR,16,16,[[-2,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0],
  3. [1,-3,0,1,1,0,0,0,0,0,0,0,0,0,0,0],
  4. [0,1,-3,1,0,1,0,0,0,0,0,0,0,0,0,0],
  5. [1,0,0,-3,0,0,1,1,0,0,0,0,0,0,0,0],
  6. [0,1,0,0,-3,0,1,0,1,0,0,0,0,0,0,0],
  7. [0,0,0,1,0,-3,1,0,0,1,0,0,0,0,0,0],
  8. [3,0,0,0,0,0,-3,0,0,0,6/(3-2*r),0,0,0,0,0],
  9. [0,3,0,0,0,0,0,-3,0,0,3/(3-2*r),3/(3-2*r),0,0,0,0],
  10. [0,0,0,3,0,0,0,0,-3,0,3/(3-2*r),0,3/(3-2*r),0,0,0],
  11. [0,0,0,0,0,0,3,0,0,-3,3/(3-2*r),0,0,3/(3-2*r),0,0],
  12. [1,0,0,0,0,0,0,0,0,0,3*(2-r)/(3-2*r),3*(1-r)/(3-2*r),0,0,0,0],
  13. [0,1,0,0,0,0,0,0,0,0,3*(1-r)/(3-2*r),-3,3*(1-r)/(3-2*r),0,0,0],
  14. [0,0,0,1,0,0,0,0,0,0,3*(1-r)/(3-2*r),0,-3,3*(1-r)/(3-2*r),0,0],
  15. [0,0,0,0,0,0,1,0,0,0,3*(1-r)/(3-2*r),0,0,-3,3*(1-r)/(3-2*r),0],
  16. [0,0,0,0,0,0,0,0,0,0,3*(2-r)/(3-2*r),0,0,0,-3,1/(3-r)],
  17. [0,0,0,0,0,0,0,0,0,0,3*(1-r)*(3-r)/(3-2*r),0,0,0,0,-3],
  18. ])
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement