John Cagnol


Paris, France

Introduction au Calcul Scientifique
ESI CS 1402 - Printemps 2010

Exercice I de la feuille d'exercices 3

> with(LinearAlgebra);
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
> M:=1/3*Matrix([[4,0,3,0,-3],[-1+3*alpha,3*alpha,-3,4-3*alpha,0],[0,0,1,3,-3],[3,0,0,4,-3],[3,0,3,3,-5]]);
*
> ?JordanForm
> J:=JordanForm(M);
*
> Q:=JordanForm(M,output='Q');
*
> P:=MatrixInverse(Q);
*
> MatrixInverse(P) . J . P;
*
> map(evalf,J)^100;
*