
Sxy_table = {{0}}
Sxy( Sxy_table ) =
	sum ( Sxy_table[ $1*$2 ] ) - ( 1 / count ( Sxy_table ) )
	* sum ( Sxy_table[ $1 ] ) * sum ( Sxy_table [ $2 ] )

beta_table = {{0}}
beta_func ( beta_table ) = 
	Sxy ( beta_table ) / Sxy ( beta_table[ $1 , $1 ] )

alpha_table = {{0}}
alpha_func ( alpha_table ) =
	average ( alpha_table[ $2 ] )
	- beta_func ( alpha_table ) * average ( alpha_table[ $1 ] )

fit_Min = 0
fit_Max = 0
fit_Alpha = 0
fit_Beta = 0
line_fit2 ( fit_Min , fit_Max , fit_Alpha , fit_Beta ) =
	{{ fit_Min , fit_Alpha + fit_Beta * fit_Min }
	{  fit_Max , fit_Alpha + fit_Beta * fit_Max }}

line_fit_table = {{0}}
line_fit ( line_fit_table ) =
	line_fit2 ( min ( line_fit_table[ $1 ] ) , max ( line_fit_table[ $1 ] ),
	alpha_func ( line_fit_table ) , beta_func ( line_fit_table ) )

