Ima3PC2

n=[1:10];
an=1./n.^2;
sum(an)

n=[1:100];
an=1./n.^2;
sum(an)

n=[1:200];
an=1./n.^2;
sum(an)

n=[1:300];
an=1./n.^2;
sum(an)

syms k
S=symsum(1/k^2,k,1,inf)
double(S)

n=[1:10];
an=1./n;
sum(an)

n=[1:100];
an=1./n;
sum(an)

n=[1:200];
an=1./n;
sum(an)

n=[1:300];
an=1./n;
sum(an)

syms x
c = 1;
n = 2;
f=x/(x^2+1)
Pn=taylor(f,n+1,c)

syms x
c = 1;
f=x/(x^2+1)

X=[0.5:0.01:1.5]
Y=subs(f,x,X)
plot(X,Y,'g')

syms x
c = 1;
n = 3;
f=x/(x^2+1)
Rn=taylor(f,n+1,c)

syms x
c = 1;
n = 4;
f=x/(x^2+1)
Qn=taylor(f,n+1,c)

X=[0.5:0.01:1.5]
Y=subs(Pn,x,X)
plot(X,Y,'k')

hold on

X=[0.5:0.01:1.5]
Y=subs(Rn,x,X)
plot(X,Y,'b')

X=[0.5:0.01:1.5]
Y=subs(Qn,x,X)
plot(X,Y,'r')