clearall;formatlonge;clc;N=50;tp1=0;tp2=0;pi1=[];pi2=[];fork=0:Ntp1=tp1+2*sqrt(2)/99/99/(factorial(k))^4*(26390*k+1103)/(396^(4*k))*factorial(4*k);tp2=tp2+12*(-1)^k/factorial(3*k)/(factorial(k))^3*(545140134*k+13591409)/(640320)^(3*k+3/2)*factorial(6*k);pi1=[pi1,1/tp1];pi2=[pi2,1/tp2];end;err1=pi-pi1;err2=pi-pi2;figure,subplot(1,2,1),plot(log10(abs(err1))),gridon;subplot(1,2,2),plot(log10(abs(err2))),gridon;