Tutti i Forum  INFORMATICA - Primo Anno  Analisi Matematica  Algoritmo di Halley

Nota: Devi essere registrato per poter inserire un messaggio.
Per registrarti, clicca qui. La Registrazione è semplice e gratuita!

Larghezza finestra: 640 x 480 800 x 600 1024 x 768 1280 x 1024 Nome Utente: Password: Modo: Base  Aiuto  Finestra  Formato: Carattere Andale Mono Arial Arial Black Book Antiqua Century Gothic Comic Sans MS Courier New Georgia Impact Lucida Console Script MT Bold Stencil Tahoma Times New Roman Trebuchet MS Verdana   Grandezza 1 2 3 4 5 6   Colore Nero Rosso Giallo Rosa Verde Arancione Viola Blu Beige Bordò Verde Oliva Blu Marino Marrone Verde Limone Icona Messaggio:                             Messaggio:   * Il codice HTML è OFF * Il Codice Forum è ON Smilies [quote][i]Messaggio inserito da Michael Gardel[/i] [br]salve tutti ! qualcuno mi può aiutare a implementare l'algoritmo del metodo di halley per il calcolo dello zero d'un polinomio in Matlab ? ci ho provato ma mi vengono sempre degli errori. mi dà errore anche la funzione feval. grazie %%% %%function [x,X] = halley(f,df,d2f,x,tol,Maxit) %%% %% %INPUT % df(x) derivata prima della function, % ddf(x) derivata seconda della function, % x punto iniziale, % e toleranza, % n max iterazione. %% OUTPUT: % x soluzione dell'equazione % X vecteur approssimazioni. %%%%% %% function [x,X] = halley(f,df,d2f,x,tol,Maxit) X=x; fx =feval(f); for k=1:Maxit if abs(fx) < tol break; end; dfx = feval(df); d2fx = feval(d2f); d=dfx^2-reshape(d2fx*fx,size(dfx))/2; h = d(dfx*fx); x = x - h; X = [X,x]; fx=feval(f); if norm(h,inf) + norm(fx,inf) < tol break; end; end; end [/quote]    Allega file   Clicca qui per inserire la tua firma nel messaggio. Clicca qui per sottoscrivere questa Discussione.