Clear[A, x]
A = {{2, -1, 0}, {0, 2, -1}, {0, -1, 2}};
MatrixForm[%]
vx = {0.5, 0.5, 1.1};
Do[vy = A.vx;
Print[k, " ", vy, " ", vy[[1]]/vx[[1]], " ", vy[[2]]/vx[[2]]];
vx = vy/Max[Abs[vy]], {k, 1, 15}]
Eigensystem[A]
精品开放课程共享系统 - 爱课程 - 首页
http://www.icourses.cn/web/sword/portal/shareDetails?cId=6302#/course/chapter
网友评论