next up previous
Next: 2 数式処理とは Up: 情報処理II 資料 Mathematica 入門 Previous: 情報処理II 資料 Mathematica 入門

1 Mathematica ってこんなもの

以下の例は、 数学科の計算機である oyabun にログインして、 Mathematica を実行してみたものです (Windows 環境での実行例ではありませんが、 本質は同じです)。 ここではプログラムなどは書かずに、 式を順次入力して計算結果を表示させています。

\includegraphics[width=12cm]{eps/torus.ps}

--------------------------------------------

oyabun% math
Mathematica 4.0 for Solaris
Copyright 1988-1999 Wolfram Research, Inc.
 -- Motif graphics initialized -- 

In[1]:= 1/2 + 1/3                   ← 分数計算

        5
Out[1]= -                           → ちょっと見難いですけどね
        6

In[2]:= a={{0,1},{6,1}}             ← 行列の入力

Out[2]= {{0, 1}, {6, 1}}

In[3]:= Eigenvalues[a]              ← 行列の固有値の計算

Out[3]= {-2, 3}

In[4]:= Eigenvectors[a]             ← 行列の固有ベクトルの計算

Out[4]= {{-1, 2}, {1, 3}}

In[5]:= Expand[(x+y)^6]             ← 式の展開

         6      5         4  2       3  3       2  4        5    6
Out[5]= x  + 6 x  y + 15 x  y  + 20 x  y  + 15 x  y  + 6 x y  + y

In[6]:= N[Pi,50]                    ← 円周率 50 桁

Out[6]= 3.1415926535897932384626433832795028841971693993751

In[7]:= Integrate[Log[x],x]         ← 不定積分

Out[7]= -x + x Log[x]

In[8]:= Plot3D[x^2 - y^2, {x,-1,1}, {y,-1,1}]  ← グラフ

Out[8]= -Graphics-                  → ここで画面に図が表示されます

In[9]:= Solve[x^3+2x==1,x]          ← 3 次方程式を解かせてみる

                                            9 + Sqrt[177] 1/3
                                           (-------------)
                           2         1/3          2
Out[1]= {{x -> -2 (-----------------)    + ------------------}, 
                   3 (9 + Sqrt[177])               2/3
                                                  3
 
                                    2         1/3
>    {x -> (1 + I Sqrt[3]) (-----------------)    - 
                            3 (9 + Sqrt[177])
 
                         9 + Sqrt[177] 1/3
        (1 - I Sqrt[3]) (-------------)
                               2
>       ----------------------------------}, 
                         2/3
                      2 3
 
                                    2         1/3
>    {x -> (1 - I Sqrt[3]) (-----------------)    - 
                            3 (9 + Sqrt[177])
 
                         9 + Sqrt[177] 1/3
        (1 + I Sqrt[3]) (-------------)
                               2
>       ----------------------------------}}
                         2/3
                      2 3

(ちょっとびっくりするでしょうか。
カルダノの公式というやつで、知らないと訳が分からないかもしれないけれど。)

In[10]:= ParametricPlot3D[{Cos[t](3+Cos[u]),Sin[t](3+Cos[u]),Sin[u]},
         {t,0,2Pi},{u,0,2Pi}]       ← トーラスを描かせる。

Out[10]:= -Graphics3D-

In[11]:= Quit                       ← 終了
oyabun% 


next up previous
Next: 2 数式処理とは Up: 情報処理II 資料 Mathematica 入門 Previous: 情報処理II 資料 Mathematica 入門
Masashi Katsurada
平成20年10月18日