in statistischer Physik soll ich für unterschiedliche N diese Funktion über die Zeit plotten.
Code: Alles auswählen
def M(t,N):
m=0
v = [[Z0] for y in range(N)]
for i in range(t):
for j in range(N):
R=random.uniform(0,1)
if R<p:
v[j]=Z1
if R>p:
v[j]=Z0
for h in range(N):
if v[h]==Z1:
m=m+1
if v[h]==Z0:
m=m-1
return(m/N);
Code: Alles auswählen
t1 = np.linspace(0, 60, 61, endpoint=True)
F=M(t1,100)
plt.plot(t1,F)
startx, endx = -10.1, 10.1
starty, endy = -10.1, 10.1
plt.axis([startx, endx, starty, endy])
plt.show()
danke für die Hilfe