Mercurial > octave-nkf

comparison scripts/specfun/primes.m @ 10549:95c3e38098bf

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Untabify .m scripts
author Rik <code@nomad.inbox5.com>
date Fri, 23 Apr 2010 11:28:50 -0700
parents 8c71a86c4bf4
children c6833d31f34e
comparison
equal deleted inserted replaced
10548:479536c5bb10 10549:95c3e38098bf
81 endif 81 endif
82 endfor 82 endfor
83 x = [2, 1+2*find(sieve)]; # primes remaining after sieve 83 x = [2, 1+2*find(sieve)]; # primes remaining after sieve
84 else 84 else
85 a = [2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, ... 85 a = [2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, ...
86 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, ... 86 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, ...
87 109, 113, 127, 131, 137, 139, 149, 151, 157, 163, 167, ... 87 109, 113, 127, 131, 137, 139, 149, 151, 157, 163, 167, ...
88 173, 179, 181, 191, 193, 197, 199, 211, 223, 227, 229, ... 88 173, 179, 181, 191, 193, 197, 199, 211, 223, 227, 229, ...
89 233, 239, 241, 251, 257, 263, 269, 271, 277, 281, 283, ... 89 233, 239, 241, 251, 257, 263, 269, 271, 277, 281, 283, ...
90 293, 307, 311, 313, 317, 331, 337, 347, 349]; 90 293, 307, 311, 313, 317, 331, 337, 347, 349];
91 x = a(a <= p); 91 x = a(a <= p);
92 endif 92 endif
93 93
94 endfunction 94 endfunction