Pythonの繰り返し演算処理を行いたい

Question

実現したいこと

dydx=x+c
(xは１から１０までの整数値をとる変数、cは定数である)
という微分方程式を定数c=1,2,3,4,5の場合それぞれ解きたい.

具体的には
dydx=x+1
dydx=x+2
dydx=x+3
dydx=x+4
dydx=x+5
の５つの微分方程式をfor文を用いて順番に解いていきそれぞれのyをリスト型の
c=1の場合はypoints_1
c=2の場合はypoints_2
c=3の場合はypoints_3
c=4の場合はypoints_4
c=5の場合はypoints_5
に代入したい

• ▲▲機能を動作するようにする

前提

ここに質問の内容を詳しく書いてください。
上記の内容を解く

発生している問題・エラーメッセージ

NameError                                 Traceback (most recent call last)
<ipython-input-18-9cb7432c946d> in <cell line: 18>()
     18 for c in crange:
     19     for x in xrange:
---> 20         ypoints_c.append(dydx)
     21         k_1=func(x,a,b)
     22         dydx +=k_1

NameError: name 'ypoints_c' is not defined

該当のソースコード

Python

import numpy as np
import matplotlib.pyplot as plt

xrange=np.arange(1,11,1)

c=0
crange=np.arange(1,6,1)
dydx=0
ypoints_1=[]
ypoints_2=[]
ypoints_3=[]
ypoints_4=[]
ypoints_5=[]

def func(x,a,b):
    dydx=x+c
    return dydx
for c in crange:
    for x in xrange:
        ypoints_c.append(dydx)
        k_1=func(x,a,b)
        dydx +=k_1

試したこと

リストypointsの名前を変えてみたりした

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

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

Answer 1

ypoints を numpy.ndarray にするとよいかと思います。

xrange = np.arange(1, 11, 1)
crange = np.arange(1, 6, 1)
ypoints = np.empty((len(crange), len(xrange)))

for c in crange:
    for x in xrange:
        ypoints[c-1, x-1] = dydx
        k_1 = func(x, a, b)
        dydx += k_1

print(ypoints)

# [[  0.   2.   5.   9.  14.  20.  27.  35.  44.  54.]
#  [ 65.  68.  72.  77.  83.  90.  98. 107. 117. 128.]
#  [140. 144. 149. 155. 162. 170. 179. 189. 200. 212.]
#  [225. 230. 236. 243. 251. 260. 270. 281. 293. 306.]
#  [320. 326. 333. 341. 350. 360. 371. 383. 396. 410.]]

Answer 2

ypoints_1のような別個のリスト変数を用意するのではなくcの値をキーとするypointsという辞書変数を利用してはどうでしょうか。
なお、以下の例示コードでは微分方程式うんぬんはおいといて、変数の使い方のみを示しています

import numpy as np

xrange=np.arange(1,11,1)
crange=np.arange(1,6,1)
dydx=0

ypoints = {}

for c in crange:
    ypoints[c] = []
    for x in xrange:
        ypoints[c].append(dydx)
        #k_1=func(x,a,b)
        k_1 = 1
        dydx +=k_1

print(ypoints)
# {1: [0, 1, 2, 3, 4, 5, 6, 7, 8, 9], 2: [10, 11, 12, 13, 14, 15, 16, 17, 18, 19], 3: [20, 21, 22, 23, 24, 25, 26, 27, 28, 29], 4: [30, 31, 32, 33, 34, 35, 36, 37, 38, 39], 5: [40, 41, 42, 43, 44, 45, 46, 47, 48, 49]}

Comments

[east-site.3]

丁寧なご回答ありがとうございました.出力なのですが５つの微分方程式を個別に解いた値をypointsに代入したいのです.
具体的には微分方程式

dydx=x+1
dydx=x+2
dydx=x+3
dydx=x+4
dydx=x+5

を積分した式のyを以下のように５つのリストまたは辞書変数に代入したいのです

y=(1/2)*x**2+x. x=1から10までを代入した値 　1.5 4 7.5 12. 17.5 24. 31.5 40. 49.5
y=(1/2)*x**2+2*x x=1から10までを代入した値　2.5 6. 10.5 16. 22.5 30. 38.5 48. 58.5
y=(1/2)*x**2+3*x x=1から10までを代入した値　3.5 8. 13.5 20. 27.5 36. 45.5 56. 67.5
y=(1/2)*x**2+4*x x=1から10までを代入した値　4.5 10. 16.5 24. 32.5 42. 52.5 64. 76.5
y=(1/2)*x**2+5*x x=1から10までを代入した値　5.5 12. 19.5 28. 37.5 48. 59.5 72. 85.5

[8524ba23]

x=1～10の順で、(1/2)*x**2+c*x の値をリストに追加すればよいかと思いますが。

[east-site.3]

度重なるご返信ありがとうございます．ご返答頂いたプログラムにそれを組み込もうと色々取り組んでみたのですがどうもうまく行きません．もしよろしければご回答頂いたプログラムに
dydx=x+1を積分したy=(1/2)*x**2+x. の結果 　1.5 4 7.5 12. 17.5 24. 31.5 40. 49.5
dydx=x+2を積分したy=(1/2)*x**2+2*x 　2.5 6. 10.5 16. 22.5 30. 38.5 48. 58.5
dydx=x+3を積分したy=(1/2)*x**2+3*x 　3.5 8. 13.5 20. 27.5 36. 45.5 56. 67.5
dydx=x+4を積分したy=(1/2)*x**2+4*x 　4.5 10. 16.5 24. 32.5 42. 52.5 64. 76.5
dydx=x+5を積分したy=(1/2)*x**2+5*x 　5.5 12. 19.5 28. 37.5 48. 59.5 72. 85.5
を出力するプログラムをお見せ頂いてもよろしいでしょうか．