Noi di ing. meccanica non abbiamo l'esercitazione sul calcolo degli autovalori, quindi, se vuoi, posso darti come confronto solo i valori che ho ricavato per la seconda esercitazione sul gradiente coniugato modificato.
Al momento vista l'ora carico solo i grafici di convergenza delle varie matrici, così possiamo discuterne e confrontarne un pò l'andamento...
Nota: tranne la ricerca dell'autovalore minimo e massimo, i grafici relativi al GCM sono presi per risolvere un problema con termine noto unitario (per capirsi `vec b=[1,1,1,....,1]^T` e con `vec x_0 = [0,0,...,0]^T`
Domani o in settimana dovrò anche riscrivere meglio il codice del CGM (devo staccare la parte delle Correzioni Residue dal GCM stesso)
Intanto:
Matrice 9x9 data a Calcolo Numerico (vedi sopra)


Vettore Soluzione:
`vec x = ((0.6875),(0.875),(0.6875),(0.875),(1.125),(0.875),(0.6875),(0.875),(0.6875))`
`lambda_n=1.17157288`
`vec v_n=((0.842954356),(1.19211748),(0.842954356),(1.19211748),(1.68590871),(1.19211748),(0.842954356),(1.19211748),(0.842954356))`
`lambda_1=6.82842712`
`vec v_1=((0.250000004),(-0.353553391),(0.250000004),(-0.353553391),(0.499999991),(-0.353553391),(0.250000004),(-0.353553391),(0.250000004))`
Matrice 48x48



Vettore soluzione: `vec x= ((0.000335667076),(3.59115768E-006),(8.77033532E-005),(3.34230027E-008),(-1.43883264E-007),(-6.38659204E-008),(0.00015592896),(2.91595248E-006),(8.77000031E-005),(3.32134783E-008),(-1.39621638E-007),(-6.05909231E-008),(0.000155944476),(-1.24297723E-006),(0.000104768961),(-2.08309665E-08),(-1.4531924E-007),(-6.37615957E-008),(0.000335657214),(1.76477397E-006),(0.000105739382),(1.35581645E-008),(-1.42420792E-007),(-6.76969885E-008),(0.000197404787),(2.1522149E-007),(6.15461883E-005),(1.00698963E-007),(-1.86127324E-007),(-1.54127365E-006),(9.54325524E-005),(2.01871354E-006),(6.1537878E-005),(1.00920125E-007),(-1.32721228E-007),(-1.27320848E-007),(9.5419131E-005),(-9.53751198E-007),(0.000101852929),(1.51484636E-007),(-1.02786275E-007),(-1.27045852E-007),(0.000197279249),(1.58274496E-006),(0.00010253161),(8.53880094E-008),(-2.95566785E-007),(-1.51046779E-006))`
`lambda_1=3.01517938E+009`
`vec v_1=((-1.20719091E-006),(-9.66460599E-007),(9.95595236E-007),(9.99304607E-007),(-0.000274951335),(-0.000277210959),(0.000545749477),(0.000324505015),(8.12435731E-006),(1.1074333E-006),(-0.00131340903),(0.0876396804),(0.000561536827),(-0.000324837593),(-8.01190621E-006),(-0.000108340114),(-0.0075078818),(0.0893120163),(-2.93801177E-006),(1.05774856E-006),(-8.50594008E-007),(-3.11969888E-005),(-0.000816881994),(-0.000246060547),(3.37376926E-008),(-2.23276421E-008),(-1.03503462E-009),(8.14633052E-006),(-1.29816732E-005),(-0.000708960702),(3.64357864E-005),(0.00323495301),(7.54599449E-007),(8.31574135E-006),(-0.000131244244),(0.701466936),(3.68313265E-005),(-0.00323515702),(-8.19581077E-007),(-7.46684525E-005),(-0.00052389614),(0.701565493),(-1.372766E-005),(4.67861497E-007),(1.25178956E-008),(-0.000187789844),(-0.00332661919),(-7.11694421E-006))`
Matrice 90x90

Matrice 100x100

Matrice 132x132

Matrice 1138x1138


PS: alcune matrici non pubblicate le ho trovate sempre nel sito delle dispense/esercitazioni di Gambolati. Basta cercare con Google:
ES: matrici 1138 site:dmsa.unipd.it
Appena posso aggiungo commenti / note e soluzioni del sistema. Dato che devo leggermente editare il programma che avevo creato spero di non giungere con la nuova versione a risultati diversi

Vi farò sapere...