-1
sur la compilation de ce code à l'aide Code :: Blocks 16.0.1 sous Windows, l'erreur suivante apparaîtErreur dans boucle do: Estimé une parenthèse droite dans l'expression à (1)
a(i) = (1/(delx)^2)-(1/x-1)
1
Error: Expected a right parenthesis in expression at (1)
b(i) = ((-2/((delx)^2)) - bi*ce)
1
Error: Expected a right parenthesis in expression at (1)
c(i) = (1/(delx)^2) + (1/(2*(delx)*(x-1)))
1
Error: Expected a right parenthesis in expression at (1)
En ajoutant au & fin de chacune des 3 déclarations, ces 2 derniers deviennent bien, mais la 1ère reste le même
program HW0
implicit none
integer, parameter :: nx=101
integer, parameter :: ce = 100
real*8, parameter :: bi = 10.d0
real*8 :: delx
integer :: i
real*8, dimension(nx) :: t, x
real*8, dimension(nx) :: a, b, c, d
i = 1
delx = 1/(nx-1)
a(1) = 0; b(1) = 1; c(1) = 0; d(1) = 0
a(nx) = -(1/(delx))
b(delx) = (1/(delx)) + bi
c(nx) = 0; d(nx) = 0;
do i = 2, (nx - 1)
a(i) = (1/(delx)^2)-(1/x-1)
b(i) = ((-2/((delx)^2)) - bi*ce)
c(i) = (1/(delx)^2) + (1/(2*(delx)*(x-1)))
d(i) = 0
enddo
do i = 2, nx
b(i) = b(i) - c(i-1)*(a(i)/b(i-1)); d(i) = d(i) - d(i-1)*(a(i)/b(i-1))
enddo
t(nx) = d(nx)/b(nx)
do i = (nx-1), 1, -1
t(i) = (d(i) - c(i)*t(i+1))/b(i)
enddo
stop
end program HW0
Veuillez également faire attention avec "delx = 1/(nx-1)", qui est évalué à 0 (parce que 1/100 -> 0 comme division entière). Nous pouvons obtenir une valeur correcte comme delx = 1.0d0/(nx-1) (ou 1.0_dp/... etc etc). – roygvib
Oh, ça m'a échappé. Merci pour l'aide! –