function [w]=winkel(px1,py1,px2,py2,px3,py3,px4,py4)

w=acos(((px1-px2)*(px3-px4)+(py1-py2)*(py3-py4))/...
        (sqrt((px1-px2)^2+(py1-py2)^2 )*sqrt((px3-px4)^2+(py3-py4)^2)) );
    
     if w>pi/2, w=pi-w;end
     
        