function res = c_approx_gamma( t )

d1=log(2)/log(3)+2*i*pi/log(3);
d=log(2)/log(3);
d2=log(2)/log(3)-2*i*pi/log(3);
c_approx_gamma=t^(-d)*(coefc(-1)*d1*exp(2*i*pi*log(t)/log(3))*(exp(i*pi/2*(d1-1))-exp(-i*pi/2*(d1-1)))/(2*i)*(-0.0001)+coefc(0)*d*sin(pi/2*(d-1))*(-3.8537)+coefc(1)*d2*exp(-2*i*pi*log(t)/log(3))*(exp(i*pi/2*(d2-1))-exp(-i*pi/2*(d2-1)))/(2*i)*(-0.0001));
res=c_approx_gamma;
end

