суббота, 19 февраля 2011 г.

CORDIC для SciLab, заодно тестируем вставку кода

Пока без подсветки, но все же.
//заполняем таблицу
i=1:10;
log10_table_p = log10(1 + 2^(-i));
log10_table_m = log10(1 - 2^(-i));

//выводим значение логарифма
x = 2;
printf("%f\n", log10(x));

//CORDIC
z = 0;
for i = 1:10
  if (x > 1)
    x = x - x*2^(-i);
    z = z - log10_table_m(i);
  else
    x = x + x*2^(-i);
    z = z - log10_table_p(i);
  end;
end;

printf("%f", z);

Комментариев нет:

Отправить комментарий