%This program calculates the function value at t given the laplace transform of 
%the function and also gives an error bound.
%Laplace transform is given in laplace2.m
% input=t, define Laplace tranform in laplace2(x)
function euler1=inversion2(which,t,Ps1,Ps2,delta,mu,p2not,p2,r,s,x,q,xnot,qnot,beta,mndpevt,littlegama)
if t==0
   euler1=1.0;
else
A=20;
m=11;
n=100;
a=zeros(m+n+1,1);
S=zeros(n+m+1,1);
anot=exp(A/2)*laplace2(which,A/(2*t),Ps1,Ps2,delta,mu,p2not,p2,r,s,x,q,xnot,qnot,beta,mndpevt,littlegama)/(2*t);
anot=anot+exp(A/2)*real(laplace2(which,(A/2+i*pi)/t,Ps1,Ps2,delta,mu,p2not,p2,r,s,x,q,xnot,qnot,beta,mndpevt,littlegama)*exp(i*pi))/t;
for k=1:(m+n+1)
   a(k)=exp(A/2)*real(laplace2(which,(A/2+(1+k)*i*pi)/t,Ps1,Ps2,delta,mu,p2not, p2,r,s,x,q,xnot,qnot,beta,mndpevt,littlegama)*exp(i*pi))/t;
end
S(1)=anot-a(1);
for k=2:n+m+1
   S(k)=S(k-1)+((-1)^k)*a(k);
end
euler1=S(n);
for k=1:m
   euler1=euler1+nchoosek(m,k)*S(n+k);
end
euler1=(2^(-m))*euler1;
end
   
%euler2=S(n+1);
%for k=1:m
%   euler2=euler2+nchoosek(m,k)*S(n+k+1);
%end
%euler2=(2^(-m))*euler2;
% euler1 
% euler2-euler1