Discussion :

[nailita]

Bonjour tout le monde,
J'ai essayé de faire un petit programme en langage C, calculant les coefficients qui nous permettent d'avoir une fonction avec l'interpolation de newton, mais ça n'a pas marché.
Ci-dessous le code du programme, et je ne sais pas ou est l'erreur.
J'obtiens comme résultat ce qui suit: A[1]=29817731, A[2]=0, et ce n'est pas le résultat que je devrai avoir

Merci d'avance pour votre aide.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
#include <stdio.h>
 
void newton(int n)
{
       int i,j;
     double X[20];
      double Y[20]; 
 
        float F[n][n];
        float A[20];
for(i=1;i<n;i++)
{                  
  F[i][0]=Y[i];
  }
 
 
          for(i=1;i<n;i++)
          {                   
                 for(j=1;j<i;j++)
                                 {
                                        F[i][j]=(F[i][j-1]-F[i-1][j-1])/(X[i]-X[i-j]);
 
                                 }
          }
 
 
 
          for(i=1;i<n;i++)
          {
                          A[i]=F[n][i]; 
                         }
 
 
                 for(i=1;i<n;i++)
          {
 
                   printf("A[%d]=%d \n" , i,A[i]);                      
}
}
 
 
 
int main() {
            int i,n,j;        
 float F[20][20];
 double X[20];
  double Y[20];
 
 
 printf("donner le degre du polynome \n");
 scanf("%d",&n);
 
 
 
 for(j=1;j<=n;j++)
{
printf("element X[%d]= ",j);
scanf("%d",&X[j]);
}
 
 
 
for(i=1;i<=n;i++)
{
printf("element Y[%d]= ",i);
scanf("%d",&Y[i]);
}
 
 
 
newton(n);
 
 system("pause");
 return 0;
 
 
}

[Invité(e)]

Bonjour,

Les fonctions main et newton travaillent sur deux tableaux différents : Dans newton, X et Y valent ne sont pas initialisé.

Il faudrait ajouter en paramètre de la fonction newton les tableaux X et Y.

EDIT :
Il faut aussi faire attention à l'utilisation de scanf; et ne pas ignorer les warnings du compilateur

PS : merci de penser aux balises code (#)

[nailita]

Bonjour,

Les fonctions main et newton travaillent sur deux tableaux différents : Dans newton, X et Y valent ne sont pas initialisé.

Il faudrait ajouter en paramètre de la fonction newton les tableaux X et Y.

EDIT :
Il faut aussi faire attention à l'utilisation de scanf; et ne pas ignorer les warnings du compilateur

PS : merci de penser aux balises code (#)

Je les ai ajouter, et j'ai toujours le même problème.

Voici ma nouvelle fonction newton:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
void newton(int n, float X[], float Y[])
{
       int i,j;
 
 
        float F[30][30];
        float A[30];
        float k;
        float v;
 
 
for(i=1;i<n;i++)
                {              
                 F[i][0]=Y[i];
                 }
 
 
          for(i=1;i<n;i++)
          {                   
                 for(j=1;j<i;j++)
                                 {
                                        k=F[i][j-1]-F[i-1][j-1];
                                        v=X[i]-X[i-j];
             F[i][j]=k/v;
 
                                 }                              
          }
 
          for(i=1;i<n;i++)
          {
                          A[i]=F[n][i]; 
           }
 
 
                 for(i=1;i<n;i++)
          {
 
                   printf("A[%d]=%d \n" , i,A[i]);                      
}
}

[Invité(e)]

Je les ai ajouter, et j'ai toujours le même problème.

Voici ma nouvelle fonction newton:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
for(i=1;i<n;i++)
{                   
    for(j=1;j<i;j++)
    {
        k=F[i][j-1]-F[i-1][j-1];
        v=X[i]-X[i-j];
        F[i][j]=k/v;
    }                              
}
 
for(i=1;i<n;i++)
{
    A[i]=F[n][i]; 
}

Dans la première boucle for, i va de 1 à n-1 pour assigner F[i][...]
Dans la dernière, tu utilises les valeurs de F[n][...], or F[n][...] ne sont jamais assignées (tu t'arrêtes à n-1).

Aussi, F[0][0] n'est jamais assigné, or il est utilisé (dans le premier calcul de k).