function result = K_tau_22(x,y,t)
temp1 = AiryKernel(x,y);

[w, lambda]=ClenshawCurtis(0,100,256);
temp_integral = w*(airy(lambda+y).*exp(t.*lambda));

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

%result = real(result*conjugation(x,y,t));
end


function s = conjugation(x,y,t)
    s=exp(t*(x-y));
end