Trying to solve a simple matrix multiplication operation but getting a single error.

M2N.m:

% 1D PSHT Matlab code
function M=M2N(n)
n=2;
B=4;
if n==1
  M=eye(2);
else
  p=2^(n-1)-2; m=2^(2*n-2);
  a=sqrt(3*m/(4*m-B)) ; b=sqrt((m-B)/(4*m-B));
  M=[[1 0; 0 b] zeros(2,p) [0 0; a 0] zeros(2,p);
      zeros(p,2) eye(p) zeros(p,2) zeros(p,p);
      [0 0 ; 0 a] zeros(2,p) [0 1 ;-b 0] zeros(2,p);
      zeros(p,2) zeros(p,p) zeros(p,2) eye(p)];
end
end

psht.m:

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function Y=psht()
n=2;
for i=1:n
  Y=kron(eye(2^(n-i)),M2N(i))* kron(eye(2^(n-i)),kron(hadamard(2),eye(2^(i-1))))*Y;
end
end

[Walter Roberson]

Why are you overwriting the input argument n, in M2N ?

[Aziuzl Hakim Chowdhury]

I was checking the codes individually so used n in both cases.Al though I escaped the n=2 in M2N but the error in the psht is still shown. :(

[Walter Roberson]

You do not initialize Y in psht . If I initialize it to 1, then the code runs for me.

[Aziuzl Hakim Chowdhury]

SIr I have added Y=1 in psht.m but still the error is showing.Would you please give me the code of psht portion?

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function Y=psht()
Y=1;
n=2;
for i=1:n
Y=kron(eye(2^(n-i)),M2N(i))* kron(eye(2^(n-i)),kron(hadamard(2),eye(2^(i-1))))*Y;
end
end

[Walter Roberson]

Your code works fine for me in R2019a.

What is the complete error message? And also which MATLAB release are you using?

In particular are you using R2016a or earlier?