Obliczanie pierwiastka funkcji nieliniowej metodą siecznych

Question

Witam mam napisać program do liczenia jak w temacie. Sugerowałem się tym opisem tej metody. Wzór i wykres mojej funkcji jest taki:

Napisałem coś takiego ale program się w ogóle nie kończy jako że to dopiero początek mojej przygody z c++ więc pewnie popełniłem jakis elementarny błąd.

#include <stdio.h>
#include <math.h>
using namespace std;

double f(double x)
{ return (1+(3*(pow(x, 2.)))); }

int main() {
    double x=-2., y, xk[99];
    do
    { for(double x=-2.; x<=3.; x+=0.1)
    {
       for(int i=0; i<=98; i++)
        {
        xk[i-1]= f(x);
        xk[i+1]=xk[i]-(f(xk[i]*(xk[i]-xk[i-1])))/(f(xk[i]-f(xk[i-1])));
        printf("x=%f \n",xk[i-1]);
        }
    }
    }while (y>=0.00001);
   printf("Wyniki: x=%f \n",x );
   return 0;
}

Z góry dzięki za pomoc :D

Answer 1

Wziąłeś pochodną do szukania pierwiastka... Ona nie ma miejsce zerowych :D

Comments

hmm po wpisaniu tej wyżej jest to samo.

#include <stdio.h>
#include <math.h>
using namespace std;

double f(double x)
{ return (pow(x,3.)+x-3); }

int main() {
    double x=-2., y, xk[99];
    do
    { for(double x=-2.; x<=3.; x+=0.1)
    {
       for(int i=0; i<=98; i++)
        {
        xk[i-1]= f(x);
        xk[i+1]=xk[i]-(f(xk[i]*(xk[i]-xk[i-1])))/(f(xk[i]-f(xk[i-1])));
        printf("x=%f \n",xk[i-1]);
        }
    }
    }while (y>=0.00001);
   printf("Wyniki: x=%f \n",x );
   return 0;
}

 

---

Masz trochę przekombinowane z pętlami i tablicami. Na ekranie wypisuje się 100 razy to samo. Usunąłem pętlę for i zobacz ile mam wyników bez powtórzeń:
 

#include <stdio.h>
#include <math.h>
#include <windows.h>
using namespace std;

double f(double x)
{ return (pow(x,3.)+x-3); }

int main() {
    int i=0;
    double x=-2., y, xk[99];
    do
    { for(double x=1.; x<=3.; x+=0.1)
    {

        xk[i-1]= f(x);
        xk[i+1]=xk[i]-(f(xk[i]*(xk[i]-xk[i-1])))/(f(xk[i]-f(xk[i-1])));
        printf("x=%f \t %f \n",xk[i-1], f(xk[i-1]));


    }
    }while (y>=0.00001);
   printf("Wyniki: x=%f \n",x );
   return 0;
}

Ile wynosi y? Jest zadeklarowany, ale nic z nim nie robisz. Dlaczego nagle x jest z przedziału od -1 do 24? 

---

rzeczywiście y został mi z poprzedniej wersji właśnie x powinien być z przedziału -2 do 3 jak na wykresie. A masz jakiś pomysł jak zrobić to żeby nie było przekombinowane z tablicami?

---

Np. tak: http://eff10.internetdsl.tpnet.pl/programowanie/mz_fun/pages/regulafalsi.htm

Najlepiej jakbyś kierował początkiem, a nie gotowym kodem.