function [D,Ke,Fe] = LocalSMRT(resolution)
E = 1; m = 0.33;
p = E / ( (1+m) * (1-2*m) );
D = [p*(1-m) p*m     p*m     0         0         0
    p*m      p*(1-m) p*m     0         0         0
    p*m      p*m     p*(1-m) 0         0         0
    0        0       0       p*(0.5-m) 0         0
    0        0       0       0         p*(0.5-m) 0
    0        0       0       0         0         p*(0.5-m)];
[Ke,Fe] = GaussintKF(1.0/resolution,D);
end