" En el fitxer TERNSP.MTH es tracta de trobar terns pitag•rics primitius. "

" Concretament, la funci˘ TERNP(n,m,k) troba els primers k terns pitag•rics p~
rimitius "

" generats per 2ab, a^2-b^2, a^2+b^2, mitjan‡ant les k primeres parelles  (a,~
b), "

" a partir de (n,m), que compleixen a>b, mcd(a,b)=1 i a,b de diferent paritat~
. "

" ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo"

" ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo"

PRIMI(v):=v SUB 1>v SUB 2 AND (GCD(v SUB 1,v SUB 2)=1 AND MOD(v SUB 1,2)/=MOD~
(v SUB 2,2))

PRIMER_PRIMI(v):=ITERATE(IF(PRIMI(w),w,IF(w SUB 1>w SUB 2+1,[w SUB 1,w SUB 2+~
1],[w SUB 1+1,1])),w,v+[0,1])

KESIM_PRIMI(v,k):=IF(k>0 AND MOD(k)=0,IF(PRIMI(v),ITERATE(PRIMER_PRIMI(w),w,v~
,k-1),ITERATE(PRIMER_PRIMI(w),w,v,k)),"Error")

TERN(x,y):=[2*x*y,x^2-y^2,x^2+y^2,x,y]

TERNP(n,m,k):=VECTOR(IF(t=0,[" 2ab "," a^2-b^2 "," a^2+b^2 "," a "," b "],TER~
N((KESIM_PRIMI([n,m],t)) SUB 1,(KESIM_PRIMI([n,m],t)) SUB 2)),t,0,k)

TERNP(4,3,11)

;Simp(#12)
[[" 2ab "," a^2-b^2 "," a^2+b^2 "," a "," b "],[24,7,25,4,3],[20,21,29,5,2],[~
40,9,41,5,4],[12,35,37,6,1],[60,11,61,6,5],[28,45,53,7,2],[56,33,65,7,4],[84,~
13,85,7,6],[16,63,65,8,1],[48,55,73,8,3],[80,39,89,8,5]]