Exercício 1.7

> restart:

> Digits:=3;

> S:=0.0;
for k from 1 to 3 do
S:=S+evalf(1/(3^k));
od;

> S_exacto:=0.499314129;
Delta:=abs(S_exacto-S); # Cálculo do erro absoluto

> delta:=Delta/abs(S_exacto); # Cálculo do erro relativo

> S_maple:=evalf(sum(1/(3^i),i=1..3),10);
Delta:=abs(evalf(S_maple-S,10));
delta:=Delta/abs(S_maple);