[delphi] problem : Bisekcijska metoda

//bisekcijska metoda
      Hmin:=R4; // je 0.05
      Hmax:=100;
      C := (Hmin + Hmax) / 2
      yA := jakost(Hmin)
      yC := jakost(C );
    while (abs(yC)>1E-10) and (Hmax-Hmin > 1E-10) do
    begin
      if yA * yC < 0 then begin
            Hmax = C;
      else begin
            Hmin := C;
            yA := yC;
      end;
      H := (Hmin+Hmax)/2;
      yC := jakost(C );
    end;
    LabelRezultat4.Caption:=('h = ' + FloatToStr(H));  

Hvala za pomoč, samo me baga kaj si mislu s tem "jakost(C)" pa to. Ne štekam.

jakost ( C ) je procedura ja

edit: sem popravil

1.) Kam je zginla enačba za Eh (A si to zamenju s temi yC, yA, ...)?
Predvidevam da s to enačbo izračunah jakost v točki x (oziroma H). Naredi si torej proceduro jakost z ustreznimi parametri, ki ti vrne izračunano jakost.
2.) Kaj sploh predstavlja yC, yA, kako to deklerirat (?), a kot funkcijo ali klasična spremenljivka?
To sta spremenljivki, ki predstavljata jakost v točkah Hmin in C. Deklariraj jih kot realno število.

Pa pojdimo čez algoritem:
Imaš:
- funkcijo, ki ti izračunba jakost v določeni točki (H oziroma rečimo ji kar x)
Se pravi: y = jakost(x)
- točki a in b, ki sta izbrani točki na nekem intervalu (x os), kjer imata vrednosti funkcije jakost (vrednosti y v točkah a in b) nasprotna predznak, oz. a je negativna, b je pozivitna ali obratno.
Se pravi če imata a in b enak predznak bisekcija ne bo delovala, ker med točkama ni ničle.
- točka c, ki leži med točkama a in b

a := leva točka
b := desna točka
c := (a+b)/2
yA := jakost(A)
yC := jakost(C )
Se pravi tukaj imamo dve točki (x, y), ki ležita na tvoji funkciji.

dokler (yC ni dovolj blizu ničle) in dokler (razlika med trenutnima krajiščema a in b ni dovolj majhna) {
če je ena od vrednosti yA in yC desno krajišče b premakni v točko c
drugače pa levo krajišče a premakni v c in zračunaj jakost v tej točki

izračunaj novo točko c med krajiščema a in b (eno od teh dveh se je premaknilo zato je potrebno spet vzeti vmesno točko)
izračunaj jakost v točki c
}

Kokr jst vem, ne morš kr napisat "jakost(C )", ker je delphi čist navdušen z errorji , dekleracija bla bla.
Proceduro jakost moraš deklarirati.

Upam da je dovolj jasno drugače pa vprašaj.

Eh malo slabo sva začela. Ti iščeš točko (x,y), kjer tvoja funkcija (jakost) seka x.
Se pravi rešuješ enačbo jakost(x) = 0.
To točko ti najde tvoj algoritem, toda moraš imeti eno točko kjer je jakost(x) manjša od nič in drugo kjer je jakost(x) večja od nič.

Nisem elektrotehnik, tako da bova težko debatirala o jakosti sami, toda tole je čista numeričnam matematika.

Hmin:=R4; // je 0.05
Hmax:=100;
C:=(Hmin + Hmax)/2;
yA:=(U4*sqrt(sqr((Hmin)/(R4))-1)/(R4*((Hmin)/(R4)-1)*ln((Hmin)/(R4)+sqrt(sqr((Hmin)/(R4))-1))));
yC:=(U4*sqrt(sqr(©/(R4))-1)/(R4*(©/(R4)-1)*ln(+sqrt(sqr(©/(R4))-1))));
while (abs(yC)>1E-10) and (Hmax-Hmin>1E-10) do
begin
if yA*yC < 0 then
Hmax:=C
else
begin
Hmin:=C;
yA:=yC;
end;
H:=(Hmin+Hmax)/2;
yC:=(U4*sqrt(sqr(©/(R4))-1)/(R4*(©/(R4)-1)*ln(+sqrt(sqr(©/(R4))-1))));
end;
LabelRezultat4.Caption:=('h = ' + FloatToStr(H));


A je to to, kot si mislu? Kaj pa un tazadnji yC (a je enkak unmu zgoraj, al maš notri končni H?)

Mislim da je namesto besede function procedure.

Nope, pascal (ali delphi) loči funkcije ki nič ne vračajo - procedure in funkcije ki nekaj vrnejo - function.

Še nekaj sem spregledal
v čisto prvem postu napišeš

if (Eh>E4) then   //če je Eh>500 potem

kar se sklada s sliko, k si jo dal ampak se ne sklada z mojim razumevanjem štroma. Namreč - ali se ne bi morala električna poljska jakost zmanjševati od vodnika proti neskončnosti?
Ker če je to res, potem samo obrni tisti 'večje' v manjše' in bi moralo delat. Torej tako:

if (Eh<E4) then   //če je Eh>500 potem

Sem obrnu un znak in mislim, da dela.
Ubistvu je bilo cel cajt sam to narobe.

Bom še mal stestiru, poročam.

EDIT: nekje dela, nekje ne :) Morm se mal poigrat.