let(data,
<<0.6,1.6>,
<1.8,1.6>,
<2.8,2.6>,
<3.6,2.0>,
<4.2,4.0>,
<5.6,3.6>>):

func(xayb,{for(<j,1,len([data]),1>,{
  ((elem(elem([data],[j]),1))^[x]) * 
  ((elem(elem([data],[j]),2))^[y]) +
})0}):

let(x1,xayb(1,0)):
let(y1,xayb(0,1)):
let(xy,xayb(1,1)):
let(x2,xayb(2,0)):

let(n,len([data])):

let(m,
  ([n]*[xy]-[x1]*[y1]) /
  ([n]*[x2]-[x1]*[x1])
):

let(b,
  ([y1]*[x2]-[x1]*[xy]) /
  ([n]*[x2]-[x1]*[x1])
):

print_<[m],[b]>