function modelmat=creatcatcardeucl45 th=[0 0 pi/2 pi/2 pi/2 pi/2 pi pi pi pi pi*1.5 pi*1.5 pi*1.5 pi*1.5 0 0]+pi/4; [x,y]=pol2cart(th',1); % draw unit circle in Cartesian coords xy=[x y]; % cat for easier indexing below for a=1:length(th) for b=1:length(th) modelmat(a,b)=norm(xy(a,:) - xy(b,:)); %euclidean end end end