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

| アクセス |

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

グラフィックギャラリー

差分方程式の部屋1

竜巻

h := 0.05 ; m := 1000

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

Do[{X[n + 1] = X[n] + h Sin[Y[n]^2], Y[n + 1] = Y[n] + h X[n]}, {n, 0, m}]

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

{1, 0}

[Graphics:HTMLFiles/gallery4_7.gif]

- Graphics -

富士山

h := 0.3 ; m := 2000

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

Do[{X[n + 1] = X[n] + h Sin[Y[n]]/(1 + Y[n]^2), Y[n + 1] = Y[n] + h X[n]}, {n, 0, m}]

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

{1, 100}

[Graphics:HTMLFiles/gallery4_15.gif]

- Graphics -

たまねぎ

h := 0.1 ; m := 2000

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

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

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

{0, 0.5}

[Graphics:HTMLFiles/gallery4_23.gif]

- Graphics -

魔球

h := 0.1 ; m := 1000

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

Do[{X[n + 1] = X[n] + h Sin[Y[n]^3], Y[n + 1] = Y[n] + h Sin[X[n] + Y[n]]}, {n, 0, m}]

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

{1, 2}

[Graphics:HTMLFiles/gallery4_31.gif]

- Graphics -

乱歩

h := 0.1 ; m := 30000

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

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

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

{0, 1}

[Graphics:HTMLFiles/gallery4_39.gif]

- Graphics -

渦

a := 0.01 ; d := 0.1 ; p := 0.03

m := 5000

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

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

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

{1, 0}

[Graphics:HTMLFiles/gallery4_48.gif]

- Graphics -

鏡の部屋

a := 0.4 ; p := 200

m := 20000

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

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

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

{1, 0}

[Graphics:HTMLFiles/gallery4_57.gif]

- Graphics -

まんじ

a := 0.0003 ; p := 15000

m := 4500

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

Do[{X[n + 1] = a X[n] - Y[n] + p Sin[X[n]], 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/gallery4_66.gif]

- Graphics -

十字架

a := 0.02 ; p := 20

m := 5000

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

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

{1, 0}

[Graphics:HTMLFiles/gallery4_75.gif]

- Graphics -

六角星

a := 0.999 ; p := 300

m := 20000

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

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

{1, 0}

[Graphics:HTMLFiles/gallery4_84.gif]

- Graphics -

ファスナ -

d := 0.1 ; p := 0.2

m := 5000

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

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

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

{1, 0}

[Graphics:HTMLFiles/gallery4b_8.gif]

- Graphics -


Created by Mathematica  (July 21, 2006) Valid XHTML 1.1!
このページのトップへ戻る