Привет Всем!

Я запрограммировал расчет одной периодической последовательности.
Алгоритм очень простой: задаем начальное число x(0) < 1, следующие числа последовательности считаются следующим образом: x(n)=2 * x(n-1) и если полученный x(n) > 1 то x(n)=x(n) - 1.

Однако при вычислении последовательности с начальным значением 0.7 получается следующее:

0,7 0,4 0,8 0,6 0,199999999999999 0,399999999999999 0,799999999999997 0,599999999999994 0,199999999999989 0,399999999999977 0,799999999999954 0,599999999999909 и т.д.

Весь этот процесс быстро сходиться к единице, чего быть не должно.

Скажите пожалуйста, откуда вместо 0,2 появляется 0,199999999999999 и как этом можно убрать???

Код следующий:
_________________________________________________________________

function bnf = bernoise(xs,l)
%
% Bernulli noise
% Функция bernoise(xs,l) используется для вычисления значений
% отображения Бернулли с начальным значением xs и длиной l
%
% Пример использования:
% bn=bernoise(0.7,100);


bn=zeros(1,5000);

bn(1)=xs;

for n=2:1:length(bn)

bn(n)=2*bn(n-1);

if bn(n)>1
bn(n)=bn(n)-1;
end

end

bnf=bn;

plot(bn);
hold on;