Jacobi and Legendre symbol

JACOBI computes the Jacobi symbol (m/n), a generalization of the Legendre symbol.

Description

For the Legendre symbol (m/p), p must be an odd prime. The Jacobi symbol (m/n) allows n to be any odd number.

• MATLAB release: MATLAB 7.5 (R2007b)

[Sean McClintock]

There are a couple of bugs in this program...

In line 36, we want to test if n=+/-1 mod8, however, since matlab will return a value from 0 to 7 for mod(n,8),

if abs(mod(n,8))==1 is not adequate,

Perhaps

if mod(n,8)==1
    j = jacobi(m/2,n);
elseif mod(n,8)==7
    j = jacobi(m/2,n);
etc would be better.

Also, I don't think this program deals with negative m.