Создать СЛАУ размерности 1000х1000, чтобы система имела только одно решение

// matrica.cpp : Defines the entry point for the console application.
//
 
#include "stdafx.h"
#include <iostream>
#include <iomanip>
#include<time.h>
#include<stdlib.h>
#include <fstream>
 
#define n 1000
 
using namespace std;
 
int _tmain(int argc, _TCHAR* argv[])
{
 
 
    srand(time(NULL));
    double tab[n][n];
 
    /* */
    double B[n];
    for(int i=0;i<n;i++) {
        B[i]=rand()%10;
    }
 
    ofstream f1;
 
        for(int i=0;i<n;i++) {
        f1.open("e:\\b_matrix_old.txt", ios_base::app);
        f1<<B[i]<<" ";
        f1.close();
        }
 
    /* */
 
    int i=0, j=0;
 
    for (i=0; i<n; i++)
    {
        for (j=0; j<n; j++) {
            if (i==j) 
            {tab[i][j]=1;}               
            else 
            {tab[i][j]=0;}
        }
    }
 
 
        double K;
        int N2,N1;
        
        for(int i=0;i<10000;i++) {
        K=(1+rand()%1);
        N1=rand()%n;
        N2=rand()%n;
 
        for(int ii = 0; ii < n; ++ii)
        {
            if(N1 != N2) {
                tab[N2][ii] += tab[N1][ii] * K;
                if(ii==0) {B[N2] = B[N2] + B[N1]*K;}
            }
        }
        }
 
 
        ofstream f;
 
        for(int i=0;i<n;i++) {
        f.open("e:\\b_matrix.txt", ios_base::app);
        f<<B[i]<<" ";
        f.close();
        }
 
 
        for(int i=0;i<n;i++) {
        for(int j=0;j<n;j++) {
            f.open("e:\\matrix.txt", ios_base::app);
            f<<tab[i][j]<<" ";
            f.close();
        }
            f.open("e:\\matrix.txt", ios_base::app);
            f<<"\n";
            f.close();
    }
 
    system("pause");
    return 0;
}