Discussion :

[mc4456]

Bonjour à tous,
j'ai un problème avec ce petit programme que j'ai créé pour un calcul de dérivée d'un polynome

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
program poly
implicit none 
real::derivpoly,x
real,allocatable::a(:),b(:)
integer::i,n
print*,'quel est le degré'
read*,n
allocate (a(n+1),b(n+1))
do i=1,n+1
  print*, 'que vaut a', i, '?'
  read*,a(i)
end do 
do i=1,n+1
  b(i)=(n+1.0-i)*a(i)*(x**(n-i))
end do
do i=1,n+1
  derivpoly=derivpoly+b(i)
end do
print*, 'la dérivée vaut', derivpoly
end program poly

Lorsque je veux éxecuter le programme, le problème vient à la ligne 14.

Merci d'avance pour votre aide.
Je débute dans le domaine.

[__dardanos__]

Salut,
Il faut mieux compiler avec l'option -Wall :

Code :

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

warning: 'x' may be used uninitialized in this function [-Wuninitialized]

[mc4456]

Salut,
merci pour ta réponse .
mais comment insérer cette instruction ?
J'utilise Silverfrost (Plato IDE), et je ne vois pas comment ajouter une option pour l'execution

[FR119492]

Salut!

Lorsque je veux éxecuter le programme, le problème vient à la ligne 14.

La raison est évidente: les valeurs de x ne sont pas définies.
Jean-Marc Blanc

[mc4456]

Je suis entièrement d'accord mais j'aimerai garder le polynôme en tant que fonction de x pour ensuite l'injecter dans un calcul d'intégrale.
C'est possible ?

Ou alors si j'ai envi que le programme affiche "2x" comme dérivée de (x²+3), en gros c'est pas possible ?!

[genteur slayer]

Attention: tu cherches à faire de la dérivation formelle.
c'est à dire de garder tes polynômes sous une forme sum_n (a_n x**n) ce que ne sait pas faire un ordinateur qui ne sait manipuler que des nombres.

Pour faire ce que tu veux il vaut mieux définir un tableau de coefficients (les a_n) tu modifie ce tableau car la dérivée d'un polynôme est un polynôme... et tu as ton polynome dérivée

tu aura une fonction d'évaluation

Code :

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

1
2
3
4
5
6
7
8
function evalpoly(A,N,X)
 real evalpoly,A(:),X
 integer N,I
 evalpoly=0.0
 DO I=1,N
  A(I)*X**(I-1)
 ENDDO
END

et la fonction qui dérive:

Code :

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

1
2
3
4
5
6
7
8
subroutine derivpoly(A,N,B)
 real A(:),B(:)
 INTEGER N,I
 B=0.0
 DO I=1,N-1
  B(I)=(I+1)*A(I+1)
 ENDDO
END