質問編集履歴
4
訂正線の追加
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -172,15 +172,13 @@
|
|
|
172
172
|
|
|
173
173
|
|
|
174
174
|
|
|
175
|
-
~~free_en(400)
|
|
175
|
+
~~free_en(400)~~
|
|
176
|
-
|
|
176
|
+
|
|
177
|
-
plt.plot(xli,free, label="t=400")
|
|
177
|
+
~~plt.plot(xli,free, label="t=400")~~
|
|
178
|
-
|
|
178
|
+
|
|
179
|
-
plt.legend()
|
|
179
|
+
~~plt.legend()~~
|
|
180
|
-
|
|
181
|
-
|
|
182
|
-
|
|
180
|
+
|
|
183
|
-
↑これを500,600,700,,1100についてやったら目標のグラフができた。~~
|
|
181
|
+
~~↑これを500,600,700,,1100についてやったら目標のグラフができた。~~~~
|
|
184
182
|
|
|
185
183
|
|
|
186
184
|
|
|
@@ -196,7 +194,7 @@
|
|
|
196
194
|
|
|
197
195
|
|
|
198
196
|
|
|
199
|
-
頂いた回答に
|
|
197
|
+
頂いた回答に凡例を加えた。
|
|
200
198
|
|
|
201
199
|
|
|
202
200
|
|
3
頂いた回答をもとに完成したグラフ等を加えた
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -66,9 +66,9 @@
|
|
|
66
66
|
|
|
67
67
|
|
|
68
68
|
|
|
69
|
-
xli=[]
|
|
69
|
+
#xli=[] 座標x用のリスト
|
|
70
|
-
|
|
70
|
+
|
|
71
|
-
free=[]
|
|
71
|
+
#free=[] 自由エネルギーF用のリスト
|
|
72
72
|
|
|
73
73
|
|
|
74
74
|
|
|
@@ -92,6 +92,14 @@
|
|
|
92
92
|
|
|
93
93
|
|
|
94
94
|
|
|
95
|
+
#関数の中でxliとfreeを宣言、に変更
|
|
96
|
+
|
|
97
|
+
xli=[]
|
|
98
|
+
|
|
99
|
+
free=[]
|
|
100
|
+
|
|
101
|
+
|
|
102
|
+
|
|
95
103
|
#x=0だとlogが使えないからx1でFを求める
|
|
96
104
|
|
|
97
105
|
F= -(Ua*(1 - x1)) -(Ub*x1) -(N*z*k*x1*(1-x1)) - N*k*t*(x1*math.log(x1)+(1-x1)*math.log(1-x1))
|
|
@@ -130,7 +138,7 @@
|
|
|
130
138
|
|
|
131
139
|
|
|
132
140
|
|
|
133
|
-
|
|
141
|
+
return xli, free
|
|
134
142
|
|
|
135
143
|
|
|
136
144
|
|
|
@@ -138,9 +146,13 @@
|
|
|
138
146
|
|
|
139
147
|
for l in range(400,1200,100): #400~1100でプロットしたい
|
|
140
148
|
|
|
141
|
-
free_en(l)
|
|
149
|
+
xli, free = free_en(l)
|
|
142
|
-
|
|
150
|
+
|
|
143
|
-
plt.plot(xli,free)
|
|
151
|
+
#plt.plot(xli,free)
|
|
152
|
+
|
|
153
|
+
plt.plot(xli,free,label="t={}".format(l))#凡例を加えた
|
|
154
|
+
|
|
155
|
+
plt.legend()
|
|
144
156
|
|
|
145
157
|
|
|
146
158
|
|
|
@@ -160,7 +172,7 @@
|
|
|
160
172
|
|
|
161
173
|
|
|
162
174
|
|
|
163
|
-
free_en(400)
|
|
175
|
+
~~free_en(400)
|
|
164
176
|
|
|
165
177
|
plt.plot(xli,free, label="t=400")
|
|
166
178
|
|
|
@@ -168,7 +180,23 @@
|
|
|
168
180
|
|
|
169
181
|
|
|
170
182
|
|
|
171
|
-
↑これを500,600,700,,1100についてやったら目標のグラフができた。
|
|
183
|
+
↑これを500,600,700,,1100についてやったら目標のグラフができた。~~
|
|
184
|
+
|
|
185
|
+
|
|
186
|
+
|
|
187
|
+
|
|
188
|
+
|
|
189
|
+
for l in range(400,1200,100): #400~1100でプロットしたい
|
|
190
|
+
|
|
191
|
+
xli, free = free_en(l)
|
|
192
|
+
|
|
193
|
+
plt.plot(xli, free,label="t={}".format(l))
|
|
194
|
+
|
|
195
|
+
plt.legend()
|
|
196
|
+
|
|
197
|
+
|
|
198
|
+
|
|
199
|
+
頂いた回答に
|
|
172
200
|
|
|
173
201
|
|
|
174
202
|
|
|
@@ -193,3 +221,11 @@
|
|
|
193
221
|
温度毎に色を変えたいし凡例も入れたいけれどまず直線を消したい。
|
|
194
222
|
|
|
195
223
|
|
|
224
|
+
|
|
225
|
+
|
|
226
|
+
|
|
227
|
+
|
|
228
|
+
|
|
229
|
+
いただいた回答をもとに作成したグラフ。
|
|
230
|
+
|
|
231
|
+
|
2
不要なリストがあった
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -100,10 +100,6 @@
|
|
|
100
100
|
|
|
101
101
|
xli.append(x)
|
|
102
102
|
|
|
103
|
-
Uli.append(-(Ua*(1 - x1)) -(Ub*x1) -(N*z*k*x1*(1-x1)))
|
|
104
|
-
|
|
105
|
-
Sli.append(- N*k*(x1*math.log(x1)+(1-x1)*math.log(1-x1)))
|
|
106
|
-
|
|
107
103
|
x+=0.01
|
|
108
104
|
|
|
109
105
|
|
1
手打ちの内容を加えた。
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -160,7 +160,19 @@
|
|
|
160
160
|
|
|
161
161
|
|
|
162
162
|
|
|
163
|
-
なにが原因か分からないため試したことがないです。
|
|
163
|
+
~~なにが原因か分からないため試したことがないです。~~
|
|
164
|
+
|
|
165
|
+
|
|
166
|
+
|
|
167
|
+
free_en(400)
|
|
168
|
+
|
|
169
|
+
plt.plot(xli,free, label="t=400")
|
|
170
|
+
|
|
171
|
+
plt.legend()
|
|
172
|
+
|
|
173
|
+
|
|
174
|
+
|
|
175
|
+
↑これを500,600,700,,1100についてやったら目標のグラフができた。
|
|
164
176
|
|
|
165
177
|
|
|
166
178
|
|