function kval = K2(u,v)

k1 = K1(u,v);
kval = k1 + k1.^2;

