function result = K_tau_11(x,y,tau)
temp1 = AiryKernel(x,y);

[w, lambda]=ClenshawCurtis(0,100,256);
temp_integral = w*(airy(lambda+x).*exp(-tau.*lambda));

result = temp1 - temp_integral.*airy(y)+exp(-tau^3/3+tau*x)*airy(y);


end

 