# restart; r1:=-Sum(x^k/k,k=1..m)-Int(t^m/(1-t),t=0..x); r2:=subs(m=10,r1); r3:=value(r2); r4:=subs(x=-3,r3); value(r4); I1:=Int(t^m/(1-t),t=0..x); I2:=subs(x=-3,I1); for n from 2 to 10 do print(n,evalf(subs(m=n,I2))); od; I3:=subs(x=3/4,I1); for n from 2 to 10 do print(n,evalf(subs(m=n,I3))); od; S1:=Sum(x^k/k,k=1..m); S2:=subs(x=3/4,S1); for n from 2 to 10 do print(n,evalf(subs(m=n,S2))); od; evalf(log(4)); evalf(subs(m=100,S2)); S3:=Sum((-1)^l*x^(2*l+1)/(2*l+1),l=0..m); S4:=subs(m=5,S3); value(S4); I4:=(-1)^(m+1)*Int(t^(2*m+2)/(1+t^2),t=0..x); I5:=subs(m=5,I4); r5:=S4+I5; r6:=subs(x=4,r5); value(r6); S4:=subs(x=4/(1+sqrt(17)),S3); for n from 2 to 10 do print(n,evalf(subs(m=n,2*S4))); od; evalf(arctan(4)); evalf(subs(m=100,2*S4)); #