Calcolo numerico - riassunti, formulario

ii

tl

ii

I il

i i

;,i i i

ij i i ì

i I r i

i^j ì

I i i i i

t *t'

p**-& ÒL, i

(xnn',Ue*)}*"

* É*i.o*to

R(x;rfi.flG")} :'e(r'rrs

!-r1 ,

l'*ry , i

i:**${e"d,+"* iuh f i

[n

* rxl I

;'L,d-tt I

j

oJ+4.

*.uiq\ie.,*Al+' "do

<r-

i"to" =i

F-dS1,-{ i

;il"ic}frF,iEi{*;T*"

ipe-r i

i;ry,fr*,H I

i"'i i"l"l*-i-'*i i

yft."{srt" I

i i.i

t€ùlj*.$ i

ri\c,'*$cidb tr'nlt{"is

I Fx* i

1 i

i

r : i.{t';i-'*-.r1*la,i',*r$

r

i I i i i

I i I I i i i i I I I i I

i.A"f--3 .,8"

j;É'f,'

o.(d iq

dr€c.rd.,fuva?t.'.rl

frii l"

' I '

,;,***

ti*, '

6; ,'.'

:rtìro.l:,&t,r 'ir. -ur-l j

i

i I I

, ì

I, i I I

, i

i i

ioiiú'r;dl$ù r

tr-c*erà*e,-i

e-l* iLìrf

I,T-=t'"i:i iI, i

i,

iii1l.irIIiiii

;,it,iitiiiitiiiiiriij;i

iii:itì

ii:iiììiriìiijtìiriiii;t;iir;;i

ttiriaiiilllrljlrillrì

ilill:i, .' ',:,

'; i

i

I i ! i

. i ì | i

! i : I

I i | ,

i. i '-.r

I

i Ì

I I :. r

i I

i i

I i ; : I tii:iii:ii,

liiiìii.:ii"i;;li i i i

I i i

"ti:'i :ti i

r i i i I

i i iI I i

r

ìilliiiit: i ii

irj, : 'lii'l

ir, i

i

i ii ì i r, i

i i

ì : i r

I L I

i I I I

i

I I

i I i

liiì;iiliiiiiìii,r',iì,Ì!ii,:r i

r

lrljiitililiiììt':li,:iliit i

,9

PPotsLFt-Ift DI DÉL OI2DINE.

CSOGIU

-Ou {.:

a6pEl-rpú.c Ppogugna sD

Drc 0l.tf,

cl2-btrF

@c_E_ o

D6L QDtl,E

dE

ryx

ddo

$a tn po[0mo deJ onC.r'ne;

9s

qB(x)=o x>o

sg'og*

\g"cx)+ {

{uc"l =

to= Q

Lq' di

tti rirTct€t€o î=[ctr].

contlt€{o ìn

r.rchieilo

vierìe n

q

€o p.phùcra dgl

deto

Qqindi uc'ortO,t'cB 4g'orcf/ne.

d;

tn

drùeuy{rcì

Sa%erbtR-.0ncld.rYvra!:bc sr-s(gcnl eA

CDt 2-

t(\ @r't

). ocg0r{o uacQ,tr't€.l.

o0CJne/

d"

dQJ ۓitero

rl\ tohe qllsìAar/e.

co

rcfnoCu cnodo

dS,he, ,Ì\

vo-t q.BS+o :

e d.

t6C<) cr,Li gttG)

<.C$

U trC<)

g\CX)

E) = =

= ,

€ $

qiÀCIir,o ig sit.kmè t-^gltruj',€.':

q"€S-le

zr8)*3?G)+zuCX)=o X>o

U(ol=!=gLo)

:

€cè1=o Uq?ò

I

cart E 69 CX)

=U

ú-cx)

.tiorcl:ùncurd,cQo

cH.enqc,

Quindd. r

c

( x>o

SUe)-+ 2trOÒ

z'Cx)+ =

(rrck)

;eG)

\ =

ògt8ì=: 'sf?s'S*o(ffi'e3cx';)

at€Poi.

frPpuco

Discu-1itìa{*crre: :g

!=[o,{ tt 0,5 =F/n n

sùpecnq- fi-

L-o =Flfr B

p 9r

= =

= q+,...,,b

r.,,fi txi) qiît

i= l,

Xi = Xo €,0 OICX")

*.,

WW

C.l,Locoanne: xi€I

-3ECxi)-2uc<i)=o

a'Gi)=

I '

u'cxÈ)

r

{úffi, o

LZtóf = :

li'ir

;i',. C<f)S'.Zun).

.* 9

iefn r-fl

'el-4

(Hìl

l

:\-

ltì

1: p

l (llr&'

,t e

t.ctccno;

d,r".

,;.rll,,ii.i,

.r t & itl

oigKl

*.T ,o

r{gbn*Ù,'-*'

où1 Si

W

jbscr''èi,ru i

j;i

.qrl1n. i

i

i i i i'i

i--r-) i i i i i

i i ìi-i

r i i i

1e,=qlq,l".*Qol l

l

ìr,.i:l.ir';Ì .i iiil :i-,iii:

.n,o'n-t

r."f+re fiedthil{e

s'crs{o

o

.tn"li<ìtl

rr*ol

I i

'''

aefe,sttr] irirlirtilijlit:

ili io i, .'i,liiii.il

b l

dr\

drÀtgmgQ'arqs. hO*o"e

\/ri,i::,iliì,ìl ]::Ìil

']tttri . :

,ii ;]l

Ì\1,

c

i${r * H

*o* (re"L +

I t-{

tea,

'-el L*-(e"t;leu*\)+ffi*]

l) ;fi,ge*!f

fe,p+,

| Usìil

I'Rp,x;,fi I

ffig,(ú)l

à*, *

fi;lóÒj=-.#,;i gts,i i

.

l Ì :

i

I : i ir, x**{

i i-;

: i

i i I i i i- i i i I

ì ; f f I

i ,

I -l

'

L\t

.:iì l

,,'€1, ,

iIEirlfìi ìiXi',,,)l iiiIil

€ "

X,'-.ri

6ì,rL; :

iìIi;',1;iiiiii,ì;iiìi

iI l,!;irilli

ir,.tf:lli

, ; gbb"l*

ErtuÉ (J

ii

ì'

llii

1i1il1;:ìiri,,:.,11liìi,iii,l

:il,i]l

;i,iilii

ii:iiiil ' .

riÌiiril

:ii,:l

i I

I

7i ::

l cno*e,$e1f

Faftu

)

l

l *.

C*;o*s rw

; i i*

\**sl

i

I ^Jt

\til

iìl;,i l

ii

Algoritmo di Thomas

# tridiagonale

# soluzione lineare

sistema

#

# DELLA

HELP THOMAS

FLINZIONE

# INPIJT:

# gli

D principale

# vettore con elernenti della diagonale della matrice coefficienti

dei

gli

# coefîcienti

vettore con elementi sulla diagonale superiore della matrice dei

S gli

A inferiore

vettore con coefficienti

elementi sulla diagonale della matrice dei

# B termini noti

vettore dei

# OUTPUT:

# X vettore della

# soluzione

# Thomas(D,S,A,B):

def import librerie

#

from import

numpy zeros, size

# inizializzazione array

N: size(B)

V: zeros(N-l)

U: zerosQll)

alpha: zeros(N-1)

Y: zeros(N)

X: zerosQll) LU

# Fattortzzazione

V:S

:

u[0] D[0]

i

for in range(N-1):

: A[i]ru[i]

alpha[i] : *V[i]

U[i+ D[i+1

]

I ]-alphafi]

# soluzione del lineare

sistema

:

Y[0] B[0]

i

for in range(N- 1) :

: *Y[i]

Y[i+1 B[i+1

] ]-alphafi]

:

x[N-1] Y[N-1]/U[N-1]

i

for in range(N-2,-1,-1):

*xIi+

xtil:(Ytil-vlil til

])/u

1

# output

X

return ree

HgfcDt S\S6{t {rrEAU

tT€QSnVt

S ft

la g€ct-o

rn &0r\è""

SccnnQccre- crs1n&Cs ú:.,

a

b

L+D+U to nnatucadr"!f,C4.

doB.

R= t ó

b )

d.ra.pnce" (-o-

d{

m-u$erp- p*rtup"0, ry

fnryt"wo lrrnr++m:rE#hi+tr'ch'e v'ìe L"-

@-u;he e"ido

@

-rgffit

-lr,.u^.q"*""S"S#"dl e\-rftp

È

MM,n#fu..ru' @

Ch* A"-

Òcr..tfiQna- eo:*e

a*

tt H. 4.-

n^lo:./-ruftcqttrc

ìn-€dJ

S$nnncs\e pfJatt-

dC