function [a] = calcType1Long (n) s=0; i= sym('i'); for i=1:n if mod(n,i)==0 s=s+feval(symengine,'numlib::moebius', n/i)*(feval(symengine,'numlib::fibonacci', i+1)+feval(symengine,'numlib::fibonacci', i-1)); end end a=s/n; end