Discussion :

[dianbobo]

bonjour tout le monde
Je veux implémenter la récurrence de bonnet
NB : <<p(0,x) lire "p zéro de x">>
on a les données; p(0,x) = 1 et p(1,x) = x
p(n+1,x) = (2n+1)*x*p(n,x)/(n+1)-n*p(n-1,x)/(n+1)
je travaille dans l'intervale ]0,1[
pour cela j'ai crée 2 fonctions
la premiere fonction doit me generer un nombre aléatoire dans l'intervalle ]0;1[

Code :

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

1
2
3
4
5
6
double float_0_1(double a, double b)
{
    double x ;
    x = ( rand()/(double)RAND_MAX ) * (b-a) + a;
    return x ;
}

la deuxième fonction pour la récurrence :

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
void PolynomeBonnet()
{
    double L[]={};
    L[0]=1;
    int const N=10;
    double a,b,x;
    cout<<"la borne inferieur vaut"<<a<<endl;
    cin>>a;
    cout<<"la borne supérieure vaut"<<b<<endl;
    cin>>b;
    for (int i = 1;i<N; i++)
        {
            x = float_0_1(a,b);
            L[1]=x;
            L[i+1] = ((2i+1)/(i+1))*x*L[i]-(i/(i+1))*L[i-1];
        }
}

enfin dans ma fonction main() :

Code :

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

1
2
3
4
5
6
7
8
9
int main()
{
    double a,b,x;
    cout<<"la borne inferieur vaut"<<a<<endl;
    cin>>a;
    cout<<"la borne supérieure vaut"<<b<<endl;
    cin>>b;
return 0;
}

mais lorsque j’exécute ,j'arrive pas déjà à bien borner mon intervalle , par contre j'ai 2 nombres aléatoires qui se génèrent
alors que je voudrais juste qu'un nombre aléatoire x , se génère dans l'intervalle ]a;b[ , et qu'ensuite j’utilise ce nombre pour calculer le vecteur L

merci de votre aide

[3DArchi]

Salut,
Ce code ne compile pas. Je ne vois pas comment il pourrait s'exécuter...
double L[]={}; est une construction invalide en C++. Préfères

Code :

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

1
2
int const N=10;
std::vector<double> L(N+1);

ou

Code :

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

1
2
int const N=10;
std::array<double,N> L;

Sinon, pourquoi

Code :

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

1
2
x = float_0_1(a,b);
L[1]=x;

sont dans la boucle for ?

[dianbobo]

Bonjour
et merci à 3DArchi , pour ces éléments de réponses . j'ai déjà présenté le code comme ci-dessous:

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
#include <iostream>
#include<vector>
#include<cstdlib>
#include<ctime>
using namespace std;
double float_a_b(double a, double b)
{
    double x ;
    x = ( rand()/(double)RAND_MAX ) * (b-a) + a;
    return x ;
}
// calcul du polynome de bonnet
 
void PolynomeBonnet()
{
    int const N=10;
    std::vector<double> L(N+1);
        L[0]=1;
    double a,b;
    double x;
    x = float_a_b(a,b);
    L[1]=x;
    for (int i = 1;i<N; i++)
        {
           L[i+1] = ((2i+1)/(i+1))*x*L[i]-(i/(i+1))*L[i-1];
        }
}
int main()
{
    double a,b,dx;
        cout<<"la borne inferieur vaut a , entrer la valeur de a"<<endl;
cin>>a;
        cout<<"la borne supérieur vaut b , entrer la valeur de b"<<endl;
cin>>b;
cout<<"le pas de temps est dx"<<endl;
cin>>dx;
int N ;
    N = (b-a)/dx;
    cout<<"le nombre de pas de discretisation est N="<<N<<endl;
return 0;
}

mais j'ai une érreur à la ligne 25 .
error : cannot convert double__complex_' to ' double 'in assignment

merci de votre aide

[3DArchi]

2i => 2*i