function [kernel_x,kernel_y]=gc_kernel(I,T,E)

  kernel_x = (0:700);
  
  kernel_y = I * exp( 1 + (E./(8.617e-5.*kernel_x)).*((kernel_x-T)./T) - (kernel_x.*kernel_x)./(T*T) .* exp((E./(8.617e-5.*kernel_x)).*((kernel_x-T)./T).*(1-(E.*8.617e-5.*kernel_x./E)))-(E*8.617e-5*T/E));
  
  