Closed gusang123 closed 5 months ago
我试了代码,没问题啊。是不是可以最后加一句 plt.show()
On Fri, Feb 2, 2024 at 3:21 AM gusang123 @.***> wrote:
姜博士您好,我在学习您的第三本书,在跑这个代码时,发现无法输出图像,想和您请教一下:
代码如下: import numpy as np import matplotlib.pyplot as plt
n_start = 6
B_6 = np.sin(np.pi/n_start)n_start A_6 = np.tan(np.pi/n_start)n_start
B_array = [] A_array = [] n_array = [6,12,24,48,96]
B_i = B_6 A_i = A_6 n_i = n_start
for i in n_array:
B_array.append(B_i) A_array.append(A_i)
updating
A_i = 2A_iB_i/(A_i + B_i) B_i = np.sqrt(A_i*B_i)
B_array = np.array(B_array) A_array = np.array(A_array) n_array = np.array(n_array)
fig, ax = plt.subplots()
plt.axhline(y=np.pi, color='r', linestyle='-') plt.plot(n_array,B_array, color = 'b', marker = 'x') plt.plot(n_array,A_array, color = 'g', marker = 'x') plt.fill_between(n_array, B_array, A_array, color = '#DEEAF6') plt.tight_layout() plt.xticks([6,12,24,48,96]) plt.xlim((6,96)) plt.xlabel('Number of sides, n') plt.ylabel('Estimate of $\pi$')
盼回复,谢谢
— Reply to this email directly, view it on GitHub https://github.com/Visualize-ML/Book3_Elements-of-Mathematics/issues/151, or unsubscribe https://github.com/notifications/unsubscribe-auth/AZHC6V2C4HJIF5SPKHMTRBTYRSOYDAVCNFSM6AAAAABCWIRDZ6VHI2DSMVQWIX3LMV43ASLTON2WKOZSGEYTIMZXGA2TSNI . You are receiving this because you are subscribed to this thread.Message ID: @.***>
姜博士您好,我在学习您的第三本书,在跑这个代码时,发现无法输出图像,想和您请教一下:
代码如下: import numpy as np import matplotlib.pyplot as plt
n_start = 6
B_6 = np.sin(np.pi/n_start)n_start A_6 = np.tan(np.pi/n_start)n_start
B_array = [] A_array = [] n_array = [6,12,24,48,96]
B_i = B_6 A_i = A_6 n_i = n_start
for i in n_array:
B_array = np.array(B_array) A_array = np.array(A_array) n_array = np.array(n_array)
fig, ax = plt.subplots()
plt.axhline(y=np.pi, color='r', linestyle='-') plt.plot(n_array,B_array, color = 'b', marker = 'x') plt.plot(n_array,A_array, color = 'g', marker = 'x') plt.fill_between(n_array, B_array, A_array, color = '#DEEAF6') plt.tight_layout() plt.xticks([6,12,24,48,96]) plt.xlim((6,96)) plt.xlabel('Number of sides, n') plt.ylabel('Estimate of $\pi$')
盼回复,谢谢