[C] Znajdowanie msc. zerowego funkcji metodą bisekcji

Dla specjalistów Programowanie
#1

Witam, czy mógłby ktoś rzucić okiem czy dobrze napisałem ten program:

#include <stdio.h>
#include <stdlib.h>
 
double f(double x) {
return x*x*x + x -9;
}
double modul(double liczba) {
if(liczba<0) {
return -liczba;
}
else {
return liczba;
}
}
 
int main() {
float p,k,epsi,xs;
int i;
printf("podaj p\n");
scanf("%f", &p);
printf("podaj k\n");
scanf("%f", &k);
 
epsi=0.001;
 
if( (f(p)*f(k)) < 0 ) {
while(modul(p - k) > epsi) {
xs=(p+k)/2;
if( (f(p)*f(xs)) < 0 ) {
k=xs;
}
else if( (f(xs)*f(k)) < 0 ) {
p=xs;
}
i++;
}
printf("\nf(xs)= %f \nxs= %f\ni=%d", modul(f(xs)),xs,i);
}
else {
printf("Funkcja nie spelnia zalozen \n");
}
}

Z góry dziękuje!

#2

Z uwagi na brak wcięć (dzięki czemu kod jest nieczytelny), odpowiadam na pytanie: program jest źle napisany.

#3

Emm, nie możesz sobie wziąć jakiś 5 funkcji obliczyć dla nich miejsca zerowe z użyciem matematyki analitycznej i potem zobaczyć czy metoda numeryczna zwraca wyniki w przybliżeniu podobny?

#4

albo forum mi obcina wcięcia albo nie kopiują się z dev c++.

Kostek robiłem tak, i działa, chodziło mi bardziej czy czegoś nie brakuje, jakiejś dodatkowej instrukcji warunkowej, to że działa dla 5 nie znaczy że będzie działać dla wszystkich.

 

Pozdro…

#5

U mnie zostawia wcięcia. Nie pamiętam już jak sam pisałem coś takiego w szkole ale co w przypadku funkcji np. kwadratowych? Oba końce wybranego przedziału mogą być po jednej stronie osi ale przy pierwszym ifie program się zakończy i wyświetli że funkcja nie spełnia założeń. Chyba że tak ma być a ja po prostu nie pamiętam, że tak miało być. Jak dla mnie kod wygląda ok jednak po paru latach programowania zdecydowanie bardziej wolę pisać np. float poczatekPrzedzialu, koniecPrzedzialu, epsilon, polowaPrzedzialu niż tajemnicze pojedyncze literki, które nic nie mówią. Oczywiście na Twój użytek może być jednak jak będziesz musiał wrócić do bardziej skomplikowanego kodu i zobaczysz a,b,c,x,y,z,p,q,k to się przestraszysz :slight_smile:

#include <stdio.h>
#include <stdlib.h>
double f(double x)
{
    return x*x*x + x -9;
}
double modul(double liczba)
{
    if(liczba<0)
    {
        return -liczba;
    }
    else
    {
        return liczba;
    }
}
int main()
{
    float p,k,epsi,xs;
    int i;
    printf("podaj p\n");
    scanf("%f", &p);
    printf("podaj k\n");
    scanf("%f", &k);
    epsi=0.001;
    if( (f(p)*f(k)) < 0 )
    {
        while(modul(p - k) > epsi)
        {
            xs=(p+k)/2;
            if( (f(p)*f(xs)) < 0 )
            {
                k=xs;
            }
            else if( (f(xs)*f(k)) < 0 )
            {
                p=xs;
            }
            i++;
        }
        printf("\nf(xs)= %f \nxs= %f\ni=%d", modul(f(xs)),xs,i);
    }
    else
    {
        printf("Funkcja nie spelnia zalozen \n");
    }
}
#6

Marcin, tak właśnie zastanawiałem się nad parabolą, ale chyba na razie może tak zostać. Dzięki za sprawdzenie :slight_smile: