kill(all)$
numer:true$
fpprintprec:6$
ratprint:false$
f(x):=sin(x)$
a:-6$ b:6$ n:7$ h:(b-a)/n$
xu:create_list(a+i*h,i,0,7)$
fz:create_list(f(a+i*h),i,0,7)$
ft:create_list([a+i*h,f(a+i*h)],i,0,7)$
m:4$
phi(i):=xu^i$
mh[i,j]:=phi(i).phi(j)$
p:genmatrix(mh,m,m,0,0);
bb:create_list(fz.phi(i),i,0,m);
s:invert(p).bb;
kill(x);mnk(x):=sum(s[i+1,1]*x^i,i,0,m);

mnk(x);

plot2d([f,mnk],[x,-6,6]);