Quer encontrar 0 em funções reais? Aí vai o primeiro algoritmo, lembrando que o algoritmo da bisseção de acordo com os matemáticos é o menos eficaz para se achar 0 de funções devido ao número de iterações que ele executa.
OBS: este código ainda não pega sua função, você tem que inserí-la manualmente no código. Neste código exemplo calculo a função F(x) = x^3 - 9*x + 3.
#include <stdio.h>
#include <stdlib.h>
int main(){
float a = 0., b = 0., E = 0., M = 0., Fa = 0., Fb = 0., Fm = 0.;
int k = 0;
printf( "Digite valor de A:" );
scanf( "%f", &a );
printf( "Digite valor de B:" );
scanf( "%f", &b );
printf( "Digite valor de E:" );
scanf( "%f", &E );
printf("\n");
do{
M = ( a + b ) / 2.;
Fa = ( a*a*a ) - (9. * a) + 3.; //insira a sua função aqui para calcular F(a)
Fm = ( M*M*M ) - (9. * M) + 3.; //insira a sua função aqui para calcular F(M)
Fb = ( b*b*b ) - (9. * b) + 3.; //insira a sua função aqui para calcular F(b)
if ( Fa * Fm < 0 )
b = M;
else
a = M;
k += 1;
printf( "iteracao = %d\n", k );
printf( "A = %f\nB = %f\n", a, b );
} while( b - a >= E || Fa >= E || Fb >= E);
//se for verdade ele continua o laço
printf("\nb - a = %f\n", b - a);
printf( "A solucao final eh: %f\n", ( a + b ) / 2. );
return 0;
}
[1] Comentário enviado por HelioCampos em 02/05/2008 - 07:56h
Dependendo da função que a pessoa colocar, o valor de float estoura antes de vc conseguir achar o ponto pois vc colocou seu E (Epsilon) igual a 0. O algoritmo serve apenas para funções MUITO bem comportadas ( sem muitos graus de liberdade ) e, mesmo assim, pode não dar certo.
[2] Comentário enviado por rafaelhenrique em 02/05/2008 - 11:20h
Nesta função novamente HelioCampos, eu não defini E constantemente como 0, mas sim inicializei E como 0 pois, como diz Deitel, quando se inicializa uma variável quando se declara, o programa se torna mais rápido e mais livre de erros de stacks overflow, o E na verdade vai ser digitado pelo usuário, portanto E não será 0, será o valor que o usuário digitar depois que E for declarado.
[3] Comentário enviado por g_s.lima em 04/05/2010 - 18:13h
Eu sugiro que você utilize um epsilon cte de acordo com o tipo de dado que você utilizar. No caso do float, temos precisão simples. De uma estudada na biblioteca float.c.
[4] Comentário enviado por manolosk em 09/09/2013 - 16:55h
Galera vou ressuscitar o topico.. Rafael vc poderia me ajudar mudando o código? è bem simples.. faço engenharia mas não manjo de programação, e um dos trabalhos que tenho que é fazer um programa para calcular a bisseção.
A função que vc usou foi F(x) = x^3 - 9*x + 3.
A minha esta na imagem anexada...
E também o programa deve perguntar quantas iterações será realizada.
Voce pode me ajudar? Vai quebrar um galhão!
Te mandie msg no face mas vc nao viu heheh.