blob: b9dae24ebd7af8aeb22b4521c5f91a2f03c62e54 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
|
//Chapter 04:Number Theory and Cryptography
clc;
clear all;
function primefactors(n)
while modulo(n,2) == 0 //To print all the 2s that divide input
disp('2')
n=n/2
end
for i=3:2:sqrt(n)//increment by 2 so as to obtain odd numbers only
while modulo(n,i)==0
disp(i)
n=n/i
end
end
if(n>2) then //to check for prime number
disp(n)
end
endfunction
n1=7007
mprintf("Prime factors of %d are:",n1)
disp(primefactors(n1))
|
