python ループ zip 等差数列

Question

ループ処理で困っています。
現在、下記のようなプログラムを実装しました。しかしこちらをループ処理で簡単にしたいと
思いここ何日か悩んでいます。ネットを見ても具体例があまりなく参考になりませんでした。よければ教えてください・
fig = plt.figure()
ax = fig.add_subplot(1, 1, 1)
ax.set_aspect('equal')
ax.grid()

p1 = (0, 0)
p2 = (0.0768, -0.133)
p3 = (0.1536, 0)
p4= (0.63, 0)
p5= (0.7068,-0.133)
p6= (0.7836,0)
p7= (1.26,0)
p8= (1.3368,-0.133)
p9= (1.4136,0)

ax.plot(*zip(p1, p2, p3, p4,p5,p6,p7,p8,p9,))
plt.xlim([0,7])
plt.ylim([-0.2,0.1])
plt.show()

－－－－－－－－－－－－－－－－－
x=0.63
z=-0.133
L=0.1536

p1 = (0, 0)
p2 = (L/2, -0.133)
p3 = (L, 0)
p4= (X, 0)
p5= (X+L/2,-0.133)
p6= (X+L,0)
p7= (2X,0)
p8= (2X+L/2,-0.133)
p9= (2X+L,0)

(p1,p4,p7),(p2,p5,p8)(p3,p6,p9)が等差数列です.

Comments

[tiitoi]

p_i = (x_i, y_i) を表す漸化式を考えてみてください

[11390]

p_iとはどう意味ですか？

Answer 1

p_nは2つの要素で成り立っているように見えます。
このそれぞれの要素について、漸化式(前のいくつかの項から次の項を計算する)や一般項の式(nからその項を計算する)を作ってみると良いと思います。
今回であれば、p_n = (x_n, y_n)とすると、(n=0,1,2...)
x_n=0.63*(n//3)+0.0768*(n%3)
で、
y_n=0 (このときn%3=0,2)
y_n=-0.133 (このときn%3=1)
のように見えます。
この式をコードにしてみてください。if文で分岐させるのでも、//や%の演算子を使うのでも実装できると思います。
ここで、割り算a/b=cあまりdとしたとき、a//b=cかつa%b=dです。

Comments

[11390]

追加なのですが、現在初期値がｐ１（０．０）となっていますが、これを任意の数にしたうえでループさせるにはどうすればよいでしょうか？

[退会済みユーザー]

初項を(a,b)として、
x_n=a+0.63*(n//3)+0.0768*(n%3)
y_n=b (このときn%3=0,2)
y_n=b-0.133 (このときn%3=1)
などとして目的の挙動になればこうすれば良いと思います。
目的の挙動にならないときは、漸化式をきちんと立ててやると良いと思います。漸化式とは、数列{a}について、a_(n+1)をa_nで表す式のことです。