Matlab: символьное вычисление членов посл-ти A006206

[maximkarimov]

kotenok gav в сообщении #1458165 писал(а):

maximkarimov в сообщении #1458162 писал(а):

s=3.980154972736918e+18

Это при $n=89$?

Да.

[kotenok gav]

maximkarimov в сообщении #1458166 писал(а):

Да.

И параллельно вы говорите, что $a006206(89)=1$.

[maximkarimov]

kotenok gav в сообщении #1458165 писал(а):

maximkarimov в сообщении #1458162 писал(а):

s=3.980154972736918e+18

Это при $n=89$?

Да.
После деления на n получаю 4.472084239030245e+16 - это неверное значение. Вот собственно и проблема.

-- 27.04.2020, 12:22 --

kotenok gav в сообщении #1458167 писал(а):

maximkarimov в сообщении #1458166 писал(а):

Да.

И параллельно вы говорите, что $a006206(89)=1$.

Сорри! Это значение, которое для любого простого числа дает A000358 минус A006206.
Это позволяет проверять верность вычисления численными методами.
Для вычисления членов A000358 получилось задействовать Symbol toolbox - там целое число на выходе сколь угодно большое без проблем получаю. А для вычисления членов A006206 использовать символьные вычисления у меня не получается. В этом собственно и прошу помощи.

[kotenok gav]

maximkarimov в сообщении #1458168 писал(а):

После деления на n получаю 4.472084239030245e+16 - это неверное значение. Вот собственно и проблема.

Что вы знаете о такой вещи, как потеря точности при делении больших чисел на малые?

-- 27 апр 2020, 19:00 --

maximkarimov в сообщении #1458168 писал(а):

4.472084239030245e+16

Вы, наверное, удивитесь, но это правильное значение.

[maximkarimov]

kotenok gav в сообщении #1458170 писал(а):

maximkarimov в сообщении #1458168 писал(а):

После деления на n получаю 4.472084239030245e+16 - это неверное значение. Вот собственно и проблема.

Что вы знаете о такой вещи, как потеря точности при делении больших чисел на малые?

Нет, не слышал. Но вижу. Хочу решить проблему вычисления членов А006206 задействуя символьные вычисления, но не получается.
Можете предложить другой путь?
Вот работающий код для точного (символьного) вычисления членов А000358 (на всякий случай):

Используется синтаксис Matlab M

function [a] = calcType1Total (n)
    % подсчитывает А000358 по формуле Ходжикостаса
    n=sym(n);
    s=0;
    for i=1:n
        if mod(n,i)==0
            s=s+totient(n/i)*(fibonacci(i+1)+fibonacci(i-1));
        end
    end
    a=s/n;
end
 

[kotenok gav]

maximkarimov в сообщении #1458171 писал(а):

Хочу решить проблемы задействуя символьные вычисления, но не получается.

Давайте по-другому. Зачем вам вообще вычислять эту последовательность?

[maximkarimov]

kotenok gav в сообщении #1458170 писал(а):

maximkarimov в сообщении #1458168 писал(а):

После деления на n получаю 4.472084239030245e+16 - это неверное значение. Вот собственно и проблема.

Что вы знаете о такой вещи, как потеря точности при делении больших чисел на малые?

-- 27 апр 2020, 19:00 --

maximkarimov в сообщении #1458168 писал(а):

4.472084239030245e+16

Вы, наверное, удивитесь, но это правильное значение.

Угу, но не точное.
Помогите получить точное.

-- 27.04.2020, 12:35 --

kotenok gav в сообщении #1458172 писал(а):

maximkarimov в сообщении #1458171 писал(а):

Хочу решить проблемы задействуя символьные вычисления, но не получается.

Давайте по-другому. Зачем вам вообще вычислять эту последовательность?

Любопытно. Вас такой ответ устроит?

[kotenok gav]

maximkarimov в сообщении #1458174 писал(а):

Угу, но не точное.

Вы, наверное, удивитесь, но это точное значение.

[maximkarimov]

kotenok gav в сообщении #1458175 писал(а):

maximkarimov в сообщении #1458174 писал(а):

Угу, но не точное.

Вы, наверное, удивитесь, но это точное значение.

То есть ошибка возникает при делении этого числа на 89?
А как ее избежать? Умоляю, скажите!!! Очень нужно!

[kotenok gav]

maximkarimov в сообщении #1458176 писал(а):

То есть ошибка возникает при делении этого числа на 89?

Вы, наверное, удивитесь, но она не возникает.

maximkarimov в сообщении #1458174 писал(а):

Любопытно. Вас такой ответ устроит?

Если вы хотите найти какую-то закономерность - есть http://oeis.org/A006206/b006206.txt.

[maximkarimov]

kotenok gav в сообщении #1458177 писал(а):

maximkarimov в сообщении #1458176 писал(а):

То есть ошибка возникает при делении этого числа на 89?

Вы, наверное, удивитесь, но она не возникает.

maximkarimov в сообщении #1458174 писал(а):

Любопытно. Вас такой ответ устроит?

Если вы хотите найти какую-то закономерность - есть http://oeis.org/A006206/b006206.txt.

Вот!!! Смотрим 89 член по этой ссылке: 44720842390302450
Далее смотрим А000358 (89): 44720842390302451
Разница равна 1 (как и для любого n=p).
Второе число я вычисляю точно, а первое - не получаецца!
Как вычислить первое число в Матлабе - подскажите, прошу Вас!

[kotenok gav]

maximkarimov в сообщении #1458181 писал(а):

44720842390302450

maximkarimov в сообщении #1458168 писал(а):

4.472084239030245e+16

Вы видите разницу между этими значениями?

[GAA]

maximkarimov, можно проще в MuPAD
Matlab R2013b::MuPAD

A006206 := proc(n)
  local sum, d;
 begin
     sum := 0;
     for d in numlib::divisors(n) do 
       sum := sum + numlib::moebius(n/d)*(numlib::fibonacci(d+1)+numlib::fibonacci(d-1))
     end;
    sum/n;
 end;

A006206(89)
    44720842390302450 

Вложение:

Комментарий к файлу: Matlab R2013b::MuPAD

A006206.PNG [ 6.18 Кб | Просмотров: 2872 ]

[maximkarimov]

kotenok gav в сообщении #1458182 писал(а):

maximkarimov в сообщении #1458181 писал(а):

44720842390302450

maximkarimov в сообщении #1458168 писал(а):

4.472084239030245e+16

Вы видите разницу между этими значениями?

У числа 4.472084239030245e+16 не видно значение в разряде, который следует за 5 в конце.
Должна быть 1, но когда я вычитаю из него 44720842390302450, то получаю то 3, то 2.

[kotenok gav]

maximkarimov в сообщении #1458184 писал(а):

Должна быть 1, но когда я вычитаю из него 44720842390302450, то получаю то 3, то 2.

Это называется "потерей точности".