## Wolfdieter Lang, Sep 25 2012 ## Maple 13 program for the radical (also known as squarefree kernel) of a product of numbers ## a,b and c=a+b. ## copy and save this under some name, e.g., radabc.txt, and write in your ## maple program a line: read "radabc.txt ", then you can use rabc(a,b) ## in your maple program. rabc:=proc(a,b) local c,sqfkn: c:=a+b: sqfkn:=proc(n) ##(adapted from OEIS A007947) local j,t1,nrt1,t2: if n=1 then 1 else t1:=ifactors(n)[2]: nrt1:=nops(t1): t2 := product('op(1,t1[j])','j'=1..nrt1): fi: end proc: ## sqfkn sqfkn(a*b*c) end proc: ## e.g., rabc(2,8) = 10