import matplotlib.pyplot as plt
import numpy as np
 
# 数据
d_squared = [16, 2.25, 6.25]  # 光阑直径的平方 (mm²)
I_saturated = [51, 1279, 18.1]  # 饱和光电流 (pA)
 
# 绘制散点图
plt.scatter(d_squared, I_saturated, color='blue', label='Data Points')
 
# 添加线性拟合线（可选）
coefficients = np.polyfit(d_squared, I_saturated, 1)
polynomial = np.poly1d(coefficients)
x_fit = np.linspace(min(d_squared), max(d_squared), 100)
y_fit = polynomial(x_fit)
plt.plot(x_fit, y_fit, color='red', label='Linear Fit')
 
# 添加标题和标签
plt.title('Saturated Photocurrent vs. Aperture Diameter Squared')
plt.xlabel('Aperture Diameter Squared (mm²)')
plt.ylabel('Saturated Photocurrent (pA)')
plt.legend()
 
# 显示图像
plt.grid(True)
plt.show()# your code goes here

aW1wb3J0IG1hdHBsb3RsaWIucHlwbG90IGFzIHBsdAppbXBvcnQgbnVtcHkgYXMgbnAKCiMg5pWw5o2uCmRfc3F1YXJlZCA9IFsxNiwgMi4yNSwgNi4yNV0gICMg5YWJ6ZiR55u05b6E55qE5bmz5pa5IChtbcKyKQpJX3NhdHVyYXRlZCA9IFs1MSwgMTI3OSwgMTguMV0gICMg6aWx5ZKM5YWJ55S15rWBIChwQSkKCiMg57uY5Yi25pWj54K55Zu+CnBsdC5zY2F0dGVyKGRfc3F1YXJlZCwgSV9zYXR1cmF0ZWQsIGNvbG9yPSdibHVlJywgbGFiZWw9J0RhdGEgUG9pbnRzJykKCiMg5re75Yqg57q/5oCn5ouf5ZCI57q/77yI5Y+v6YCJ77yJCmNvZWZmaWNpZW50cyA9IG5wLnBvbHlmaXQoZF9zcXVhcmVkLCBJX3NhdHVyYXRlZCwgMSkKcG9seW5vbWlhbCA9IG5wLnBvbHkxZChjb2VmZmljaWVudHMpCnhfZml0ID0gbnAubGluc3BhY2UobWluKGRfc3F1YXJlZCksIG1heChkX3NxdWFyZWQpLCAxMDApCnlfZml0ID0gcG9seW5vbWlhbCh4X2ZpdCkKcGx0LnBsb3QoeF9maXQsIHlfZml0LCBjb2xvcj0ncmVkJywgbGFiZWw9J0xpbmVhciBGaXQnKQoKIyDmt7vliqDmoIfpopjlkozmoIfnrb4KcGx0LnRpdGxlKCdTYXR1cmF0ZWQgUGhvdG9jdXJyZW50IHZzLiBBcGVydHVyZSBEaWFtZXRlciBTcXVhcmVkJykKcGx0LnhsYWJlbCgnQXBlcnR1cmUgRGlhbWV0ZXIgU3F1YXJlZCAobW3CsiknKQpwbHQueWxhYmVsKCdTYXR1cmF0ZWQgUGhvdG9jdXJyZW50IChwQSknKQpwbHQubGVnZW5kKCkKCiMg5pi+56S65Zu+5YOPCnBsdC5ncmlkKFRydWUpCnBsdC5zaG93KCkjIHlvdXIgY29kZSBnb2VzIGhlcmU=