First, some examples :
17 * 29 * 41 * 149 * 113 = 340326281
= 2984^2 18205^2 = 18220^2 2891^2
509 * 853 * 769 = 333882113
= 11528^2 14177^2 = 18272^2 127^2
The Main programme is the subfunction Test_Random_2_Squares within the function : Runs_Of_Test_Random_2_Squares_1.m which tries to find 2 sets of arrays of positive integers such that the corresponding sums of squares of the numbers in the 2 arrays equals the product of r no of Prime Factors.
The function : find_the_2_Squares_from_the_Sum ( Sum ) forms the heart of this suite, which can also be independently used.
Whereas, currently (Aug 2005), the top function : Runs_Of_Test_Random_2_Squares_1.m
deals only with prime factors of the form (4*n 1), the heart-function : fnd_the_2_Squares_from_the_Sum.m can still be independently used to find sum of squares - even if the factors include primes of type (4*n-1), or even 2.
The various steps of the subfunction Test_Random_2_Squares.m are :
1) This overall pgm first calls Prime_Fctrs_Of_Sum_Of_2_Sqrs(N)
to obtain the list of primes which are <= N, and which are of the form (4n 1).
2) Next, it randomly chooses "r" nos from the output primes_N_4nPlus1 of the above call.
3) Next, it calls [ Match_Rnded_Array_1, Match_Rnded_Array_2 ] =
to get the 2 sets of arrays.