For the Legendre symbol (m/p), p must be an odd prime. The Jacobi symbol (m/n) allows n to be any odd number.
Petter (2021). Jacobi and Legendre symbol (https://www.mathworks.com/matlabcentral/fileexchange/24672-jacobi-and-legendre-symbol), MATLAB Central File Exchange. Retrieved .
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!Create scripts with code, output, and formatted text in a single executable document.
There is a serious bug in your implementation of the Jacobi symbol, which was already pointed out by Sean. Please fix
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.