Metoda secantelor.

Fie dată ecuaţia f(x) = 0

Pentru utilizarea metodei secantelor vom cere îndeplinirea aceloraşi condiţii asupra funcţiei f(x) ca şi pentru metoda precedentă
Ca şi la metoda precedentă în calitate de soluţie vom căuta un punct interior al intervalului, dar de această dată el nu va fi mijlocul segmentului, ci punctul intersecţiei dreptei ce trece prin punctele (a, f(a)) şi (b,f(b)) cu axa 0X.

Ideea generală a metodei este următoarea: prin punctele (a, f(a)) şi (b,f(b)) se duce o dreaptă. Se determină punctul c în care ea intersectează axa 0X. Apoi se determină semnul f(c). Extremitatea în care semnul funcţiei coincide cu semnul f(c), trece în c. Procesul se repetă, până nu obţinem o apropiere suficientă de soluţia exactă. Această interpretare a metodei o face să fie tot atât de lentă ca şi metoda bisecţiei.

Dacă însă suplimentar, vom cere existenţa şi semnul constant al derivatei de ordin 2 pentru funcţia f(x) pe tot segmentul [a,b]. ( pe segmentul dat funcţia să fie sau concavă sau convexă.), putem aplica o variaţie mai puternica a acestei metode:

Ecuaţia dreptei ce trece prin 2 puncte date este determinată de relaţia:


În cazul intersecţiei cu axa X y=0, deci putem calcula poziţia punctului c pe axa X:

Pentru a realiza calculul soluţiei vom crea şirul {xn} conform regulei:

La paşii următori valorile elementelor şirului se vor forma, ţinând cont de următoarele situaţii:

a) dacă funcţia este convexă (f’’(x) > 0 ) (secanta este plasată mai sus de graficul funcţiei) atunci extremitatea segmentului iniţial, în care valoarea funcţiei este pozitivă rămâne nemişcată (fixă), iar cealaltă extremitate trece de fiecare dată în punctul nou calculat.

b) dacă funcţia este concavă (f’’(x) < 0 ) (secanta este plasată mai jos de graficul funcţiei) atunci extremitatea segmentului iniţial, în care valoarea funcţiei este negativă rămâne nemişcată (fixă), iar cealaltă extremitate trece de fiecare dată în punctul nou calculat.

Cu alte cuvinte, metoda secantelor are un punct fix, celălalt se modifică de fiecare dată. Observăm, că fixă este extremitatea, pentru care are loc relaţia:
f(x) f’’(x) > 0.
Pentru a calcula cealaltă extremitate a noului segment, trebuie să ţinem cont de următoarele cazuri posibile:

1) f’(x) > 0 f’’(x) > 0 2) f’(x) > 0 f’’(x) < 0 3) f’(x) < 0 f’’(x) > 0 4) f’(x) < 0 f’’(x) < 0

Vom cerceta cazul 1 – funcţia este convexă, crescătoare. (Celelalte cazuri se cercetează la fel)
Deoarece a doua derivată e mai mare ca 0, graficul funcţiei este inferior coardei. Pornind din extremitatea în care funcţia este pozitivă vom obţine intersecţia cu axa 0x într-un punct interior al [a,b] (mai întâi se intersectează axa, apoi se ajunge la extremitatea opusă.) Repetând procedura vom obţine un şir de puncte interioare, care tind spre soluţia exactă. Dacă însă vom fixa extremitatea negativă, coarda generată la pasul 2 va intersecta mai întâi graficul, apoi axa 0X, ceea ce poate genera părăsirea segmentului [a,b].
Deoarece funcţia este crescătoare, vom fixa extremitatea B, luând în calitate de x0 extremitatea A.
Conform formulei pentru secantă

după care a  x1
Repetând procedura, obţinem:

Vom demonstra acum, că şirul de aproximări a= x0 , x1 , …, xn … converge către soluţia ecuaţiei f(x) = 0
Din formula recurentă

rezultă că xn+1 > xn  n (numitorul fracţiei este pozitiv deoarece funcţia este crescătoare, numărătorul – negativ, deoarece xn, este punctul de intersecţie a coardei cu 0X, iar graficul funcţiei în acest punct este inferior coardei, deci + negativ, (b - xn) – pozitiv, deci tot produsul este negativ, cu – în faţă, trece în +).
Aşa dar avem şirul x0 < x1 < … Şirul este crescător şi mărginit, deci există limita lui 1. Trecând la limită în formula recurentă, obţinem

sau, după înlocuire

de unde rezultă 0 = f(1). (numitorul fracţiei e diferit de 0, al doilea factor – la fel). Deci 1 este soluţia ecuaţiei iniţiale. 1 = , ceea ce şi trebuia de demonstrat.
Celelalte trei cazuri se demonstrează analogic.

Algoritmizarea metodei

1. Determinăm, care va fi extremitatea fixă a segmentului. Pentru aceasta calculăm punctul de intersecţie (x1) a coardei ce trece prin f(a)) şi (b,f(b)) cu axa 0X. Drept extremitate fixă vom lua acel punct pentru care semnul funcţiei diferă de semnul f(x1).
2. Extremitatea “flotantă” primeşte valoarea calculată x1.
3. Calculăm următoarea aproximaţie conform formulei:

4. Repetăm pasul 3 până nu obţinem soluţia cu exactitatea cerută.
Estimarea erorii.

Faptul că şirul aproximărilor succesive prin metoda coardelor converge către soluţia exactă implică următoarea concluzie: cu cât mai multe iteraţii ale metodei vor fi realizate, cu atât mai bine va fi aproximată soluţia. Totuşi, am putea determina o formulă, care permite estimarea exactă a erorii de calcul, şi, prin urmare, exactitatea soluţiei obţinute.

Vom considera, că prima derivată a funcţiei f(x) este continuă pe segmentul [a,b].
Fie atunci m1 şi M1 două numere pozitive, pentru care are loc relaţia

0 < m1  | f’(x) |  M1 < + 

Pentru simplitate, vom cerceta cazul, când aproximările succesive se fac după formula (2) (cazul, când e folosită formula (1) se demonstrează analogic). Aşa dar,


(exprimăm xn prin xn-1)

(Ţinând cont de faptul că 0 = f(), îl adăugăm la numărătorul fracţiei, apoi estimăm diferenţa între soluţia exactă şi cea aproximativă.


Pentru a deduce formula finală vom folosi următoarea teoremă (se demonstrează în cursul de analiză matematică)

Teorema Lagrange Fie f:[a,b] R, continuă şi derivabilă pe [a,b] . Atunci  c(a, b) astfel încât:

f(b) –f(a) = (b -a) f’(c)

 fără demonstraţie 

Vom aplica formula din teoremă aparte pentru partea stângă şi partea dreaptă a ultimei egalităţi.

Înlocuind, obţinem :

Deci dacă e necesar să obţinem soluţia cu exactitate dată , vom repeta aproximările până la obţinerea inegalităţii

Aici xn şi xn+1 sunt două aproximaţii succesive iar M1 şi m1 corespunzător marginea superioară şi inferioară a f’(x) pe [a, b].