%  short  M-file (solution by a student)

for i=10^6+1:8:10^7
	for j=primes(50)
	    	if 1 ~= jac(j,i)
	  		break
	end
	if j == 47
			y = sqrt(i);
			if floor(y) ~= y
				disp ('!!!!!!!!!!')
				disp( i)
			end
		end
	end
end

%  not an efficient algorithm  (for Matlab)
%  but is doing the job  (in an hour or so)