そういや日曜に手計算*1の前にMathematicaで計算していたので貼っておこう． Mathematicaでは，例えばSeries[f[x,y,z],{x,0,6},{y,0,6},{z,0,6}]でfを各文字につき0を中心にして6次まで級数展開できる．しかしこれではx^6 y^6 z^6とかの項が含まれてしまう．…

a**b**cはNonCommutativeMultiply[a,b,c]の略記であり，またNonCommutativeMultiplyの属性Flatのためにa**(b**c)等とも同値である．NonCommutativeMultiplyには属性OneIdentityも設定されているため，NonCommutativeMultiplyの引数をパタンマッチで取り出す…

DSolve[{x''[t] == -Sign[ x[t] ], x[0] == 1, x'[0] == 0}, x, t] が返ってこない件について．Reduceに変えたりAssumingで に制限したりしても駄目．がっかりさせないでくれよ... NDSolve（数値解法）を使ってPlotすると放物線らしきものが出て，√2くらいに…

Mathematicaで単に Table[k, {k, 1000000}] を評価させるとガリガリ言って戻ってこないけど Total[Table[k, {k, 1000000}]] とか Take[Table[k, {k, 1000000}]] は一瞬で答えが返ってくる．つまり単純にeagerな評価をしているのではない．まぁあれだけ成熟し…