Metodo di doppia falsa posizione in Fibonacci

Da testwiki.
Vai alla navigazione Vai alla ricerca

Template:F Il metodo di doppia falsa posizione di origine indiana, noto anche come metodo elchataym, consente di affrontare problemi riconducibili ad equazioni lineari della forma ax=b o della forma ax+b=c. In questo articolo vediamo come questo metodo viene esposto nel Liber abbaci di Leonardo Fibonacci; come l'autore consideriamo solo casi con a,b,c>0).

A differenza del metodo di falsa posizione, o regula falsi nell'elchataym sono scelte arbitrariamente due 'false posizioni' da cui ricavare due approssimazioni distinte della condizione fissata che la soluzione esatta deve soddisfare.

Per maggiore leggibilità, come in genere si usa fare, utilizziamo il simbolismo algebrico; indichiamo con fp1, fp2, a1, a2, s, v, rispettivamente le due false posizioni, le due approssimazioni, la soluzione cercata e il valore noto della condizione che la soluzione soddisfa.

Supponiamo, senza ledere in generalità, che fp1<fp2 e quindi a1<a2 (essendo, nell'equazione ax+b=c, a>0 per i casi considerati).

Volendo riassumere in formule, l'elchataym è basato sulla seguente proporzione:

(a2a1):(fp2fp1)=(va2):(sfp2)
(1.1)

da cui segue, per la regola del quarto proporzionale, che

s=fp2+(fp2fp1)(va2)(a2a1)
(1.2)

La (1.1) può essere ricavata direttamente dalla regula falsi.
Infatti, poiché per ogni falsa posizione si ha a1:v=fp1:s,a2:v=fp2:s, per le proprietà delle proporzioni, si ottiene che

a2:fp2=a1:fp1
(a2a1):(fp2fp1)=a2:fp2=v:s

e similmente

v:s=(va2):(sfp2) ;

da qui, per la proprietà transitiva di rapporti uguali, si ha la (1.1).

Mostriamo come il metodo elchataym possa essere applicato alla risoluzione di equazioni del tipo ax=b, riconsiderando come esempio il problema dell'albero, cioè l'equazione (7/12)h=21.

Se scegliamo fp1=h=12 come prima falsa posizione, la prima approssimazione risulta essere a1=(7/12)12=7, da cui segue il primo errore e1=ba1=217=14 . Dalla seconda falsa posizione fp2=h=24 segue invece a2=(7/12)24=14, e il secondo errore e2=ba2=2114=7 .

Osserviamo che, aumentando di 12 la falsa posizione, nel passare da fp1 a fp2, per linearità dell'equazione, l'errore diminuisce per un fattore 7 da e1 a e2.

Ci chiediamo di quanto occorre aumentare fp2 per ottenere la soluzione e ridurre così e2 a zero. In altre parole, vogliamo determinare (sfp2) che soddisfi

7:12=7:(sfp2)
(1.3)

Basta applicare la regola del quarto proporzionale per ricavare

s=fp2+(sfp2)=24+1277=36

Si noti che la (1.3) è la proporzione (1.2) applicata all'esempio considerato, in quanto a2a1=(va1)(va2)=e1e2 .

Fissate arbitrariamente due false posizioni, indichiamo con e1, e2 gli errori che si ottengono come differenza fra il valore noto v e le rispettive approssimazioni (vale a dire e1=va1, e2=va2); la relazione (1.2) può essere riscritta come

s=fp2+e2(fp2fp1)(e1e2)
(1.4)

oppure come

s=(fp2e1fp1e2)(e1e2)
(1.5)

Quest'ultima espressione rappresenta in formula una variazione del metodo elchataym, ovvero quello che Fibonacci presenta come metodo di ‘aumento e diminuzione'. Osserviamo che, a seconda della scelta di fp1 e fp2, e1 ed e2 possono essere errori per eccesso o per difetto, cioè, algebricamente e1, e2 possono essere rispettivamente quantità negative o positive.

Questo non crea alcun problema dal punto di vista algebrico. Le formule (1.4), (1.5) hanno validità generale se consideriamo e1, e2 con i loro rispettivi segni. È importante tuttavia sottolineare che Fibonacci nel Liber abbaci lavora unicamente con numeri positivi (anche nei primi capitoli, introducendo l'operazione di sottrazione, tratta solo differenze positive tra un numero ed un secondo minore); gli errori di approssimazione sono considerati solo come quantità positive e rappresentano, a seconda della scelta delle false posizioni, l'eccesso o il difetto dell'approssimazione rispetto al valore noto.

Per questo motivo, nella trattazione dell'elchataym, fissate le false posizioni fp1, fp2, Fibonacci distingue tre casi:

  1. entrambi gli errori sono per difetto (in simboli e1,e2>0; quando fp1,fp2<s)
  2. entrambi gli errori sono per eccesso (in simboli e1,e2<0; quando fp1,fp2>s)
  3. un errore è per difetto e l'altro per eccesso (in simboli e1>0,e2<0; quando fp1<s<fp2)

Tale distinzione assume maggior rilevanza quando viene introdotto il metodo di 'aumento e diminuzione'.

Mostriamo come quest'ultimo possa essere dedotto dall'elchataym, lavorando, come proposto da Fibonacci, con dei segmenti.

Distinguendo i tre casi precedenti, sia .ab. (indichiamo i segmenti con la notazione usata nella traduzione in lingua inglese del Liber abbaci curata da Laurence E. Sigler) la soluzione esatta di un qualunque problema risolvibile con il metodo di doppia falsa posizione.

Caso 1. fp1,fp2<s

Siano .ag. e .ad. la prima e la seconda falsa posizione, e siano .ez. e .iz. come in figura (1.1) gli errori d'approssimazione corrispondenti.

Poiché .ei.=e1e2 e .gd.=fp2fp1, applicando la proporzione (1.1),

.ei.:.gd.=.iz.:.db.
(1.6)

da cui segue .db..ei.=.gd..iz.. Allora:

.ez..ad.=(.ei.+.iz.).ad. =.ei..ad.+.iz..ad. =.ei..ad.+.iz.(.ag.+.gd.) =.ei..ad.+.iz..ag.+.db..ei. =.ei.(.ad.+.db.)+.iz..ag. =.ei..ab.+.iz..ag.

cioè

.ab.=.ez..ad..iz..ag..ei.=e1fp2e2fp1e1e2
(1.7)

Caso 2. fp1,fp2>s

Sia .af. la prima falsa posizione e .ae. sia la seconda. Siano .gi. e .gk. il primo e il secondo errore di approssimazione. Poiché è possibile applicare l'elchataym, deve valere che

.ki.:.ef.=.gk.:.be.
(1.8)

da cui si ha .ef..gk.=.ki..be.

Allora:

.gi..ae.=(.gk.+.ki.).ae. =.gk..ae.+.ki.(.ab.+.be.) =.gk..ae.+.ki..ab.+.ef..gk. =.gk.(.ae.+.ef.)+.ki..ab. =.gk..af.+.ki..ab.

cioè

.ab.=.gi..ae..gk..af..ki.=e1fp2e2fp1e1e2
(1.9)

Caso 3. fp1<s<fp2

Siano .ag., .ad., .ez., .zi., rispettivamente la prima e la seconda falsa posizione e il primo e il secondo errore di approssimazione, come in figura. Osserviamo che in questo caso, per il metodo "elchataym", deve valere la seguente proporzione

.ei.:.gd.=.zi.:.bd.
(1.10)

ove .ei. è la somma dei due errori (il primo per difetto, il secondo per eccesso); analogamente, se si considera (sfp1),

.ei.:.gd.=.ez.:.gb.

da cui si ricava che .ez..gd.=.gb..ei..

Allora

.ez..ad.+.zi..ag.=.ez..ag.+.ez..gd.+.zi..ag. =.ez..gd.+(.ez.+.zi.).ag. =.ei..gb.+.ei..ag. =.ei.(.gb.+.ag.) =.ei..ab.

cioè vale effettivamente

.ab.=.ez..ad.+.zi..ag..ei.=e1fp2+e2fp1e1+e2
(1.11)

Come abbiamo visto, anche il metodo elchataym, come la regula falsi, è basato solo su un argomento di proporzione (vale a dire la (1.1)).

L'elchataym può inoltre essere utilizzato per la risoluzione di problemi riconducibili ad equazioni lineari della forma

ax+b=c
(1.12)

come procedimento alternativo alla risoluzione algebrica diretta dell'equazione.

Se infatti consideriamo due false posizioni x1, x2 e le sostituiamo nella (1.12), ottenendo così le due approssimazioni c1=ax1+b, c2=ax2+b, osserviamo che c2c1=a(x2x1), cioè a=(c2c1)/(x2x1) e b=c2x2(c2c1)/(x2x1). Riportando nella (1.12) i valori trovati per a e b, segue che

x=x2+(cc2)(x2x1c2c1)
(1.13)

del tutto analoga alla (1.2).

Rappresentazione grafica del metodo elchataym

L'applicazione del metodo elchataym per la risoluzione dell'eq. (1.12) può essere rappresentata graficamente con gli strumenti della geometria analitica.

Tracciamo il grafico della retta y=ax+b (prendiamo in esame solo il caso a,b>0 e consideriamo per comodità solo i punti di ascissa x positiva); vogliamo determinare il valore esatto x=s tale che as+b=c, con c fissato.

Scelte arbitrariamente le due false posizioni x1, x2 e calcolate le rispettive approssimazioni, ovvero le corrispondenti ordinate c1=y(x1), c2=y(x2), distinguiamo i tre casi:

  1. se x1<x2<s
  2. se s<x1<x2
  3. se x1<s<x2

Consideriamo il caso 1.

Come si può vedere dalla figura, poiché il segmento BE è parallelo al segmento AF e BC è parallelo a DF, i triangoli ABC e BDE hanno angoli congruenti e quindi sono simili. Questo implica che

BC:AC=DE:BE

da cui segue

s=x2+BE=x2+DEACBC

che coincide con la (1.2), in quanto

DE=cc2va2AC=x2x1fp2fp1BC=c2c1a2a1

Nel caso 2. consideriamo invece la similitudine tra i triangoli BDE e ADF.

Poiché i lati corrispondenti sono in proporzione,

DE:BE=DF:AF

da cui si ricava

s=x2AF=x2BEDFDE,

che coincide con la (1.2), poiché

BE=x2x1fp2fp1DF=c2ca2v=(va2)DE=c2c1a2a1

La proporzione

DF:AF=DE:BE

giustifica anche il caso 3. in cui una falsa posizione è maggiore e l'altra minore della soluzione esatta s.

Dalla figura si può ricavare

s=x2BE=x2AFDEDF

che coincide con l'espressione (1.2), in quanto

AF=x2x1fp2fp1DE=c2ca2v=(va2)DF=c2c1a2a1