弘前大学 理工学部 数物科学科 数理科学コース

| アクセス |

トップページグラフィックギャラリー > 差分方程式の部屋2

グラフィックギャラリー

差分方程式の部屋2

つむじ風

d := 0.1 ; p := 0.1

m := 5000

{X[0], Y[0]} = {1, 0}

Do[{X[n + 1] = -Y[n] + p Sin[X[n] + n], Y[n + 1] = X[n] + d Y[n]}, {n, 0, m}]

ListPlot[Table[{X[n], Y[n]}, {n, 0, m}], AspectRatio→Automatic, PlotStyle→RGBColor[1, 0, 1]]

{1, 0}

[Graphics:HTMLFiles/gallery9_8.gif]

-Graphics -

ヒトデ

d := 0.5 ; p := 0.1

m := 10000

{X[0], Y[0]} = {1, 0}

Do[{X[n + 1] = -Y[n] + p Sin[X[n] + n], Y[n + 1] = X[n] + d Y[n]}, {n, 0, m}]

ListPlot[Table[{X[n], Y[n]}, {n, 0, m}], AspectRatio→Automatic, PlotStyle→RGBColor[0, 0, 1]]

{1, 0}

[Graphics:HTMLFiles/gallery9_17.gif]

-Graphics -

団子

d := 0.25 ; p := 0.05

m := 8000

{X[0], Y[0]} = {1, 0}

Do[{X[n + 1] = -Y[n] + p Sin[X[n]] + 0.005n, Y[n + 1] = X[n] + d Y[n]}, {n, 0, m}]

ListPlot[Table[{X[n], Y[n]}, {n, 0, m}], AspectRatio→Automatic, PlotStyle→RGBColor[0, 1, 0]]

{1, 0}

[Graphics:HTMLFiles/gallery9_26.gif]

-Graphics -

ボ - リング

d := 0.1 ; p := 0.04

m := 2500

{X[0], Y[0]} = {1, 0}

Do[{X[n + 1] = -Y[n] + p Sin[X[n]] + 0.005n, Y[n + 1] = X[n] + d Y[n]}, {n, 0, m}]

ListPlot[Table[{X[n], Y[n]}, {n, 0, m}], AspectRatio→Automatic, PlotStyle→RGBColor[0, 0.5, 1]]

{1, 0}

[Graphics:HTMLFiles/gallery9_35.gif]

-Graphics -

六角星

d := 0.9 ; p := 0.2

m := 2000

{X[0], Y[0]} = {1, 0}

Do[{X[n + 1] = -Y[n] + p Sin[X[n]] + Sin[n], Y[n + 1] = X[n] + d Y[n]}, {n, 0, m}]

ListPlot[Table[{X[n], Y[n]}, {n, 0, m}], AspectRatio→Automatic, PlotStyle→RGBColor[1, 0, 1]]

{1, 0}

[Graphics:HTMLFiles/gallery9_44.gif]

-Graphics -

4つのリング

d := 0.01 ; p := 0.1

m := 10000

{X[0], Y[0]} = {1, 0}

Do[{X[n + 1] = -Y[n] + p Sin[X[n]] + 3Cos[n], Y[n + 1] = X[n] + d Y[n]}, {n, 0, m}]

ListPlot[Table[{X[n], Y[n]}, {n, 0, m}], AspectRatio→Automatic, PlotStyle→RGBColor[1, 0, 0]]

{1, 0}

[Graphics:HTMLFiles/gallery9_53.gif]

-Graphics -

4つのリング

p := 0.05

m := 10000

{X[0], Y[0]} = {1, 0}

Do[{X[n + 1] = -Y[n] + p Sin[X[n]], Y[n + 1] = X[n] + 3Cos[n]}, {n, 0, m}]

ListPlot[Table[{X[n], Y[n]}, {n, 0, m}], AspectRatio→Automatic, PlotStyle→RGBColor[0, 0.5, 1]]

{1, 0}

[Graphics:HTMLFiles/gallery9_62.gif]

-Graphics -

四方

a := 0.01 ; p := 2

m := 20000

{X[0], Y[0]} = {1, 2}

Do[{X[n + 1] = a X[n] - Y[n] + p n/(1 + X[n]^2), Y[n + 1] = X[n]}, {n, 0, m}]

ListPlot[Table[{X[n], Y[n]}, {n, 0, m}], AspectRatio→Automatic, PlotStyle→RGBColor[0, 0.5, 0]]

{1, 2}

[Graphics:HTMLFiles/gallery9_71.gif]

-Graphics -

五方

a := 0.6 ; p := 300

m := 2000

{X[0], Y[0]} = {1, 0}

Do[{X[n + 1] = a X[n] - Y[n] + p n X[n]/(1 + X[n]^2), Y[n + 1] = X[n]}, {n, 0, m}]

ListPlot[Table[{X[n], Y[n]}, {n, 0, m}], AspectRatio→Automatic, PlotStyle→RGBColor[0, 0, 1]]

{1, 0}

[Graphics:HTMLFiles/gallery9_80.gif]

-Graphics -

六方

a := 0.9 ; p := 3000

m := 2000

{X[0], Y[0]} = {1, 0}

Do[{X[n + 1] = a X[n] - Y[n] + p n X[n]/(1 + X[n]^2), Y[n + 1] = X[n]}, {n, 0, m}]

ListPlot[Table[{X[n], Y[n]}, {n, 0, m}], AspectRatio→Automatic, PlotStyle→RGBColor[0.5, 0.5, 0]]

{1, 0}

[Graphics:HTMLFiles/gallery9_89.gif]

-Graphics -

座布団

a := 0.01 ; p := 20000

m := 15000

{X[0], Y[0]} = {1, 0}

Do[{X[n + 1] = a X[n] - Y[n] + p ArcTan[X[n] + 0.6n], Y[n + 1] = X[n]}, {n, 0, m}]

ListPlot[Table[{X[n], Y[n]}, {n, 0, m}], AspectRatio→Automatic, PlotStyle→RGBColor[1, 0, 1]]

{1, 0}

[Graphics:HTMLFiles/gallery9_98.gif]

-Graphics -


Created by Mathematica  (April 1, 2009) Valid XHTML 1.1!
このページのトップへ戻る