線形代数



  • 線形代数の計算にはパッケージの読み込みが必要です。
    それをしないと以下のようになります。

    > a:=vector([2,sin(x),4,5.3,beta]);
    > A:=matrix([[1,x,y],[0,1,z],[0,0,1]]);



  • パッケージを読み込みます。

    > with(linalg);



  • ベクトルを読み込みます。

    > v:=vector([2,sin(omega),cos(omega),Pi,alpha]);



  • マトリックスを読み込みます。

    > A:=matrix([[1,alpha,beta],[0,1,gamma],[0,0,1]]);



  • もう一つベクトルを読み込みます。

    > u:=vector([1,cos(omega),cos(omega),2.5,beta]);



  • ベクトルを加算します。

    > u+v;



  • この場合 evalm による評価が必要です。

    > evalm(u+v);



  • 関数 add を使えば evalm は必要ありません。

    > add(u,v);



  • 以下、いくつかの例です。

    > 2*u+5*v;



    > evalm(2*u+5*v);



    > add(u,v,2,5);



  • ベクトルの乗算は以下の方法ではできません。

    > u*v;



    > evalm(u*v);