我們利用Newton Method (綠色部分的程式) 即可求得方程式的解:

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

#define F(x) (double) (pow(x,2)-2*x-15)
#define Del_X (double) 0.00000001

int main(void)
{
   int i;
   double dF,x;
   char buf[12],*c;

   printf("\nInput Start X value : ");
   c=gets(buf);
   x=atof(c);

   for(i=0;i<100;i++)
{
     dF=(F(x+Del_X)-F(x))/Del_X;
     x=x-(F(x)/dF);
   }

   printf("\nRoot : %f",x);

   return 0;
}

對於不同的函數只要更改粉紅色部分定義函數 F(x) 的程式碼即可