代码1

作者: xuqiuhao | 来源:发表于2017-01-02 13:44 被阅读0次
from __future__ import division
import numpy as np
import matplotlib.pyplot as plt

steps = np.linspace(0, 100, 101)
x_ave = np.zeros(101)
x_y0 = np.zeros(101)
x_now = np.zeros(500)

for i in range(100):
    for j in range(500):
        ruler = np.random.rand()
        if ruler<=0.5:
            x_now[j] = x_now[j] + 1
        else:
            x_now[j] = x_now[j] - 1
    average = sum(x_now)/500
    x_ave[i+1] = average
    
plt.scatter(steps, x_ave)
plt.plot(steps, x_y0)
plt.xlim(0,100)
plt.ylim(-1,1)
plt.grid(True)
plt.xlabel('step number(= time)')
plt.ylabel('<x>')
plt.title('<x> of 500 walkers')

plt.show()

相关文章

  • 代码1

  • 代码1

  • 代码代码day1

    Mar.23 题目 Write a program that concatenates two linked li...

  • 优雅代码(1)

    异常 自定义异常,设定状态码和相应的异常信息,在关键代码处打日志,并且抛出异常信息。什么是关键代码处,比如调用接口...

  • jQuery代码 1

  • Vue代码 1

    效果: 代码:

  • 代码记录1

    1 .控制器的基类 - (void)viewDidLoad { [super viewDidLoad];...

  • 代码神教1

    //代码段001H int main(void) { int God = 1; /...

  • 代码规范-1

    JS/TS部分: 函数与函数之间空一行 js/ts 函数的大括号与函数名再同一行 注意行末的分号4 注意空格 函数...

  • 代码规范(1)

    来源:https://www.cnblogs.com/Zachary-Fan/p/CodeStandards.ht...

网友评论

      本文标题:代码1

      本文链接:https://www.haomeiwen.com/subject/qzywvttx.html