négyzetösszeges feladat - szerencsére erre is van beépített függvényünk.

negyeshf[q_,k_]:=PowersRepresentations[q^2,k,2];

két adatsor, egyik: első 50 négyzetszám hányféleképpen bontható fel 2-5 négyzetszám összegére

data=Table[Length[negyeshf[q,i]],{i,1,5},{q,50}];

i=1??

másik pedig első 100, de most 4 négyzetszám összegére bontunk, így is sokáig fut.

data2=Table[Length[negyeshf[q,i]],{i,2,4},{q,100}];

i=2,3,4!!

esetleg nem csak a négyzetszámokra is megnézhetjük

(* data3=Table[Length[PowersRepresentations[w,i,2]],{i,2,5},{w,100,500}]; *)

rajzoljuk ki őket érdekességképpen

ListDensityPlot[data,Mesh-> None, InterpolationOrder-> 0,ColorFunctionScaling -> True, ColorFunction->"GrayYellowTones",AspectRatio -> 1/2]

ListDensityPlot[data2,Mesh-> None, InterpolationOrder-> 0,ColorFunction->"GrayYellowTones",AspectRatio -> 1/2]

ListDensityPlot[data3,Mesh-> None, InterpolationOrder-> 1,ColorFunction->"GrayYellowTones",AspectRatio -> 1/2]

ColorData["Gradients"]