2015年8月1日土曜日

150801

Numpy


多項式の展開

昨日行った多項式の展開を(次数の順に気をつけながら)Numpyでやってみた。

# -*- coding: cp932 -*-
import numpy

f0 = numpy.poly1d([1, 1])
f1 = [3, 2, 1]
print f0 * f1 # 片方がリストでも答えが出る
print f1 * f0 # 片方がリストでも答えが出る
print f0 ** 3
print list(reversed((f0 ** 3).c))

出力結果
   3     2
3 x + 5 x + 3 x + 1
   3     2
3 x + 5 x + 3 x + 1
   3     2
1 x + 3 x + 3 x + 1
[1, 3, 3, 1]

0 件のコメント:

コメントを投稿

注: コメントを投稿できるのは、このブログのメンバーだけです。