Αριθμητικη Ανάλυση I (Εργαστήριο)

[fatness]

Ευχαριστώ πολύ kostas!

[coyote]

Η συναρτηση nr ειναι μια συναρτηση της βιβλιοθηκης του matlab ή μας δινεται κάπου?
Το matlab 7.3 που έχω δεν ξερει τι είναι το nr.m

[O kanenas]

Η συνάρτηση nr δεν βρίσκεται στην βιβλιοθήκη του matlab. Την φτιάχνεις εσύ και ο κώδικας της συνάρτησης σου δίνεται έτοιμος από το φυλλάδιο του 2ου εργαστηρίου ή μπορείς να βρεις έτοιμο το αρχείο σε ένα ftp του οποίου την διεύθυνση δεν θυμάμαι.

Ρώτησα τον Φιλλιπάκη τι εννοεί με αυτό το h και μου είπε ότι το δίνει για να χρησιμοποιήσουμε τον σχετικό τύπο μεταβολής. Αυτός ο τύπος είναι αυτός που ανέφερε ο dimk? Μη βγει λάθος ολόκληρη η εργασία...

[fatness]

Last favour..
Μπορεί κάποιος να ποστάρει τον κώδικα της Newton-Raphson όπως δόθηκε στα εργαστήρια; (Το nr.m δλδ).. Θανξ

[O kanenas]

Αυτός είναι ο κώδικας του nr.m

Κώδικας: Επιλογή όλων

function [x, xlist, iter] = nr (f, df, x0, tol, maxit)
%Root of f(x)=0 by Newton-Raphson's method.
%Method terminates when relative change > tol
%or maximum number exceeds of iterations.
if nargin < 3, error('Insufficient input \n'); end;
if nargin < 4, tol=eps; end;
if nargin < 5, maxit=50; end;
x = x0;
iter = 0;
xdiff = inf;
xlist = x;
while xdiff >= tol
    iter = iter + 1;
    xold = x;
    x = x-feval(f,x)/feval(df,x);
    xdiff = abs(x-xold)/abs(x);
    xlist = [xlist;x];
    if iter >= maxit, error('Not converged after maxit iterations.');
    end;
end;

[rainel]

Help me!
Συμφωνα με τον τυπο : lim ( (f(x+h) - f(x)) / h ) = f'(x)
Μπορω να υπολογισω τη 2η παραγωγο??????????
Αν ναι , ποιος ειναι???????????????

[nikos]

lim((f'(x+h)-f'(x))/h)=f''(x)

[sunday]

Στην εργασία στο γ ερώτημα αν βάλω χ0 ισο με 5 δεν μου βγάζει αποτέλεσμα μου βγάζει NaN.Πρέπει να αλλάξω το χ0?γιατι δοκίμασα μερικά αλλά και πάλι δεν βγαίνει...τι αλλαγή πρέπει να γίνει?

[nikos]

Αν χρησιμοποιήσες το σωστό nr, τότε έχεις ορίσει λάθος την f ή την df!!!

[sunday]

ναι πράγματι είχα κάνει λάθος το df.Σ'ευχαριστώ πολύ αν δεν μου το έλεγες δεν θα το έβλεπα Όμως εξακολουθώ να μην καταλαβαίνω τι εννοεί στην εκφώνηση όταν λέει "λαμβάνοντας υπόψη την γραφική παράσταση κάντε τις κατάλληλες αλλαγές"

[nikos]

Πρέπει να κάνεις τη γραφική παράσταση της f σε ένα διάστημα (π.χ. [0,5]), έτσι ώστε να δεις γραφικά που κοντά είναι η ρίζα της f. Έπειτα βάλε x0 έναν αριθμό κοντά στη ρίζα και θα δεις οτι η μέθοδος συγκλίνει πιο γρήγορα!!! Ουσιαστικά αυτό ζητάει λέγοντας για απαραίτητες αλλαγές.

[sunday]

ok σ'ευχαριστώ πολύ!!!!!

[ian]

paidia mhpws mporei kapoios na mou ekshghsei ti tropopoihseis prepei na kanw sto nr gia to erwthma D? Oti kai an exw dokimasei den trexei.

[O kanenas]

Θυμάται κανείς την διεύθυνση του ftp που σε μερικά εργαστήρια πάμε για να κατεβάσουμε διάφορα m-files, όπως το gauss.m και το solvelu.m?

[Ulm]

Η σελίδα είναι:
ftp://147.102.205.20