亚洲欧美第一页_禁久久精品乱码_粉嫩av一区二区三区免费野_久草精品视频

? 歡迎來到蟲蟲下載站! | ?? 資源下載 ?? 資源專輯 ?? 關于我們
? 蟲蟲下載站

?? test.cpp

?? 牛頓法解多項式的根 包含重根的判斷和多項式的輸入
?? CPP
字號:
#include <stdio.h>

#include <malloc.h>





double Polynomial_Root(double c[], int n, double a, double b, double eps);

 

int main()
{

    double eps = 0.00005;

    double *c;

    int n;

    double a, b;

    int i;

 

    while (scanf("%d", &n)!= EOF){

       c = (double *)calloc((n+1), sizeof(double));

       for (i=n; i>=0; i--) 

           scanf("%lf", &c[i]);

       scanf("%lf %lf", &a, &b);

       printf("%.4lf\n", Polynomial_Root(c, n, a, b, eps));

       free(c);
    }

 
    return 0;

}



/*
double ff(double c[], int n, double x)
{
       double sum = 0;
	   int i;
			
	   for ( i = n; i >= 0; i-- )
			sum = sum * x + c[i];

	   return sum;
}

double dff(double c[], int n, double x)
{
     double sum = 0;
	 int i;

     for ( i = n; i >= 0; i-- )
		 sum = sum * x + c[i] * i;
	     
	 return sum;
}
*/

#include <math.h>

void HornerMethod(int n, double a[], double x0, double *f, double *df)
{
	double y, z;
	int j;

	y = a[n];
	z = a[n];

    for(j = n-1; j >= 1; j--)
	{
		y = x0 * y + a[j];	
		z = x0 * z + y;
	}

	y = x0 *y + a[0];
	
	*f = y;
	*df = z;
}

/*5 1 -5 10 -10 5 -1
-100 20
二分法結果時是1.0001,牛頓法就是1.0000*/


int NewtonMethod(double c[], int n, double eps, double p0, double *p, double a, double b)
{
   double f, df;
   int i;
   int MAX_N = 1000;


   for(i = 0; i < MAX_N; i++)
   {
	   //f = ff(c, n, p0);
	   //df = dff(c, n, p0);
	   HornerMethod(n, c, p0, &f, &df);
	   
	   if((df==0)&&(f==0))
	   {
           if(p0 < a || p0 > b)
      		   return -1;

	       if(fabs(*p) <= (eps*1000))
		   {
			   *p = +0.0;
			   return 0;
		   }
			
		   *p = p0;
		   return 0;
	   }
	   else
		   if((df==0)&&(f!=0))
		      return -1;

       *p = p0 - (f) / (df);
	   
	   if(fabs(*p - p0) <= eps)
	   {
		   if(p0 < a || p0 > b)
      		   return -1;

	       if(fabs(*p) <= (eps*1000))
		   {
			   *p = +0.0;
			   return 0;
		   }
			
		   *p = p0;
		   return 0;
	   }

	   p0 = *p;
   }

   return -1;
}

double Polynomial_Root(double c[], int n, double a, double b, double eps)
{
     double p0, p, h;
	 int MAX_N = 20;
	 int i = 0, m, mark;
	 
     if(a > b)
	 {
	    p = a;
		a = b;
		b = p;
	 }

     m = 20;
	 h = (b - a) / (double)m;
	 for(i = 0; i <= m; i++)
	 {
		p0 = a + h * (double)i;
		mark = NewtonMethod(c, n, eps/1000, p0, &p, a, b);

		if(mark == 0)
		{
			  return p;
		}
		else 
			  continue;
	 }

	 printf("Newton's method failed!!\n");
	 return 0; 
}

?? 快捷鍵說明

