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

| アクセス |

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

グラフィックギャラリー

差分方程式の部屋3

しめ縄リング

d := 0.1 ; p := 0.1

m := 10000

{X[0], Y[0]} = {5, -5}

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.7, 0.3, 0]]

{5, -5}

[Graphics:HTMLFiles/gallery13_8.gif]

-Graphics -

ジグソ - パズル

d := 0.001 ; p := 0.2

m := 10000

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

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.5, 1]]

{1, -5}

[Graphics:HTMLFiles/gallery13_17.gif]

-Graphics -

牛

d := 0.001 ; p := 0.05

m := 6000

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

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]

{10, 10}

[Graphics:HTMLFiles/gallery13_26.gif]

-Graphics -

怪文字

d := 0.0001 ; p := 0.03

m := 10000

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

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]

{5, 5}

[Graphics:HTMLFiles/gallery13_35.gif]

-Graphics -

クッション

d := 0.02 ; p := 0.002

m := 10000

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

Do[{X[n + 1] = -Y[n] + p Sin[X[n]] + 3Sin[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, 0.7]]

{3, 3}

[Graphics:HTMLFiles/gallery13_44.gif]

-Graphics -

十字

d := 0.001 ; p := 0.01

m := 10000

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

Do[{X[n + 1] = -Y[n] + p Sin[X[n]] + 0.55Cos[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.5, 0]]

{3, 3}

[Graphics:HTMLFiles/gallery13_53.gif]

-Graphics -

輝き

a := 0.1 ; p := 200

m := 10000

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

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.5, 0, 1]]

{1, 1}

[Graphics:HTMLFiles/gallery13_62.gif]

-Graphics -

星A

a := 0.001 ; p := 10000

m := 5000

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

Do[{X[n + 1] = a X[n] - Y[n] + p 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[1, 0.3, 0]]

{1, 0}

[Graphics:HTMLFiles/gallery13_71.gif]

-Graphics -

星B

a := 0.03 ; p := 30

m := 5000

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

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

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

{10, 10}

[Graphics:HTMLFiles/gallery13_80.gif]

-Graphics -

的

a := 0.01 ; p := 1000

m := 5000

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

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

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

{2, 2}

[Graphics:HTMLFiles/gallery13_89.gif]

-Graphics -

歯車

a := 0.1 ; p := 3

m := 5000

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

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

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

{10, 10}

[Graphics:HTMLFiles/gallery13_98.gif]

-Graphics -

スト - ンサ - クル

a := 0.1 ; p := 3.9

m := 5000

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

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

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

{5, 5}

[Graphics:HTMLFiles/gallery13_107.gif]

-Graphics -

あぶく

a := 0.002 ; p := 12000

m := 30000

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

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

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

{5, 0}

[Graphics:HTMLFiles/gallery13_116.gif]

-Graphics -


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