## 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