Forum de mathématiques - Bibm@th.net

SALUT TOUT LE MONDE
dans ce programme ,j ai une erreur dans l'affichage de resultat de determinant
je vois pas ou est le probleme

# include <iostream>
# include <cmath>
using namespace std;
/*lecture de deux matrices carrées*/
void lecture (double X[][10],int a)
{
     int I,J;   
        for (I=0;I<a;I++)
        for (J=0;J<a;J++)
    {
    cout<<"donner votre element   "<<I<<" , "<<J<<endl;
    cin>>X[i][J];
    }
}
/* affichage de la matrice*/
void affichage(double X[][10],int a)
{
     int I,J;
        cout<<"votre matrice est la suivante:"<<endl;
        for (I=0;I<a;I++)
{
        for (J=0;J<a;J++)
        cout<<X[i][J]<<"    ";
        cout<<endl;
}    
}   
void produit(double X[][10],double Y[][10],double Z[][10],int a)
{
int I,J,K;     
//initialisation de la matrice C
        for (I=0;I<a;I++)
        for (J=0;J<a;J++)
         Z[i][J]=0;      
// Calcul de la matrice C=A*B
        for (I=0;I<a;I++)
        for (J=0;J<a;J++)
        for (K=0;K<a;K++)
         Z[i][J]=Z[i][J]+X[i][K]*Y[K][J];   
}
void puissance (double X[][10],double CB[][10],int a)
{
     double Ck[10][10];
     int i,j,k,L;
     cout<<"donner la valeur de k pour caculer C^k"<<endl;
     cin>>k;
    for(i=0;i<a;i++)
    for(j=0;j<a;j++)
        Ck[i][j]=X[i][j];

     for(L=1;L<k;L++)
      {produit(X,Ck,CB,a);
        for(i=0;i<a;i++)
        for(j=0;j<a;j++)
        Ck[i][j]=CB[i][j];   
} 
}

/*Cette fonction nous est utile dans le calcul du determinant, elle supprime dans
  la matrice ma les lignes l et c puis enregistre les resultats dans mb*/
void det_aux(double ma[][10],double mb[][10],int l,int c)
{
int i,j,d,e=0;
    for(i=0;i<10;i++)
    {
   d=0;
        if(i!=l)
        {
      for(j=0;j<10;j++)
                if(j!=c)
                {
            mb[e][d]=ma[i][j];
            d++;
                }
            e++;
        }
    }
}
/*******************************************************************************/
/*Cette fonction teste si un entier est multiple de 2,
cette fonction nous sert dans le calcul du determinant*/
double expo(int n)
{
if(!(n%2)){return (1);}
return (-1);
}
/********************************************************************************/
/* Cette fonction calcul le determinant d'une matrice de maniere recursive*/
double determinant(double m[][10],int l)
{
int i;
double m2[10][10],x=0;
if(l==1)
{
        return (m[0][0]);
}
for(i=0;i<l;i++)
{
det_aux(m,m2,i,0);
x=x+(expo(i)*m[i][0]*determinant(m2,(l-1)));
}
return (x);
}
 
int main()
{
    int n;
    double A[10][10],B[10][10],C[10][10],CP[10][10];
    cout<<"donnez la taille des matrices"<<endl;
    cin>>n;
    cout<<"lecture de la matrice A"<<endl;
    lecture(A,n);
    cout<<"lecture de la matrice B"<<endl;
    lecture(B,n);
    produit(A,B,C,n); 
    cout<<"affichage des elements de la matrice A "<<endl;
     affichage( A,n);
     cout<<"affichage des elements de la matrice B "<<endl;
     affichage( B,n);
     cout<<"affichage des elements de la matrice C "<<endl;
     affichage(C,n);
     puissance (C,CP,n);
     cout<<"la matrice de puissance k "<<endl;
     affichage(CP,n);
cout<<"Le determinant de notre matrice C est :"<<endl;
determinant(C,n);
system("PAUSE");
return 0;
}

Bonsoir
j'ai travaller sur un prgramme qui fait des opérations sur une fraction mais il y a un défauts dans ce dernier que je trouve pas
merci de m aider

/* on utilise la notion de structure et fonction pour faire des opérations sur une fraction*/
# include <iostream>
# include <cmath>
using namespace std;
struct fract
{
       double N,D;};
       void saisie_fract (fract &F)
       {
            cout<< "donnez la valeur de nominateur de la fraction "<<endl;
            cin>>F.N;
            cout<< "donnez la valeur de dénominateur de la fraction "<<endl;
            cin>>F.D;
            cout<<endl;
            }
       void affichage_fract (fract F)
       {
            cout<<" votre fraction finale  est "<<F.N<<"/"<<F.D<<endl;
            cout<<endl;
            }     
       double add_fract (fract A,fract B)
       {
              fract C;
              C.N=A.N*B.D+B.N*A.D;
              C.D=A.D*B.D;
              return(C.N,C.D);
              }
       double diff_fract (fract A,fract B)
       {
              fract C;
              C.N=((A.N*B.D)-(B.N*A.D));
              C.D=A.D*B.D;
              return(C.N,C.D);
              } 
       double pro_fract (fract A,fract B)
       {
              fract C;
              C.N=A.N*B.N;
              C.D=A.D*B.D;
              return(C.N,C.D);
              }   
        double rapp_fract (fract A,fract B)
       {
              fract C;
              C.N=A.N*B.D;
              C.D=A.D*B.N;
              return(C.N,C.D);
              } 
        int main()
        {
            fract a,b;
            fract plus,mois,fois,sur;
            cout<<" saisie de la fraction F1 "<<endl;
            saisie_fract (a);
            cout<<" saisie de la fraction F2 "<<endl;
            saisie_fract (b);
            plus  = add_fract (a,b);
            moins = diff_fract (a,b);
            fois = pro_fract (a,b);
            sur =  rapp_fract (a,b);
            cout<<"votre fraction addition vaut   "<<"="<<endl;
            affichage_fract (plus);
            cout<< " FIN "<<endl;
cout<<"votre fraction difference vaut   "<<"="<<endl;
            affichage_fract (moins);
            cout<< " FIN "<<endl;
cout<<"votre fraction produit vaut   "<<"="<<endl;
            affichage_fract (fois);
            cout<< " FIN "<<endl;
cout<<"votre fraction rapport vaut   "<<"="<<endl;
            affichage_fract (sur);
            cout<< " FIN "<<endl;
system("PAUSE");
return 0;
}

Bonjour tous monde
je viens juste de commencer à étudier la méthode de Jacobi et j'arrive à toucher son principe mas je trouve des exemples pour bien assimiler la méthode et fare face aux difficultés de calcul;

alors j 'ai besoin de vos aides
merci

Bonsoir tout le monde
j ai besoin d'aide pour:
ecrire un progrmme en C++ qui calcul
1_______le produit de ces deux matrices A(N,N) et B(N,N) : C=A*B
2_______ C^k tel que k>0
3_______det C
4_______C^(-1) c à d l'inverse de C
et tout cela en utilisant le principe de fonction pour mettre le tout dans un seul programme
merci d'avance