Matlab - Ayuda con ejercicio

clear all

fprintf('Ingresar una Matriz Año \n \n' )

f=input('Ingrese el numero de Filas para la Matriz \n');

c=input('Ingrese el numero de Columnas para la Matriz \n');

for i=1:f

    for j=1:c

 disp(['El elemento Año(',num2str(i),',',num2str(j),')'])

 A(i,j)=input('');

%% Interfaz

fprintf('Ingresar una Matriz Año\n\n')

f=input('Ingrese el número de Filas para la Matriz: ');

c=input('Ingrese el número de Columnas para la Matriz: ');

fprintf('\n');

%% Matriz

ini = 2021;

step = 1;

for i = 1:f

    for j = 1:c

        M(i,j) = ini - step;

        step = step + 1;

    end

end

%% Programa

idx_bis = es_bisiesto(M);

Mbis = M(idx_bis);

Mbis_orden = sort(Mbis,'descend');

%% Resultados

fprintf('Matriz lógica de años bisiestos = \n\n');

disp(idx_bis);

fprintf('Años bisisestos = \n\n');

disp(Mbis_orden);

function idx_bis = es_bisiesto(M)

% Un año es bisiesto si es divisible entre 4 a menos que sea divisible

% entre 100, pero si es divisible por 100 y además por 400 también resulta bisiesto

[row, col] = size(M);

for i = 1:row

    for j = 1:col

        if (mod(M(i,j),4)==0)

            idx_bis(i,j) = logical(1);

        end

        if (mod(M(i,j),100)==0)

            if (mod(M(i,j),400)==0)

                idx_bis(i,j) = logical(1);

            else

                idx_bis(i,j) = logical(0);

            end

        end

    end

end