Найти сумму ряда с заданной точностью

#include "stdafx.h"
#include <iostream>
using namespace std;
 
int _tmain(int argc, _TCHAR* argv[])
{
    int kol=0,k,j,l,p,j1[10],k1[10],p1[10],l1[10];
    for(int i=0;i<100;i++)
    {
        kol=0;
        for(j=0;j<i;j++)
        {
            for(k=j+1;k<i;k++)
            {
                for(l=k+1;l<i;l++)
                {
                 for(p=l+1;p<i;p++)
                 {
                     if(j*j+k*k+l*l+p*p==i)
                     {
                         j1[kol]=j;
                         k1[kol]=k;
                         p1[kol]=p;
                         l1[kol]=l;
                        kol++;
                     }
                 }
                }
            }
        }
        if(kol>1)
        {
            for (j = 0; j<kol; j++)
                for (k = 0; k<kol; k++)
                    for (l = 0; l<kol; l++)
                        for (p = 0; p<kol; p++)
                            cout<<i<<" "<<j1[j]<<" "<<k1[k]<<" "<<l1[l]<<" "<<p1[p]<<endl;
            break;
        }
    }
    return 0;
}