# -*- coding: utf-8 -*-
"""
Created on Fri Dec 24 14:32:09 2021

@author: Sim
"""
import scipy.stats as st
import numpy as np

a = 10 
b = 4
x = st.norm.rvs(size=100)
y = a*x + b
mx = np.mean(x)
vx = np.var(x, ddof=1)
sx = pow(vx, 0.5)
print(mx, vx, sx)

my1 = np.mean(y)
vy1 = np.var(y, ddof=1)
sy1 = pow(vy1, 0.5)
print(my1, vy1, sy1)

my2 = a*mx + b
vy2 = a**2*vx
sy2 = abs(a)*sx

print(my1-my2, vy1-vy2, sy1-sy2)