114北模數學A 第16題

1
import numpy as np
2
import matplotlib.pyplot as plt
3
from matplotlib.animation import FuncAnimation
4
circle_radius = 5
5
slope = -3 / 4
6
a_values = np.arange(1, 8)
7
b_values = np.arange(1, 8)
8
fig, ax = plt.subplots(figsize=(8, 10))
9
ax.set_xlim(-10, 10)
10
ax.set_ylim(-15, 10)
11
ax.set_aspect('equal', adjustable='box')
12
theta = np.linspace(0, 2 * np.pi, 100)
13
x_circle = circle_radius * np.cos(theta)
14
y_circle = circle_radius * np.sin(theta)
15
ax.plot(x_circle, y_circle, color="blue", label="Circle C: $x^2 + y^2 = 25$")
16
ax.grid(True, linestyle='--', alpha=0.7)
17
line, = ax.plot([], [], 'r', lw=1.5, alpha=0.7)
18
dice_a_text = ax.text(-7, 8, '', fontsize=15, ha='center')
19
dice_b_text = ax.text(7, 8, '', fontsize=15, ha='center')
20

21
def init():
22
    line.set_data([], [])
23
    dice_a_text.set_text('')
24
    dice_b_text.set_text('')
25
    return line, dice_a_text, dice_b_text
26

27
def uppubDate(frame):
28
    a, b = frame
29
    intercept = (3 * a + 4 * b) / 4
30
    x_vals = np.linspace(-10, 10, 100)
31
    y_vals = slope * x_vals + intercept
32
    line.set_data(x_vals, y_vals)
33
    dice_a_text.set_text(f'dice A: {a}')
34
    dice_b_text.set_text(f'dice B: {b}')
35
    return line, dice_a_text, dice_b_text
36

37
frames = [(a, b) for a in a_values for b in b_values]
38
ani = FuncAnimation(fig, uppubDate, frames=frames, init_func=init, blit=True, interval=300)
39
plt.xlabel("x")
40
plt.ylabel("y")
41
plt.legend()
42
plt.show()