Fortran

Bonjour,
j'ai essayé de faire un petit test sur fortran "résolution d'une equation de 2ieme degré" mais ça marche pas, voici ce que j'ai écrit

program equation
real a , b , c , x1 , x2 , d
do
print*, 'entrez la valeur de a'
read*, a
print*, 'entrez la valeur de b'
read*, b
print*, 'entrez la valeur de c'
read*, c
end do
d = b*b - 4*a*c
print*, 'la valeur de delta est égale :' ,d
if (d=0)
then
x1 = -b/(2*a)
print*, 'léquation admet une seule solution x=' ,x1
endif
if (d>0)
then
x1 = (-b-sqrt(d))/(2*a)
x2 = (-b+sqrt(d))/(2*a)
print*, 'léquation admet deux solutions' ,x1 ,x2
endif
if (d<0)
then
print*, 'léquation na pas de solution dans R'
endif
end

Merci de m'aidez.