f = @(n)sum(gcd(n,1:n)<2)

f(1)      % 1
f(2)      % 1
f(3)      % 2
f(8)      % 4
f(9)      % 6
f(26)     % 12
f(44)     % 20
f(105)    % 48