復制代碼 Ctrl + C
搜索代碼 Ctrl + F
全屏模式 F11
切換主題 Ctrl + Shift + D
顯示快捷鍵 ?
增大字號 Ctrl + =
減小字號 Ctrl + -
亚洲欧美第一页_禁久久精品乱码_粉嫩av一区二区三区免费野_久草精品视频
亚洲视频在线观看三级| 91原创在线视频| 97久久超碰国产精品电影| 欧美精三区欧美精三区| 日本一二三四高清不卡| 日韩精品欧美精品| 91色九色蝌蚪| 国产日韩欧美一区二区三区综合| 亚洲va天堂va国产va久| www.色综合.com| 久久久久久电影| 日本不卡中文字幕| 欧美图片一区二区三区| 自拍偷拍亚洲综合| 国产 日韩 欧美大片| 精品国产露脸精彩对白| 亚洲一区二区三区四区在线免费观看 | 亚洲风情在线资源站| 国产精品538一区二区在线| 欧美精品第一页| 伊人性伊人情综合网| av中文字幕一区| 国产欧美日韩亚州综合| 九九精品视频在线看| 欧美一区二区三区视频在线| 亚洲一区二区三区视频在线播放| thepron国产精品| 国产精品国产三级国产a| 国产成人av一区二区| 久久综合999| 国内精品嫩模私拍在线| 精品日韩av一区二区| 青青草精品视频| 欧美一区国产二区| 日本伊人午夜精品| 日韩午夜激情av| 蜜桃视频一区二区三区| 日韩视频在线观看一区二区| 蜜臀av性久久久久蜜臀aⅴ四虎| 欧美精品tushy高清| 亚洲bt欧美bt精品777| 欧美日韩国产小视频| 日韩激情视频在线观看| 91精品国产入口| 精品一区二区av| 国产亚洲精久久久久久| av激情亚洲男人天堂| 亚洲黄色在线视频| 91麻豆精品国产91久久久资源速度| 爽好多水快深点欧美视频| 欧美一区二区精品| 久草在线在线精品观看| 中文字幕av在线一区二区三区| av成人免费在线| 亚洲va欧美va人人爽| 日韩欧美一级在线播放| 丁香亚洲综合激情啪啪综合| 亚洲女同ⅹxx女同tv| 欧美日韩国产一级二级| 国产一区二区福利| 亚洲男人的天堂在线观看| 欧美日本国产视频| 国产成人8x视频一区二区| 国产精品美女久久久久久| 在线亚洲+欧美+日本专区| 蜜桃av一区二区| 国产精品入口麻豆九色| 欧美美女bb生活片| 国产精品一区二区久久不卡| 亚洲丝袜制服诱惑| 日韩欧美一级二级| av电影天堂一区二区在线观看| 日韩经典一区二区| 日韩理论片一区二区| 日韩一级二级三级精品视频| 99久久婷婷国产综合精品| 免费成人在线观看| 亚洲一区在线看| 国产精品久久久久久久久免费桃花| 在线观看不卡一区| 国v精品久久久网| 人妖欧美一区二区| 亚洲精品中文字幕在线观看| 日韩精品一区二区三区四区视频 | 91麻豆精品国产| 成人免费看片app下载| 男男视频亚洲欧美| 亚洲黄色性网站| 国产女人18毛片水真多成人如厕| 欧美精品亚洲一区二区在线播放| 成人a免费在线看| 国产九色sp调教91| 蜜臀久久99精品久久久久久9| 亚洲精品午夜久久久| 国产精品网站在线播放| 精品国产免费久久| 日韩一区二区在线看| 在线观看免费视频综合| 91片黄在线观看| 成人三级伦理片| 国产成人在线免费| 国产在线精品视频| 精品在线免费观看| 男人的j进女人的j一区| 亚洲成av人片在线| 夜夜精品视频一区二区| 亚洲人成人一区二区在线观看| 久久精品人人做人人综合 | 正在播放亚洲一区| 在线观看欧美精品| 欧美色男人天堂| 欧美午夜影院一区| 欧美视频一区二区三区| 欧美日韩一区久久| 欧美一区二区三区不卡| 欧美一级理论性理论a| 欧美一级理论片| 精品国产乱码久久| 久久久国产一区二区三区四区小说 | 国产乱码精品一品二品| 国产河南妇女毛片精品久久久| 国产一区在线看| 国产98色在线|日韩| 成人一道本在线| 99国内精品久久| 欧美性猛交一区二区三区精品| 精品视频免费在线| 日韩一区二区中文字幕| 精品动漫一区二区三区在线观看| 精品国产凹凸成av人网站| 日本一区二区三区在线观看| 国产精品福利电影一区二区三区四区| 国产精品水嫩水嫩| 亚洲精品欧美综合四区| 天天综合色天天| 精品亚洲国内自在自线福利| 国产一区二区视频在线播放| 国产91精品露脸国语对白| 一本色道久久综合亚洲91| 欧美日韩免费观看一区二区三区| 91精品国产综合久久福利| 精品日韩一区二区| 最新国产精品久久精品| 午夜视频久久久久久| 精彩视频一区二区三区| av一二三不卡影片| 欧美精品一二三| 欧美—级在线免费片| 一区二区三区精品久久久| 日本视频免费一区| 国产999精品久久| 欧美高清视频一二三区| 国产欧美精品区一区二区三区| 亚洲精品成人天堂一二三| 麻豆91免费看| 色网站国产精品| 精品噜噜噜噜久久久久久久久试看| 国产精品丝袜在线| 日韩精品电影在线观看| 成人av网站在线观看免费| 欧美日韩精品系列| 亚洲国产精品传媒在线观看| 亚洲bt欧美bt精品| 成人在线综合网| 欧美一区二区观看视频| 综合婷婷亚洲小说| 国产一区二区成人久久免费影院| 欧美伊人久久久久久午夜久久久久| 欧美成人一区二区三区| 亚洲国产中文字幕| jizz一区二区| 日韩精品一区二区三区老鸭窝| 一区二区三区鲁丝不卡| 国产精品一区二区三区乱码| 欧美一级理论性理论a| 亚洲精品国产a久久久久久| 国产精品1区二区.| 欧美成人一区二区三区| 天天色综合成人网| 色综合欧美在线| 国产精品久久久久久一区二区三区 | 亚洲一区免费观看| 成人av网址在线观看| 国产午夜精品一区二区三区四区| 日本不卡一区二区三区高清视频| 97国产精品videossex| 中文字幕欧美激情一区| 国精品**一区二区三区在线蜜桃| 欧美电影一区二区| 亚洲国产欧美一区二区三区丁香婷| eeuss鲁片一区二区三区在线看| 久久久久国产精品厨房| 国产综合色产在线精品| 日韩精品一区二区三区视频播放 | 欧美午夜精品久久久久久超碰| 欧美激情一区三区| 国产在线精品一区二区夜色| 91精品免费在线| 天天av天天翘天天综合网| 欧美性生活影院|