You cannot select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
This repo is archived. You can view files and clone it, but cannot push or open issues/pull-requests.

47 lines
981 B
Python

import matplotlib.pyplot as plt
import numpy as np
#x = [146.66, 293.33, 439.99, 573.33, 693.33, 773.33, 839.99, 879.99]
x = [0.3, 0.4, 0.5, 0.6, 0.8, 1.0]
#y = [-37, -79, -122, -161, -194, -218, -234, -245]
#y=[-53, -110, -165, -214, -261, -292, -314, -327]
#y = [-67, -137, -209, -274, -328, -368, -395, -411]
#y = [-83, -166, -249, -327, -393, -441, -474, -496]
#y = [-113, -229, -336, -442, -531, -595, -638, -664]
#y = [-142, -286, -424, -553, -663, -745, -800, -831]
#y=[135, 277, 417, 547, 660, 736, 791, 821]
y = [-0.28, -0.38, -0.47, -0.56, -0.75, -0.94]
u0 = 0
if len(x) != len(y):
print("size mismatch")
exit()
y = [i - u0 for i in y]
n = len(x)
mean_xy = 0
mean_x_sq = 0
for i in range(n):
mean_xy += x[i] * y[i]
mean_x_sq += x[i] ** 2
mean_xy /= n
mean_x_sq /= n
k = mean_xy / mean_x_sq
print("k = ", k)
plt.scatter(x,y)
#plt.show()
x_axis = np.array(range(0, 1000))
y_axis = k * x_axis
plt.plot(x_axis, y_axis)
plt.show()
#`plot(x, y))