function [u,v]= fake(n) %don't need this function-doesn't seem to work anyway w=rand(4,n); w=sort(w); h=w(4,:); l=w(1,:); v=floor(2*rand(1,n)); u=w(4,:)-w(1,:); o=w(2+v,:); c=w(3-v,:); u=(h-o)./(2*h-o-c); v=(l-o)./(2*l-o-c);