Перенос цикла из matlab в c ++

Я написал это для цикла в matlab:

b1=zeros(size(raw));
a1=zeros(size(raw));
filtered = zeros(size(raw));

raw=ones(124,1);
[b,a]=butter(4,6/64,'high');

for n = 5:124

    b1(n)=(b(1)* raw(n))   + (b(2)* raw(n-1) )+ (b(3)* raw(n-2))+ (b(4)* raw(n-3) );
    a1(n)=(b(5)* raw(n-4))-    (a(2)*filtered(n-1)) - (a(3)*filtered(n-2))-(a(4)*filtered(n-3) )- (a(5)*filtered(n-4));

    filtered(n)=b1(n)+a1(n);
end

и это в c++

for(int i=0;i<124;i++)
{
    temp[i]=1;
}
long double f_temp[124]={0};

long double a[5]={1.0000,-3.2317,3.9766,-2.2014,0.4617};    
long double b[5]={ 0.6795,-2.7178,4.0768,-2.7178,0.6795};

for(int n=4;n<124;n++)
{
    b1=(b[0]*temp[n])+(b[1]*temp[n-1])+(b[2]*temp[n-2]) +(b[3]*temp[n-3]);
    a1=(b[4]*temp[n-4])-(a[1]*f_temp[n-1])-(a[2]*f_temp[n-2])-(a[3]*f_temp[n-3])-(a[4]*f_temp[n-4]);

    f_temp[n]=b1+a1;
}

Может ли кто-нибудь объяснить, почему они не дают того же ответа? Входные данные те же, и выходной массив инициализируется как ноль в обоих случаях.