salut
comment ouvrir un fechier dans un programme de fortran![]()
salut
comment ouvrir un fechier dans un programme de fortran![]()
Salut!
En utilisant l'instruction open.comment ouvrir un fechier dans un programme de fortran
Jean-Marc Blanc
Tout est là, il suffit d'aller lire la FAQ ou n'importe quel poly de fortran qu'il soit 77/90/95/2003 disponible gratuitement sur le net :
http://fortran.developpez.com/faq/
(ou encore à peut près n'importe quel post du forum fortran traitant de l'ouverture de fichier(s))
Bonne journée,
Marlan
salut j'ai essayéé de faire ce probleme à 2D equation de chaleur mais le probleme ilya des erreurs et que il faut ilya instrument de sujet resultat mais je ne peut pas faire un fichier
voicii le probleme
" program master 2D program msmf
implicit none
integer i,j,n,nmax,k,m,mmax
parameter (nmax=30,mmax=30)
real a(nmax,mmax),b(nmax,mmax),c(nmax,mmax),d(nmax,mmax),p(nmax,mmax),t0(nmax,mmax)
real q(nmax,mmax),t(nmax,mmax),tan(nmax,mmax),tb
do I =1,4
T(I,J)=00.00
enddo
n=4
a(1,1)=20.0
a(2,2)=30.0
a(3,3)=30.0
a(4,4)=40.0
b(1,1)=10.0
b(2,2)=10.0
b(3,3)=10.0
b(4,4)=0.0
c(1,1)=0.0
c(2,2)=10.0
c(3,3)=10.0
c(4,4)=10.0
d(1,1)=500+10*T(5,5)
d(2,2)=10*T(6,6)+500
d(3,3)=10*T(7,7)+500
d(4,4)=2500+T(8,8)*10
do j=1,4
c print*,a(j),b(j),c(j),d(j)
enddo
p(1,1)=b(1,1)/a(1,1)
q(1,1)=d(1,1)/a(1,1)
do j=2,4
p(i,j)=b(i,j)/(a(i,j)-(c(i,j)*p(i,j-1)))
q(i,j)=(d(i,j)+(c(i,j)*q(i,j-1)))/(a(i,j)-(c(i,j)*p(i,j-1)))
enddo
do j=4,1,-1
t(i,j)=(p(i,j)*t(i,j+1))+q(i,j)
print*,j,t(i,j)
enddo
a(5,5)= 30.0
a(6,6)= 40.0
a(7,7)= 40.0
a(8,8)= 50.0
a(9,9)= 20.0
b(5,5)= 10.0
b(6,6)=10.0
b(7,7)=10.0
b(8,8)=0.0
b(9,9)= 10.0
c(5,5)=0.0
c(6,6)=10.0
c(7,7)=10.0
c(8,8)=10.0
c(9,9)=0.0
d(5,5)=10*T(9,9)+10*T(1,1)
d(6,6)=10*T(2,2)
d(7,7)=10*T(3,3)
d(8,8)=2000+10*T(4,4)
d(9,9)=10*T(5,5)
do j=5,9
c print*,a(i,j),b(i,j),c(i,j),d(i,j)
enddo
p(5,5)=b(5,5)/a(5,5)
q(5,5)=d(5,5)/a(5,5)
do j=5,9
p(i,j)=b(i,j)/(a(i,j)-(c(i,j)*p(i,j-1)))
q(i,j)=(d(i,j)+(c(i,j)*q(i,j-1)))/(a(i,j)-(c(i,j)*p(i,j-1)))
enddo
do j=9,5,-1
t(i,j)=(p(i,j)*t(i,j+1))+q(i,j)
print*,j,t(i,j)
enddo
do I=1,9
S(1)=0.0
S(2)=0.0
IF(I.EQ.1)
DO I=I+1,N
S(2)=S(2)+a(I,J)*T(J)
ENDDO
T(I,J)=b(1,1)/a(1,I)-S(2)/a(1,I)
ELSE
DO I=I,I-1
S(1)=S(1)+a(I,J)T(J)
DO J=I+1,N
S(2)=S(2)+a(I,J)*T(J)
ENDDO
T(I)=(1/a(I,I))*(b(I)-S(1)-S(2))
END
"
ok merciii pour votre aide ^_^
Partager