弦截法求解方程

/*弦截法求解方程  Micro_lee*/

#include <stdio.h>
#include <string.h>
#include <math.h>

void asterisk_triangle(int n);

float f(float x)
{
 return x*( x*(x-5)+16 )-80;
}

float xpoint(float x1, float x2)
{
 return (x1*f(x2) - x2*f(x1))/( f(x2)-f(x1) );
}

float root(float x1, float x2)
{
 float x,y,y1,y2;

 y1=f(x1);
 y2=f(x2);
 do{
  x=xpoint(x1,x2);
  y=f(x);

  if(y*y1 > 0)
  {
   y1=y;
   x1=x;
  }
  else
  {
   y2=y;
   x2=x;
  }
 
 }while( fabs(y)>=0.0001);

 return (x);

}

int main()
{
 float x1,x2,f1,f2,x;

 do{
  printf("input x1,x2:/n");
  scanf("%f,%f", &x1, &x2 );
  f1=f(x1);
  f2=f(x2);
 }while(f1*f2>=0 );

 x=root(x1,x2);

 printf("A root is %8.4f/n/n", x);

 return 0;

}

 

 

 

 

 

 

void asterisk_triangle(int n)
{
 int i,j,k, num=0;

 for(i=1; i<=n;  i++)
 {
  
        num= (2*i-1) > n ? 2*(n+1-i)-1: 2*i-1 ;

  for(k=1; k<=(n-num)/2; k++)
   printf(" ");
  for(j=1; j<=num; j++)
   printf("*");
  
     printf("/n");
 }

}