clear all; 
close all;
c=1;
a=1
h=a/10;
[xx yy]=meshgrid(-a:h:a,-a:h:a);
xxx=xx(:);
yyy=yy(:);
em=kermatgauss(xxx,yyy,c);
lem=length(em(:,1))
figure(1)
u=em(1:lem/2,:);
quiver(xxx,yyy,u(:,1),u(:,2))
figure(2)
v=em(lem/2+1:end,:);
quiver(xxx,yyy,v(:,1),v(:,2))
figure(3)
quiver(xxx,yyy,u(:,1)-v(:,1),u(:,2)-v(:,2))
axis([-a a -a a])