Problème de programmation sur visual C++

Bonjour,

j'ai un problème je dois écrire la résolution d'un système d'équation du second degré et je dois utilisé la dichotomie, je ne suis pas très expérimentée en langage C donc j'aurais vraiment besoin d'aide, j'ai les variables x1,x2 qui doivent etre relles ainsi que delta mais je dois aussi définir une variable PRECISION valant 0,005 et tant que (borne_sup - borne_inf>=PRECISION), r2= (borne_sup + borne_inf)/2 et en comparant r3(donc sensé être la carré de r2)à delta, je dois modifier la borne adéquate? Moi, j'ai écrit ce programme déjà, je ne sais même pas si c'est juste et en plus, il me met comme erreurs:
C:\Program Files\Microsoft Visual Studio\MyProjects\stephanieentrainement\entrainement\steph.c(72) : error C2059: syntax error : 'return'
C:\Program Files\Microsoft Visual Studio\MyProjects\stephanieentrainement\entrainement\steph.c(74) : fatal error C1004: unexpected end of file found
Error executing cl.exe.
pouvez vous m'aidez svp, je dois rendre ce devoir demain, si quelqu'un pouvait m'aider à voir si je ne fais pas fausse route et ou sont mes erreurs voilà bien deux jours que j'essaie de trouver la réponse au problème et ou sont mmes erreurs mais je ne trouve pas
.Merci d'avance de vos aides.

steph.obj - 2 error(s), 0 warning(s)
# include <stdio.h>
# include <stdlib.h>
# include <math.h>
# include <float.h>
# include <string.h>
# define PRECISION 0.005

float main ()

{
float a,b,c,delta,u,borne_inf,borne_sup,r3,r2;
double x1,x2;
borne_inf = 0;
borne_sup = FLT_MAX;
delta = (b*b) - (4*a*c);

printf("resolution d'un systeme d'equation du second degre\n");
printf("ax^2+ bx + c = 0\n");

printf("saisie du coefficient a\n");
scanf("%f", &a);
printf("saisie du coefficient b\n");
scanf("%f", & b);
printf("saisie du coefficient c\n");
scanf("%f", &c);
printf("%.2f est le discriminant de cette equation\n", delta);

if (delta<0)
{
printf("pas de solution reelle possible\n");
}
if (delta==0)
{
printf("une seule solution reelle possible, u\n");
u=(-b)/(2*a);
printf("%.2d\n", u);
}
if (delta > 0)
{
do
{
while (borne_sup - borne_inf>=PRECISION);
{
r2=(borne_sup + borne_inf)/2;
r3 = r2*r2;
if (r3<delta)
{
borne_sup = (float)r2;
break;
}
else
{ borne_inf = (float)r2;
break;
}

printf("deux solutions reelles possibles,x1 et x2\n");
x1=(float)((-b - (r2))/(2*a));
x2=(float)((-b + (r2))/(2*a));
printf("%.2f\n", x1);
printf("%.2f\n", x2);
break;
}
}

return borne_inf;
}