import numpy as np, konstanter as k

def usikkerhet_x(T_points):
    T_points = np.array(T_points)
    return (1/(k.R*(T_points**2)))*k.s_T

def usikkerhet_y(p_points):
    p_points = np.array(p_points)
    return (1/p_points)*k.s_p

def varians(y_points=list, y_line=list):
    y_points = np.array(y_points)
    y_line = np.array(y_line)
    dev = y_points - y_line
    mean_dev = sum(abs(dev))/len(dev)

    var = sum((dev-mean_dev)**2)/(len(dev)-1)
    return np.sqrt(var), mean_dev