function [ Area,Omkrets ] = huron( x,y )
xx=[x;x(1)]; yy=[y;y(1)];
for i=1:3
 s(i)=sqrt((xx(i+1)-xx(i))^2+(yy(i+1)-yy(i))^2);
end
Omkrets=0;
for i=1:3
Omkrets=Omkrets+s(i);
end
p=Omkrets/2;
Area=sqrt(p*(p-s(1))*(p-s(2))*(p-s(3)));
end

