Resolution d'un systeme d'equation

#include <stdlib.h>
#include <iostream>
#include <conio.h>
#include <iomanip>
using namespace std;

int main()
{

double matrice[2][3]={2,3,6,7,6,2};
double coefficient,x,y;
int i;
int j;
i=0;
j=0;

cout << "Résolution d'un système de 2 équations a trois inconnues\n";
cout << "Première équation, entrez en ordre respectif x,y et la constante\n";
cin >> matrice[0][0]>>matrice[0][1]>>matrice[0][2];
cout << "Seconde équation, entrez en ordre respectif x,y et la constante\n";
cin >> matrice[1][0]>>matrice[1][1]>>matrice[1][2];

coefficient=(-1*matrice[1][0]/matrice[0][0]);

i=0;
for(i;i<=2;i++)
{
matrice[1][i]=(coefficient*matrice[0][i])+matrice[1][i];               
}

y=matrice[1][2]/matrice[1][1];
x=(matrice[0][2]-(matrice[0][1]*y))/matrice[0][0];
cout << "X est egal a " << setw(10) << x << "\n";
cout << "Y est egal a " << setw(10) << y << "\n";
system("PAUSE");
return 0;
}