Discussion :

[o2pubcy9]

Bonjour, je fais actuellement un programme mais je rencontre des problèmes lorsque je l’exécute.

bizarrement, je veux calculer l'intégrale entre 1 à 2 de la fonction carrée avec la méthode des rectangles alors j'ai écrit:

Code :

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

1
2
3
4
5
6
7
8
def carre (x) :
    return(x**2)
print(carre(2),carre(5),carre(14))
 
j=0
for(i) in range(10) :
    j=carre(1+i*(2-1)/10)*0.1+j
print(j)

mais le résultat est 2,185 alors que je devrais avoir 2.3333... soit 7/3
est ce normal svp ?

ensuite lorsque je veux calculer une autre intégrale on m'écris :

Code :

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

1
2
3
line 89,fonction
    y=(((y+v(0+i*h))**2))**1/2
OverflowError: (34, 'Result too large')

sur python j'ai écrit

Code :

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

1
2
3
4
5
6
def fonction(v,T,n) :
    h=(T)/n
    y=0.5*(v(0)+v(T))
    for i in range (1,n) :
        y=((1/T)*((y+v(0+i*h))**2))**1/2
    return (h*y)

Pouvez vous m'aider svp j'ai essayer de mettre import decimal mais ça ne change rien...

[VinsS]

Salut,

Tu es sur de ton algo ?

Ceci carre(1+i*(2-1)/10)*0.1+j c'est la même chose que cela carre(1+i*0.1) * 0.1 + j ce qui n'était peut-être pas ce que tu voulais faire.


Utilise les balises [code] pour ton code.

[o2pubcy9]

oui normalement c'est la méthode des rectangles pour calculer une intégrale entre deux bornes

savez-vous comment est-ce-que je pourrais faire pour résoudre le problème lié au code erreur overflow svp ?

[wiztricks]

Salut,

mais le résultat est 2,185 alors que je devrais avoir 2.3333... soit 7/3
est ce normal svp ?

7 / 3 est le résultat exact.
2,185 est une approximation obtenue avec un découpage en 10 intervalles.
Normalement, si vous voulez approcher la valeur exacte, il faut augmenter le nombre d'intervalles.

savez-vous comment est-ce-que je pourrais faire pour résoudre le problème lié au code erreur overflow svp ?

Je n'ai pas le temps d'essayer de comprendre ce que peut être y=(((y+v(0+i*h))**2))**1/2 par contre, vous pourriez réécrire cela avec des variables intermédiaires, afficher leurs valeurs avec "print" et vérifier que c'est cohérent avec ce que vous attendez (c'est ce qui s'appelle mettre au point son programme... et çà fait partie de l'exercice).

- W

[o2pubcy9]

bonjour, j'ai essayé de découper l'expression mais cela ne résoud toujours pas le problème

[o2pubcy9]

je dois faire le programme avec 10 intervalle donc c'est peut-être l'expression qui n'est pas bonne

Code :

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

1
2
3
for(i) in range(10) :
    j=carre(1+i*(2-1)/10)*0.1+j
print(j)

c'est bien l'intégrale entre 1 et 2 de la fonction carre avec la méthode des rectangles svp ???

[o2pubcy9]

Salut,

Tu es sur de ton algo ?

Ceci carre(1+i*(2-1)/10)*0.1+j c'est la même chose que cela carre(1+i*0.1) * 0.1 + j ce qui n'était peut-être pas ce que tu voulais faire.


Utilise les balises [code] pour ton code.

c'est bien ça

[o2pubcy9]

j'ai réussi à résoudre le problème lié à l'overflow...

[wiztricks]

Salut,

c'est bien l'intégrale entre 1 et 2 de la fonction carre avec la méthode des rectangles svp ???

Si vous écrivez:

Code :

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

1
2
3
for(i) in range(10) :
    j=carre(1+i*(2-1)/10)*0.1+j
print(j)

pas facile d'ajouter des intervalles pour obtenir un résultat plus précis.
Ne serait-il pas plus simple d'écrire:

Code :

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

1
2
3
4
5
6
7
8
9
a, b = 1, 2 # l'intervalle
h = 0.1 # la base des rectangles.
count = int( (b - a) / h)  # le nombre d'intervalles.
aire = 0
x = a
for i in range(count):
     aire += carre(x) * h
     x += h
print(aire)

- W

[o2pubcy9]

bonjour,
ce programme me donne approximativement le même résultat soit
2.1850000000000014 et 2.1850000000000005
du coup ce problème devrait plutôt venir de ma version de python...

[wiztricks]

Ce programme me donne approximativement le même résultat soit
2.1850000000000014 et 2.1850000000000005
du coup ce problème devrait plutôt venir de ma version de python...

Ce programme vous permet surtout de modifier la variable h pour avoir des intervalles plus petits et augmenter la précision du résultat. Essayez avec h = 0.01 puis 0.001 puis 0.0001,...

- W