`y=ax^b':

let(data,<
<0.50,(1+4)        /2>,
<0.75,(4.1+9)      /2>,
<1.00,(9.1+16.5)   /2>,
<1.25,(16.6+28)    /2>,
<1.50,(28.1+55)    /2>,
<2.00,(55.1+107.5) /2>,
<2.50,(107.6+182.5)/2>,
<3.00,(182.6+287.5)/2>,
<3.50,(287.6+425)  /2>,
<4.00,(425.1+700)  /2>,
<5.00,(700.1+1100) /2>,
<6.00,(1100.1+1300)/2>
>):

func(xayb,{for(<j,1,len([data]),1>,{
  ((ln(elem(elem([data],[j]),1)))^[x]) * 
  ((ln(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],eul([b])>