Euler's totient function

[Fotios Kasolis]

Does Octave have Euler’s totient function? If not, what do you think about including it? It could be as simple as sth like

 

function x = phi (n)

x = sum (gcd ([1:n], repmat (n, 1, n)) == 1)

endfunction

 

or slightly different to allow vector argin.

 

http://en.wikipedia.org/wiki/Euler's_totient_function

 

/Fotis