http://comp.chem.tohoku.ac.jp/mediawiki/index.php?action=history&feed=atom&title=Polynomial2_struct 2026-05-26T19:55:06Z このウィキのこのページに関する変更履歴 MediaWiki 1.36.2 http://comp.chem.tohoku.ac.jp/mediawiki/index.php?title=Polynomial2_struct&diff=1108&oldid=prev 2026-05-26T04:14:49Z

<p>ページの作成:「多項式の計算を行う構造体。 和・差・積、微分・積分、シフト・反転 等の操作が行える。 == サンプルコード == &lt;source lang=&quot;for…」</p> <p><b>新規ページ</b></p><div>多項式の計算を行う構造体。 和・差・積、微分・積分、シフト・反転 等の操作が行える。<br /> <br /> == サンプルコード ==<br /> <br /> &lt;source lang=&quot;fortran&quot;&gt;<br /> program main<br /> use polynomial2_module,only: polynomial2_struct<br /> use polynomial2_module,only: polynomial2 !!! コンストラクタ<br /> implicit none<br /> type(polynomial2_struct) :: p1, p2, p3<br /> REAL8 :: c1(0:2), c2(0:3)<br /> <br /> <br /> !!! 多項式の係数の設定 !!!<br /> c1(0:2) = [ 1.0d0, 2.0d0, 1.0d0 ] !!! 1 + 2 x + x**2<br /> c2(0:3) = [ 3.0d0, 2.0d0, 0.0d0, 2.0d0 ] !!! 3 + 2 x + + 2 x**3<br /> <br /> !!! 構造体の作成 !!!<br /> p1 = polynomial2(c1)<br /> p2 = polynomial2(c2)<br /> p3 = 3.0d0 !!! 定数関数は直接代入も可<br /> <br /> write(*,*)<br /> <br /> !!! 値の出力 !!!<br /> <br /> !!! 多項式の次数<br /> write(*,*) p1%degree !!! 2<br /> <br /> !!! x を指定したときの値<br /> write(*,*) p1%substitute(1.0d0) !!! 4<br /> <br /> !!! 文字列として出力<br /> write(*,*) p1%string() !!! 1 + 2 x + x**2<br /> <br /> <br /> write(*,*)<br /> write(*,*) &quot;###########################&quot;<br /> write(*,*)<br /> <br /> <br /> !!! 和・差・積、定数倍、代入 !!!<br /> <br /> p3 = p1 + p2<br /> write(*,*) p3%string() !!! 4 + 4 x + x**2 + 2 x**3<br /> <br /> p3 = p1 - p2<br /> write(*,*) p3%string() !!! - 2 + x**2 - 2 x**3<br /> <br /> p3 = p1 * p2<br /> write(*,*) p3%string() !!! 3 + 8 x + 7 x**2 + 4 x**3 + 4 x**4 + 2 x**5<br /> <br /> !!! 定数は右からでも左からでも掛けることができる。<br /> p3 = p1 * 2.0d0<br /> write(*,*) p3%string() !!! 2 + 4 x + 2 x**2<br /> p3 = 2.0d0 * p1<br /> write(*,*) p3%string() !!! 2 + 4 x + 2 x**2<br /> <br /> !!! 多項式の代入も可<br /> p3 = p1%substitute(p2)<br /> write(*,*) p3%string() !!! 16 + 16 x + 4 x**2 + 16 x**3 + 8 x**4 + 4 x**6<br /> <br /> <br /> write(*,*)<br /> write(*,*) &quot;###########################&quot;<br /> write(*,*)<br /> <br /> <br /> !!! 微分・積分 !!!<br /> <br /> !!! 微分<br /> p3 = p1%differentiate()<br /> write(*,*) p3%string() !!! 2 + 2 x<br /> <br /> !!! 積分<br /> !!! x = 0 のとき y = 0 になるように積分定数が決められる。<br /> p3 = p1%integrate()<br /> write(*,*) p3%string() !!! x + x**2 + (1/3) x**3<br /> <br /> <br /> write(*,*)<br /> write(*,*) &quot;###########################&quot;<br /> write(*,*)<br /> <br /> <br /> !!! シフト・反転操作 !!!<br /> <br /> !!! 多項式のグラフを右方向に 1.0d0 だけシフトさせる。<br /> !!! 負の値を指定すれば左にシフトする。<br /> p3 = p1%x_shift(1.0d0)<br /> write(*,*) p3%string() !!! x**2<br /> <br /> !!! 多項式のグラフを上方向に 1.0d0 だけシフトさせる。<br /> !!! 負の値を指定すれば下にシフトする。<br /> p3 = p1%y_shift(1.0d0)<br /> write(*,*) p3%string() !!! 2 + 2 x + x**2<br /> <br /> !!! 上下反転<br /> p3 = p1%x_reflect()<br /> write(*,*) p3%string() !!! - 1 - 2 x - x**2<br /> <br /> !!! 左右反転<br /> p3 = p1%y_reflect()<br /> write(*,*) p3%string() !!! 1 - 2 x + x**2<br /> <br /> !!! 点対称に反転<br /> p3 = p1%point_reflect()<br /> write(*,*) p3%string() !!! - 1 + 2 x - x**2<br /> <br /> <br /> write(*,*)<br /> write(*,*) &quot;###########################&quot;<br /> write(*,*)<br /> <br /> <br /> endprogram<br /> &lt;/source&gt;</div>

Hirano