• Lenguaje

    C

  • Descripción

    Encuentra la solución de f(x) = f(x)=x³+2x²+10x-20 con el método de la secante

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
#include <stdio.h>
#include <stdlib.h>
#include <math.h>

double f (double x)
{
    return x*x*x + 2.0*x*x + 10.0*x - 20;
}

const int ITERACIONES_MAXIMAS = 100;
const double TOLERANCIA = 0;

int main (void)
{
    double x0, x1, temp, err;
    int n = 1;
    printf ("M\202todo de la Secante para el c\240lculo de la funci\242n: f(x)=x\374+2x\375+10x-20");
    printf ("\n\nIngrese la aproximaci\242n inicial x0: ");
    scanf  ("%lf", &x0);
    printf ("Ingrese la aproximaci\242n inicial x1: ");
    scanf  ("%lf", &x1);
    printf ("\n%-20s%-20s%-20s\n", "n", "Xn", "Error");
    printf ("0                   %g\n", x0);
    do {
        err = fabs(x0-x1);
        printf ("%-20d%-20g%-20g\n", n, x1, err);
        if (err!=0)
        {
            temp = x1;
            x1 = x1 - (x1 - x0) * f(x1) / (f(x1) - f(x0));
            x0 = temp;
        }
        n++;
    } while (err>TOLERANCIA && n<=ITERACIONES_MAXIMAS);
    putchar ('\n');
    if (n<ITERACIONES_MAXIMAS)
        printf ("La soluci\242n es: %g", x1);
    else
        printf ("No se encontr\242 la ra\241z: cambiar aproximaciones iniciales o aumentar ITERACIONES_MAXIMAS");
    printf ("\n\n");
    system ("pause");
    return EXIT_SUCCESS;
}