Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- r=var('r')
- M=Matrix(SR,16,16,[[-2,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0],
- [1,-3,0,1,1,0,0,0,0,0,0,0,0,0,0,0],
- [0,1,-3,1,0,1,0,0,0,0,0,0,0,0,0,0],
- [1,0,0,-3,0,0,1,1,0,0,0,0,0,0,0,0],
- [0,1,0,0,-3,0,1,0,1,0,0,0,0,0,0,0],
- [0,0,0,1,0,-3,1,0,0,1,0,0,0,0,0,0],
- [3,0,0,0,0,0,-3,0,0,0,6/(3-2*r),0,0,0,0,0],
- [0,3,0,0,0,0,0,-3,0,0,3/(3-2*r),3/(3-2*r),0,0,0,0],
- [0,0,0,3,0,0,0,0,-3,0,3/(3-2*r),0,3/(3-2*r),0,0,0],
- [0,0,0,0,0,0,3,0,0,-3,3/(3-2*r),0,0,3/(3-2*r),0,0],
- [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],
- [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],
- [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],
- [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],
- [0,0,0,0,0,0,0,0,0,0,3*(2-r)/(3-2*r),0,0,0,-3,1/(3-r)],
- [0,0,0,0,0,0,0,0,0,0,3*(1-r)*(3-r)/(3-2*r),0,0,0,0,-3],
- ])
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement