Хемминг в маткаде (перевод программы из Matlab'a в Mathcad)

function Y = KOD_H
fid = fopen('a.txt', 'r');     % открытие файла на запись 
if fid == -1                     % проверка корректности открытия 
    error('File is not opened'); 
end 
Tmp = fread (fid, 'char');
fclose(fid);                % закрытие файла 
Array = Tmp;
S = size (Array);
s = S(1);
 
%ЗАНЕСЕНИЕ ИСХОДНЫХ ДАННЫХ В ДВОИЧНУЮ МАТРИЦУ
Rassh=zeros(s,12);
for i=1:S(1)
    ra = Array(i);
    step = 0;
    for j=12:-1:1
        if (j==1 || j==2 || j==4 || j==8)
            continue
        end
        sd = bitshift (ra, -step, 1);
        Rassh(i,j) = sd;
        step = step+1;
    end
end
%ВЫЧИСЛЕНИЕ КОНТРОЛЬНЫХ БИТОВ
for i=1:S(1)
    for j=1:4
        n = 2^(j-1);
        sum1 = 0;
        for k = 1:n
            ind = n-1+k;
            while (ind<=12)
                sum1 = bitxor(sum1, Rassh(i,ind));
                ind = ind+(2^j);
            end
        end
        Rassh(i,n)=sum1;
    end
end
%ЗАПИСЬ МАТРИЦЫ В ФАЙЛ
fid = fopen('b.txt','w');
fprintf(fid,'%d',Rassh');
fclose(fid);
 
Y = Rassh;
end

function Y = DEKOD_H
s = fopen('b.txt','r+'); 
a = fread(s,'logical'); 
fclose(s);
%ЗАПИСЬ В МАССИВ ПО 12 ЭЛЕМЕНТОВ В СТРОКЕ
size1=size(a);
str1=size1(1);
b=zeros(str1/12,12);
sk1=1;
sb1=1;
for l=1:str1    
    b(sk1,sb1) = bitshift (a(l), 0, 1);
    sb1=sb1+1;
    if sb1==13;
       sk1=sk1+1;
       sb1=1;
    end
end
%ВЫЧИСЛЕНИЕ КОНТРОЛЬНЫХ БИТОВ
size1=size(b);
S=size1(1);
NaN1=zeros(S,4);
NaN2=NaN1;
%сохраненеи контрольных битов в матрицу
for i=1:4
    NaN1(:,i)=b(:,2^(i-1));
    b(:,2^(i-1))=zeros;
end
 
%ПОВТОРНОЕ ВЫЧИСЛЕНИЕ КОНТРОЛЬНЫХ БИТОВ ДЛЯ ПРОВЕРКИ
for i=1:S(1)
    for j=1:4
        n = 2^(j-1);
        sum1 = 0;
        for k = 1:n
            ind = n-1+k;
            while (ind<=12)
                sum1 = bitxor(sum1, b(i,ind));
                ind = ind+(2^j);
            end
        end
        NaN2(i,j)=sum1;
    end
end
 
%ПРОВЕРКА И ИСПРАВЛЕНИЕ ОШИБКИ
 
for i=1:S(1)
    if NaN2(i,:)==NaN1(i,:)
        continue
    end
    ind = 0;
    sm1 = 0;
    sm2 = 0;
    for j=4:-1:1
       d = 2^ind;
       ind = ind+1;
       sm1 = d*NaN1(i,j)+sm1;
       sm2 = d*NaN2(i,j)+sm2;
    end
    sm1=bitxor(sm1,sm2);
    'Найденна и исправлена ошибка'
    b(i,sm1)=not(b(i,sm1));
end
 
b (:,1)=[ ];
b (:,1)=[ ];
b (:,2)=[ ];
b (:,5)=[ ];
 
G=zeros(S(1),1);
for i=1:S(1)
    ind = 0;
    sm1 = 0;
    for j=8:-1:1
       d = 2^ind;
       ind = ind+1;
       sm1 = d*b(i,j)+sm1;
    end
    G(i) = sm1;
end
 
fidl = fopen ('c.txt', 'wb+');
fwrite (fidl, G, 'char');
fclose(fidl);
 
Y = G;
 
end