Υπολογισμός ρίζας σε C με την μέθοδο Newton-Raphson (Λύθηκε)

[giannis_mathiou]

Καλησπέρα παιδιά, χρειάζομαι την βοήθειά σας. Πρέπει να κάνουμε ένα πρόγραμμα σε C που να υπολογίζει την τετραγωνική ρίζα ενός αριθμού με την μέθοδο Newton-Raphson.
Μας είπε ότι γίνεται Χj+1=1/2(Xj+a/Xj) μέχρι να είναι δύο διαδοχικά Χ ίσα. Δηλαδή για παράδειγμα η ρίζα του 4 είναι Xo=1, X1=1/2(1+1/4)=2,5 , X2=1/2(2,5+4/2,5)=2,05 ,X3=2, X4=4
Αυτό που έχω κάνει εγώ είναι το εξής:

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

#include <stdio.h>

float sqNR(float);

int main(){
float num;


printf ("\nΠληκτρολογήστε τον αριθμό για να βρώ την τετραγωνική ρίζα.\n");
scanf ("%f",&num);

if (num>0)
printf ("\nΗ τετραγωνική ρίζα με την μέθοδο Newton-Raphson είναι ίση με: %.1f",sqNR(num));

else
printf ("\nΔεν υπολογίζεται η τετραγωνική ρίζα αρνητικού αριθμού.\n");

return 0;
}

float sqNR(float N){
float x=1,xn=1/2*(x+(N/x));
while (x!=xn){
x=xn;
xn=1/2*(x+(N/x));
}

return xn;
};

Περνάει από το compile κανονικά και μόλις δώσω τον αριθμό φαίνεται σαν να έχω σταματήσει να τρέχει και δεν κάνει τίποτα. Μπορείτε να μου πείτε τι κάνω λάθος; Και πως πρέπει να το κάνω;

Σας ευχαριστώ πολύ!

[Dimitris]

Βάλε μερικά printf σε διάφορα σημεία του προγράμματος να δεις που διακλαδίζεται, τι τιμές έχουν οι μεταβλητές.

[giannis_mathiou]

Δηλαδή τι ακριβώς εννοείς; Που να βάλω; Μέσα στην συνάρτηση δεν μου τα εμφανίζει.

[Dimitris]

Αφου δεν εμφανίζει ένα printf στην αρχή της sqNR τι μπορεί να σημαίνει αυτό; Αν βάλεις ένα printf μέσα στο while τι θα δείχνει;

edit: Βασικά διάβασε ξανά τη βασική θεωρία αριθμητικών μεθόδων. Δηλαδή πότε συγκλίνει μια αριθμητική μέθοδος; Πότε σταματάει; Ποιά κριτήρια τα κριτηρία τερματισμού; Μπορεί να ικανοποιηθεί ένα κριτήριο τερματισμού; Μήπως να προσθέσουμε κι ένα δεύτερο για να είμαστε σίγουροι ότι θα τερματίσει ο αλγόριθμος;

[giannis_mathiou]

Το ότι δεν καλείται σωστά η συνάρτηση ή ότι δεν παίρνει τιμές από την main;

Edit: Θα πρέπει να βάλω αν είναι έναν μετρητή για τα loops την συνάρτησης να τερματίζεται στα 100 ας πούμε;

[Dimitris]

Για να το κάνω πιο απλό τι σημαίνει ισότητα στους floating point numbers;

[giannis_mathiou]

Δύσκολο, σχεδόν αδύνατο να είναι ίσοι γιατί υπάρχουν πολλά δεκαδικά; Άρα πρέπει να περιορίσω τα δεκαδικά ή να του θεωρήσω ακεραίους;

[simosx]

Το ότι δεν καλείται σωστά η συνάρτηση ή ότι δεν παίρνει τιμές από την main;

Edit: Θα πρέπει να βάλω αν είναι έναν μετρητή για τα loops την συνάρτησης να τερματίζεται στα 100 ας πούμε;

Είσαι στο σημείο όπου πρέπει, για κάθε απορία που μπορείς να έχεις για το πρόγραμμα, να βρίσκεις τρόπο να καταλαβαίνεις τι συμβαίνει.
Έτσι, για τη τιμή που δίνεις κατά την εκτέλεση του προγράμματος, μπορείς να τυπώνεις τη μεταβλητή μετά από κάθε εντολή μόνο και μόνο για να καταλάβεις πως μεταβάλεται.

Είναι σημαντικό να μπεις σε αυτήν τη διαδικασία όπου είσαι σε θέση να καταλάβεις πλήρως το πρόγραμμα, αντί να περιμένεις να «παίξει» από μόνο του. Πρέπει να είσαι σε θέση να κάνεις αλλαγές/προσθήκες, να δοκιμάζεις τι αποτέλεσμα έχουν, και πάλι από την αρχή. Αν δε μπεις στη διαδικασία αυτή (όπου για τώρα μπορεί να σου πάρει αρκετό χρόνο), τότε θα δυσκολευτείς πάρα πολύ στο μέλλον.

[simosx]

Δύσκολο, σχεδόν αδύνατο να είναι ίσοι γιατί υπάρχουν πολλά δεκαδικά; Άρα πρέπει να περιορίσω τα δεκαδικά ή να του θεωρήσω ακεραίους;

Δοκίμασε και το ένα, και το άλλο.

[giannis_mathiou]

Λοιπόν δοκίμασα με ακεραίους και τρέχει αλλά δεν βγάζει σωστό αποτέλεσμα. Δηλαδή συνέχεια 0 ή -nan. Όπότε λείπει κάτι άλλο πρίν.

Edit: Με αυτή την αλλαγή

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

float sqNR(float N){
float x=1,xn=0;
int i=0;

while (x!=xn){
xn=x+(N/x)/2;
x=xn;
i++;
printf ("\t%.1f\t%.1f %10d\n\n",xn,x,i);
}

return xn;
};

μου βγάζει το εξής:
Πληκτρολογήστε τον αριθμό για να βρώ την τετραγωνική ρίζα.
4
3.0 3.0 1


Η τετραγωνική ρίζα με την μέθοδο Newton-Raphson είναι: 3.0