python:自分の考えがあっているか確認したいです。iとｊに関するfor構文

Question

前提・実現したいこと

a[str(i)]はiの値が大きくなればその値も大きくなる増加関数のような形をしています。
ここで、a[str(i)]/a[str(i-1)]<1.001 ー※
となるとき、その条件※を満たすa[str(i)]をm=a[str(i)]としました。
ただし、a[str(i)]の中にはh[j]という数列が含まれており、
各h[j]に対するm[j]の値を求めたいです。
例えば、
※を満たすのは
h=0の時、m[j=0]=a[str(i=50)]
h=0.01の時m[j=1]=a[str(i=52)]
h=0.02の時m[j=2]=a[str(i=55)]
……
といった具合です。

最終的に
(h,m)=(h[1],m[1]),(h[2],m[2])…(h[100],m[100])して、
hとmの関数にしてプロットしたいです。

i,jに対するfor構文を用いて画像にプロットしてみた
のですが、自分の考えが正しいか確認していただきたいです。

該当のソースコード

a = {'0': 0.1}
m = a['0'],
h=np.arange(0,1,0.01)
for j in range(0,99):
    for i in range(1, 100):
        a[str(i)] = np.tanh(a[str(i - 1)]+h[j])
        if a[str(i)]/a[str(i-1)]<1.001:
            break
    m = m + ( a[ str(i) ], )
plt.plot(h,m)
plt.show()

補足情報（FW/ツールのバージョンなど）

ここにより詳細な情報を記載してください。

Answer 1

レビューの丸投げになっています。
プログラムが意図通り動いてるのかどうか、を確認するのもプログラマの責任です。
テストしましょう。

ここで全く関係ない第三者が「あってる」「間違ってる」と言ったところで鵜呑みにしていいモノでもないです。

Comments

[kunikuni96]

通知を付けておらず、返答が遅れて申し訳ありません。gentaroさんのおっしゃる通りです。自分で見直してやりたかった数値計算を別の方法で見直してみた結果どうやらこのプログラムで正しいようです。自分で考えるきっかけになりました。ありがとうございました。