OliForum Matematico

Salve ragazzi io dovevo svolgere questo programma in mathematica:

e = 1;
a = 2;
n = 3;
b = Table[2n, {i, 1, n}];
A = Table[1 + n*KroneckerDelta[i, j], {i, 1, n}, {j, 1, n}];
x0 = Table[2, {i, 1, n}];
r0 = b - A.x0;
x = x0;
r = r0;
i = 0;
While[Norm[r]>e*Norm ,
{
i++;
x = x + a*r;
r = b ...

errore ho ripetuto due volte la risposta

Per gli scalari quale si usa? e per i vettori?

Ok su table cè l'errore invece di 1 è 2n
b = Table[2n, {i, 1, n}];

Mi potresti fare un esempio...cosa devo correggere su questo programma per farlo funzionare correttamente non riesco a venirne a capo

e = 0.1;
α = 2;
n = 5;
b = Table[1, {i, 1, n}];
A = Table[1 + n*KroneckerDelta[i, j], {i, 1, n}, {j, 1, n}];
x0 = Table[1.5, {i, 1, n}];
r0 = b - A.x0;
x = x0;
r = r0;
i = 0;
While[Norm[r] ≥ ϵ*Norm ,
{
i++;
x = x + α*r;
r = b - A.x;
}
]
Print["Numero di iterazioni:"];
Print
Print ...

si xo sinceramente non è del tutto chiaro anche perchè la matrice A non era simmetrica..
Grazie Comunque

Salve ragazzi...dovrei svolgere Il sistema lineare Ax=b o meglio (P-N)x=b attraverso richardson....prendendo A simmetrica e P matrice identita.

Premetto conosco poco richardson e pochissimo mathematica...

Non so qualsiasi aiuto è buono...io parrtirei:

Quale la sintassi in matematica per fare una ...