2 réponses
c ******************resolution du system par gauss-siedle *****
do 50 kt=1,itmax
do 30 i=2,nx-1
smemb=ti(i-1)+fn*ti(i)+ti(i+1)
t(i)=(t(i-1)+t(i+1)+smemb)/f
30 continue
cccccccccccccccccc test de convergence cccccccccc48cccccccc
st=0.0
do 40 i=2,nx-1
st=st+t(i)
40 continue
sto=st
ecat=dabs(1.-sto/st)
if(ecat.lt.eps)go to 60
50 continue
ccccccc reinitialisation ccccccccccccccccccc60
60 print*,itemp,kt,ecat
do 70 i=2,nx-1
ti=t(i)
70 continue
ccccccccccccccccccccccccccccccccccccccccccc
100 continue
do 150 i=1,nx
write(1,*)x(i),t(i)
150 continue
print machine_data[i] - C syntax
stop
end
do 50 kt=1,itmax
do 30 i=2,nx-1
smemb=ti(i-1)+fn*ti(i)+ti(i+1)
t(i)=(t(i-1)+t(i+1)+smemb)/f
30 continue
cccccccccccccccccc test de convergence cccccccccc48cccccccc
st=0.0
do 40 i=2,nx-1
st=st+t(i)
40 continue
sto=st
ecat=dabs(1.-sto/st)
if(ecat.lt.eps)go to 60
50 continue
ccccccc reinitialisation ccccccccccccccccccc60
60 print*,itemp,kt,ecat
do 70 i=2,nx-1
ti=t(i)
70 continue
ccccccccccccccccccccccccccccccccccccccccccc
100 continue
do 150 i=1,nx
write(1,*)x(i),t(i)
150 continue
print machine_data[i] - C syntax
stop
end
Bonjour,
Le but des exercice c'est que tu apprenne quelque chose. Cours et tutoriels Fortran : https://fortran.developpez.com/
Le but des exercice c'est que tu apprenne quelque chose. Cours et tutoriels Fortran : https://fortran.developpez.com/