Discussion :

[feynman]

Bonjour,
J'ai à écrire un programme pour résoudre un système d'équations linéaires à deux inconnues; la seule chose que je n'arrive pas à faire est de faire apparaître le système (avant de le résoudre) sur l'écran.
Voilà mon programme, merci de me signaler quoi ajouter

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
implicit none
real u1,u2,v1,v2,x,y,w1,w2,delta_1,delta_2,delta_3,delta_4
print* ,'entrer le nombre u1:'
read*, u1
print*, 'entrer le nombre u2:'
read*, u2
print* ,'entrer le nombre v1:'
read*, v1
print*, 'entrer le nombre v2'
read*, v2
print*,' entrer le nombre w1'
read*, w1
print* ,'entrer le nombre w2:'
read*, w2
delta_1=v1*w2-v2*w1
delta_2=u2*v1-u1*v2
delta_3=w2*u1-u2*w1
delta_4=v2*u1-v1*u2
x=delta_1/delta_2
y=delta_3/delta_4
 
print*,x 
print*,y
end program jjjj

[FR119492]

Salut!
Quelques commentaires sur ton problème:

1) Tu utilises la méthode de Cramer que, malheureusement on continue d'enseigner dans de nombreuses écoles. Pour un système de deux équations à deux inconnues, ça peut encore aller, mais pour des systèmes plus gros, le temps de calcul devient prohibitif (n factorielle pour Cramer et n**3 pour les "bonnes" méthodes. Si tu tombes sur des systèmes plus gros, je te recommande la bibliothèque gratuite LinPack, que tu trouveras sur le site www.netlib.org .

2) Tu sembles ignorer que le Fortran permet de traiter des tableaux, correspondant aux matrices en mathématiques.

3) Tu as calculé deux fois la même chose: delta_2 et delta_4 sont égaux au signe près.

4) Pourquoi n'utilises-tu pas le paramètre Advance='No' dans ton dialogue interactif?

5) Pourquoi n'utilises-tu pas les formats pour l'affichage du système et de la solution?

Voici comment je programmerais ton problème (sous réserve des inévitables fautes de frappe et autres, car je n'ai pas eu le temps de tester!) :

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
Implicit None
Real*8 A(2,2),B(2),X(2)
Write (*,'(A11)'),Advance='No') 'Tapez A11: '
Read (*,*) A(1,1)
Write (*,'(A11)'),Advance='No') 'Tapez A12: '
Read (*,*) A(1,2)
Write (*,'(A10)'),Advance='No') 'Tapez B1: '
Read (*,*) B(1)
Write (*,'(A11)'),Advance='No') 'Tapez A21: '
Read (*,*) A(2,1)
Write (*,'(A11)'),Advance='No') 'Tapez A22: '
Read (*,*) A(2,2)
Write (*,'(A10)'),Advance='No') 'Tapez B2: '
Read (*,*) B(2)
 
Write (*,'(/A18/)') 'Système à résoudre'
Do I=1,2
  Write (*,'(2EN12.3,5X,EN12.3)') (A(I,J),J=1,2),B(I)
End Do
 
Delta=A(1,1)*A(2,2)-A(1,2)*A(2,1)
X(1)=(B(1)*A(1,2)-B(2)*A(2,2))/Delta
X(2)=(A(1,1)*B(2)-A(2,1)*B(1))/Delta
 
Write (*,'(/A19/)') 'Solution du système'
Write (*,'(2EN12.3)') (B(J),J=1,2)
 
Stop
End

Bonne chance
Jean-Marc Blanc

[feynman]

Merci pour ta reponse,
mais excuse moi, je ne suis qu un debutant en fortran, et j apprends les choses pas a pas

[feynman]

D'autre part tu écris :

Code :

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

Write (*,'(A11)'),Advance='No') 'Tapez A11: '

Je vois 3 parenthèses fermées et 2 ouvertes, alors où est-ce que je dois mettre la 3ème parenthèse pour ne plus avoir d'erreurs ?

[FR119492]

Je t'avais bien dit de te méfier de mes fautes de frappe. L'instruction devrait être:

Code :

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

Write (*,'(A11)',Advance='No') 'Tapez A11: '