function [e1x,e1y]=behelp(a1x,a1y,l1,c1);

e1x=zeros(length(a1x),1); e1y=zeros(length(a1x),1);
for j=1:length(a1x),
	dmin=10000;
	for k=1:l1,
		d=(a1x(j)-c1(k,1))^2+(a1y(j)-c1(k,2))^2;
		if d<dmin, kmin=k; dmin=d; end;
	end
	e1x(j)=c1(kmin,1)-a1x(j); e1y(j)=c1(kmin,2)-a1y(j);
end
