c程序设计 8.12 用牛顿迭代法求根。方程为:ax^3+bx^2+cx+d=0 ,系数a,b,c,d由主函数输入。求X在1附近的一个实根。求出后由主函数输出.

//https://baike.baidu.com/item/%E7%89%9B%E9%A1%BF%E8%BF%AD%E4%BB%A3%E6%B3%95/10887580?fr=aladdin#4
//百度牛顿迭代法
    #include <stdio.h>
    #include <math.h>

    double solut(double a,double b,double c,double d)
    {
        double x1=1,x,f,f1;

        //迭代
        do
        {   x=x1;
            f=((a*x+b)*x+c)*x+d;//原函数
            f1=(3*a*x+2*b)*x+c;//求导之后的函数
               x1=x-f/f1;
        } while(fabs(x1-x)>=0.00000001);
        return x1;
    }

    int main()
    {
        double solut(double ,double ,double ,double );
        double a,b,c,d;
        scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
        printf("%.2f\n",solut(a,b,c,d));
        return 0;
    }

结果:

时间: 2024-10-06 02:36:27

c程序设计 8.12 用牛顿迭代法求根。方程为:ax^3+bx^2+cx+d=0 ,系数a,b,c,d由主函数输入。求X在1附近的一个实根。求出后由主函数输出.的相关文章

【c语言】为下面的函数原型编写函数定义,这个字符串参数必须包含一个或者多个数字,函数应该把这些数字字符转换为整数并返回这个整数。

/*为下面的函数原型编写函数定义: int ascii_to_integer(char *str); 这个字符串参数必须包含一个或者多个数字,函数应该把这些数字字符转换为整数并返回这个整数. 如果字符串参数包含了任何非数字字符,函数就返回零.请不必担心算数溢出. 提示:这个技巧很简单:你每发现一个数字,把当前值乘以10,并把这个值和新的数字所代表的值相加*/ #include <stdio.h> int ascii_to_integer(char const *str) { int sum =

利用牛顿迭代法求平方根

求n的平方根,先如果一推測值X0 = 1,然后依据下面公式求出X1,再将X1代入公式右边,继续求出X2…通过有效次迭代后就可以求出n的平方根,Xk+1 先让我们来验证下这个巧妙的方法准确性,来算下2的平方根 (Computed by Mathomatic) 1-> x_new = ( x_old + y/x_old )/2 y (x_old + -----) x_old #1: x_new = --------------- 2 1-> calculate x_old 1 Enter y: 2

【leetcode】【二分 | 牛顿迭代法】69_Sqrt(x)

题目链接:传送门 题目描述: 求Sqrt(x),返回整数值即可. [代码]: 1 #include<bits/stdc++.h> 2 using namespace std; 3 const int N = 1e6+10 ; 4 /* 5 int mySqrt ( int x ){ 6 int L = 1 , R = N , mid , ans = 0 ; 7 while ( L <= R ){ 8 mid = ( L + R ) >> 1 ; 9 if( mid <=

蓝桥杯练习系统 矩阵翻硬币 大数,牛顿迭代法 难度:2

http://lx.lanqiao.org/problem.page?gpid=T126 明显,对于一个格子(i,j),设f(i)为i的约数个数,则(i,j)的翻转次数为(f(i)-1)*(f(j)-1)+1, 而只有翻转次数为奇数,也就是f(i),f(j)都为奇数的格子开始才是反面, 又因为f(i)为奇数当且仅当i为完全平方数,所以只需统计n,m中各有多少个完全平方数,然后相乘即可, 也就是sqrt(n)*sqrt(m), 但是因为n,m是大数,必须要用大数方法解决,这里采用了java的Big

经典算法:牛顿迭代法求平方根

//牛顿迭代法求平方根 1 double mysqrt(double num) 2 { 3 double x = num/2; 4 double y = 0; 5 do{ 6 x = x/2+num/(2*x); 7 y = x*x-num; 8 if(y<0) y = -y; 9 } while(y>0.0001); 10 return x; 11 } 12 int main(int argc, char* argv[]) 13 { 14 printf("%.3f",my

OJ刷题之《牛顿迭代法求根》

题目描述 用牛顿迭代法求根.方程为ax3+bx2+cx+d=0.系数a,b,c,d的值一次为1,2,3,4,由主函数输入.求x在1附近的一个实根.求出根后由主函数输出.结果保留两位小数. 输入 系数a,b,c,d的值 输出 x在1附近的一个实根 样例输入 1 2 3 4 样例输出 -1.65 提示 主函数已给定如下,提交时不需要包含下述主函数 /* C代码 */ int main() { double solut(double ,double ,double ,double ); double

牛顿迭代法 求方程根

牛顿迭代法 牛顿迭代法(Newton's method)又称为牛顿-拉夫逊方法(Newton-Raphson method),它是牛顿在17世纪提出的一种在实数域和复数域上近似求解方程的方法.多数方程不存在求根公式,因此求精确根非常困难,甚至不可能,从而寻找方程的近似根就显得特别重要. 方法使用函数f(x)的泰勒级数的前面几项来寻找方程f(x) = 0的根.牛顿迭代法是求方程根的重要方法之一,其最大优点是在方程f(x) = 0的单根附近具有平方收敛,而且该法还可以用来求方程的重根.复根,此时线性

牛顿迭代法应用——求数的平方根和立方根

牛顿迭代法,从一个值开始,用无限逼近的方式得出结果. #include<stdio.h> #include<math.h> int main() { double a; double x; scanf("%lf",&a);//求a的平方根和立方根 x = a/2; ////平方根///// while( fabs(x*x-a) > (1e-6) ) { x = (x+a/x)/2; } printf("%lf\n",x); //

牛顿迭代法求Logistic回归

接着上次的一篇文章:http://blog.csdn.net/acdreamers/article/details/27365941 在上次这篇文章中,对于Logistic回归问题,我们已经写出它的最大似然函数,现在来求最大似然估计.所以对似 然函数求偏导数,得到了个方程,即 由于我们只要根据这个方程解出所有的即可,但是这不是一件容易的事,还有Logistic回归求的是最大似 然估计,我们在多元函数求极值问题中也说过,导数等于零的点可能是极大值,极小值或者非极值.所以还要靠一个 叫Hessian