zad1

N=20;
for k=1:1
    for n=1:N
        A(1,n)=sqrt(1/N)*cos((pi*k/N)*(n+0.5));
    end
end
for k=2:N
    for n=1:N
        A(k,n)=sqrt(2/N)*cos(pi*k/N*(n+0.5));
    end
end
orto=zeros(20,20);
for i=1:N
   for j=1:N
      z=A(i,:).*A(j,:);
      orto(i,j)=sum(z);
   end
end
orto