- /
-
Symmetric pattern 1
on 19 Oct 2021
- 3
- 24
- 2
- 0
- 228
k=0;
for m=1:50
for n=1:50
k=k+1;
sp=sqrt((m+1)*m-n*(n+1))*KD(m,n+1);
sm=sqrt((m+1)*m+n*(n-1))*KD(m,n-1);
sz(k)=m*KD(m,n);
sy(k)=real(sp-sm);
end
end
plot(sy,sz,-sy,sz)
hold on
plot(sy,-sz,-sy,-sz)
axis off
function v=KD(m,n)
if m==n
v=0;
else
v=1;
end
end