Exercício 4.27 (c)

> X:=matrix(5,3,[1,1,1,1,1/2,1/4,1,1/3,1/9,1,1/4,1/15,1,1/5,1/25]);
b:=vector(5,[1.2,0.8,2.4,2.9,3.5]);

> XtX:=evalm(transpose(X)&*X);
Xtb:=evalm(transpose(X)&*b);
sol:=vector(2,0):
sol:=Gauss(XtX,Xtb);

> p:=sol[1]+sol[2]/x+sol[3]/x^2;
with(student):
px:=makeproc(p,x);
with(plots):
pplot:=plot(px(x),x=1..5):
A:=pointplot({[1,b[1]],[2,b[2]],[3,b[3]],[4,b[4]],[5,b[5]]},axes=BOXED,symbol=circle,title="Exercício 4.27 (b)"):
display(pplot,A);

Warning, new definition for D

>