Méthode de l'analyse numérique en pascal

Bonjour hakou28, je crois que tu postais ce programme pour aider dino (merci pour lui), cependant il m'a semblé bon d'y apporter quelques modifications...

J'ai principalement aérer le code pour bien voir les structures, et après je me suis apercu de nombreuses répétitions dans le code (surtout dans la partie affichage de ton programme principal) j'ai donc mis ensemble ce qui pouvait l'être pour finalement réussir à enlever plus de 100 lignes de code inutile !

Après j'ai apporté quelques améliorations ici ou là, notamment dans l'utilisation des variables globales qu'il faut à tout prix éviter...
Et puis tant que j'y étais j'ai aussi corriger les fautes de français et en particulier lors de l'utilisation du vocabulaire mathématiques.

Par contre je n'ai pas vérifié le cœur du programme (les calculs) j'ai considéré que c'était juste, enfin bref, voici ce que ça donne...

program ziouche_interface;

uses crt;

function f(x:real):real;
begin
f:=(sqr(x)+2)*x+1;
end;

function newt(x:real):real;
begin
newt:=(2*sqr(x)*x-1)/(3*sqr(x)+2);
end;

function lag(x:real):real;
begin
lag:=(-sqr(x)+x-1)/(sqr(x)-x+3);
end;

function pf(x:real):real;
begin
pf:=(-sqr(x)*x-1)/2;
end;

procedure entree(var xmin,xmax,seuil:real);
begin
textcolor(14); gotoxy(1,10);
write('Donnez une valeur < -1 pour xmin  : '); readln(xmin);
write('Donnez une valeur > 0  pour xmax  : '); readln(xmax);
write('Donnez une valeur pour le seuil   : '); readln(seuil);
End;

procedure baly(a,xmax,seuil:real);
var y1,y2,xsuit:real;
begin
while a<xmax do
      begin
      y1:=f(a);
      xsuit:=a+seuil;
      y2:=f(xsuit);
      if (y1*y2)<0 then writeln('La solution reelle est dans l''intervalle [',a:0:5,',',xsuit:0:5,']');
      a:=a+Seuil;
      end;
end;

procedure dicho(a,b,seuil:real);
var mil:real;
begin
writeln ('La solution reelle est dans l''intervalle ');
writeln('[',a:0:5,',',b:0:5,']');
while abs(b-a)>=seuil do
      begin
      mil:=(a+b)/2 ;
      if (f(a)*f(mil))<=0 then b:=mil
                          else a:=mil;
      writeln('[',a:0:5,',',b:0:5,']');
      end;
end;

procedure Newton(x1,seuil:real);
var x0,y:real;
begin
x0:=x1;
while abs(x0-newt(x0))>=seuil do
      begin
      y:=newt(x0);
      x0:=y;
      end;
writeln('La solution reelle est environ ',x0:0:9);
end;

procedure Lagrange(x1,seuil:real);
var x0,y:real;
begin
x0:=x1;
while abs(x0-lag(x0))>=seuil do
      begin
      y:=lag(x0);
      x0:=y;
      end;
writeln('La solution reelle est environ ',x0:0:9);
end;

procedure Pointfixe(x1,seuil:real);
var x0,y:real;
begin
x0:=x1;
while abs(x0-pf(x0))>=seuil do
      begin
      y:=pf(x0);
      x0:=y;
      end;
writeln('La solution reelle est environ ',x0:0:9);
End;

procedure fausep(a,b,seuil:real);
var mil:real;
begin
writeln('La solution reelle est dans l intervalle ');
writeln('[',a:0:5,',',b:0:5,']');
while abs(b-a)>=seuil do
      begin
      mil:=a-((f(a)/(f(b)-f(a)))*(b-a)) ;
      if (f(a)*f(mil))<=0 then b:=mil
                          else a:=mil;
      writeln('[',a:0:5,',',b:0:5,']');
      end;
end;

procedure Menu;
begin
clrscr;
textcolor(10); gotoxy(27,2);  writeln('ZIOUCHE Abdelhak D0307010') ;
textcolor(30); gotoxy(22,4);  writeln('Le mini projet d''Analyse Numerique ');
textcolor(10); gotoxy(12,6);  writeln('Recherche de la solution reelle de l''equation xü+2x+1=0');
textcolor(9);  gotoxy(18,8);  writeln('*****************************************');
textcolor(2);  gotoxy(20,9);  writeln('Choisissez comment faire le calcul');
textcolor(1);  gotoxy(20,11); writeln('-1- par balayage');
textcolor(8);  gotoxy(20,12); writeln('-2- par dichotomie');
textcolor(3);  gotoxy(20,13); writeln('-3- par la methode de Newton');
textcolor(4);  gotoxy(20,14); writeln('-4- par la methode du point fixe');
textcolor(5);  gotoxy(20,15); writeln('-5- par la methode de fausse position');
textcolor(6);  gotoxy(20,16); writeln('-6- par la methode de Lagrange');
textcolor(7);  gotoxy(20,17); writeln('-7- sortir');
textcolor(9);  gotoxy(18,18); writeln('*****************************************');
textcolor(15);
end;

procedure EnTete;
begin
clrscr;
textcolor(26); gotoxy(27,4); writeln('ZIOUCHE Abdelhak D0307010') ;
textcolor(14); gotoxy(22,6); writeln('Le mini projet d''Analyse Numerique ');
textcolor(26); gotoxy(15,8); writeln('Recherche de la solution reelle de l''equation xü+2x+1=0');
end;

var s:integer; xmin,xmax,seuil:real;
begin repeat
      Menu;
      readln(s);
      EnTete;
      entree(xmin,xmax,seuil);
      writeln;
      case s of
           1:baly(xmin,xmax,seuil);
           2:dicho(xmin,xmax,seuil);
           3:newton(xmin,seuil);
           4:pointfixe(xmin,seuil);
           5:Fausep(xmin,xmax,seuil);
           6:Lagrange(xmin,seuil);
           end;
      textcolor(4); writeln;
      writeln('Veuillez taper sur la touche Entree pour choisir une autre methode');
      readln;
      until s=7;
end.

Je viens seulement de commence le cours de maths sur la réduction de matrice, alors j'ai pas le code qui correspond, et si je le fais (c'est pas sûr) ce ne sera probablement pas en Pascal, mais en C++