a = [1.92306958582241e+15, 3.20449986572221e+24, 1.34271290136344e+32, ...
     2.32739765751038e+38]
b = [7.33727670161595e+27, 1.05919311870816e+36, 4.56169848520627e+42]
[div, rem] = deconv (a, b)
[p, q] = filter (a, b, [1, 0])
disp('notice that q * b(1) is similar to rem')
disp('now I add an extra digit to a(2)')
a = [1.92306958582241e+15, 3.204499865722215e+24, 1.34271290136344e+32, ...
     2.32739765751038e+38];
[div, rem] = deconv (a, b)
[p, q] = filter (a, b, [1, 0])
disp('rem has a new non-zero element, while q seems fine')
disp('in Matlab, no difference between the two rem')
