| Autore |
 Discussione  |
|
|
kanc
Utente medio
Regione: Puglia
Prov.: Bari
|
 Inserito il - 29/10/2008 : 13:24:32
|
function [y]=laplace (A)
[m,n]=size (A)
if m~=n
error ('La matrice non è quadrata')
end
if n==1
y=A
end
else
y=0
for j=1:n
A1j=[2:n,(1:j-1,j+1,n)]
y=(y+[(-1)^(j+1)]*A(1,j)*laplace (A1j))
end
endfunction
perchè non funziona???
l'ho copiata pari pari da iavernaro...
grazie in anticipo
|
Modificato da - kanc in Data 29/10/2008 13:25:40
|
|
|
Rio
Utente medio
|
Inserito il - 29/10/2008 : 13:57:12
|
Su quale riga dà errore?
|
www.festamaggiore.it
Festa Maggiore - Scoprire e diffondere le tradizioni di Terlizzi |
 |
|
|
Dr House
Utente medio
Città: Bari
|
Inserito il - 29/10/2008 : 14:07:41
|
evidentemente non l'hai copiata bene
Citazione:
Messaggio inserito da kanc
function [y]=laplace (A)
[m,n]=size (A)
if m~=n
error ('La matrice non è quadrata')
end
if n==1
y=A
end <--- NON SERVE
else
y=0
for j=1:n
A1j=[2:n,(1:j-1,j+1,n)] <--- CI SONO PARECCHIE IMPRECISIONI
y=(y+[(-1)^(j+1)]*A(1,j)*laplace (A1j))
end
endfunction
Se ci ragioni l'end dopo l if non serve, sopratutto se ci metti il ramo else.
Poi dove ho scritto ci sono parecchie imprecisioni stai utilizzando l'operatore : per fare una sottomatrice ma non specifichi da che matrice te la deve fare.
Poi nella formula dove ce la chiamata ricorsiva hai usato parentesi quadre dove dovevi usare parentesi tonde..
Spero di esserti stato utile, ciao!
|
"La realtà è fatta di cose vere e cose supposte.
Le cose vere le mettiamo da parte, ma le supposte dove le mettiamo?" |
|
|
| |
Discussione |
|
|
|
dove ho sbagliato???
