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 30 31 32 33 34 35 36 37 38 39 40
| program pic
implicit none
integer :: i,j,k
real :: energie(40),energie1(40),tab(28,28)
real:: somme(40),a(40),b(40),c(40),d(40)
open(1,file='energie.txt')
open(2,file='pic.txt')
open(3,file='pic1.txt')
do i=1,29
read(1,*) energie(i)
c(i)=energie(i) -2
d(i)=energie(i) +2
end do
do i=1,28
do j=1,28
tab(i,j)=energie(j) + energie(i+1)
write(2,*) i,j,energie(i+1), energie(j),tab(i,j)
end do
end do
close(2)
open(2,file='pic.txt',position='rewind')
do i=1,28
read(2,*) a(i),b(i),energie(i), energie1(i),somme(i)
if (energie(i) .GE. c(i)) .AND. (energie(i) .LE. d(i)) then
write(3,*) somme(i)
end if
end do
close(1)
close(2)
close(3)
END |
Partager