高分求一段简单的mathematica程序,计算一个高阶多项式的解
多项式形式为ax^10+bx^9+cx^8+dx^7+ex^6+fx^5+gx^4+hx^3+ix^2+jx^1+k其中a,b,c,d,f,等为可输入的参数,对这个方程求...
多项式形式为ax^10+bx^9+cx^8+dx^7+ex^6+fx^5+gx^4+hx^3+ix^2+jx^1+k 其中a,b,c,d,f,等为可输入的参数,对这个方程求二阶导数,并让这个二阶导数等于0,计算x的值,多谢 。毕业论文计算数据要用。导师说这是个很简单的问题,分数不是问题。
展开
1个回答
展开全部
Mathematica代码如下:
co
={a,
b,
c,
d,
e,
f,
g,
h,
i,
j,
k};
expr
=
co.x^Range[10,
0,
-1]
NSolve[D[expr,
{x,
2}]
==
0,
x]
因为方程次数高于5,所以一般只有数值解。
第一句是给系数赋初值,对应x次数由高到低(一定要先赋值,否则没法计算)。
第二句是产生表达式,这里用了矢量的标量积来缩减表达式。
第三句是求导,求数值解。
作为例子,比如:
co
=
Range[11];
expr
=
co.x^Range[10,
0,
-1]
NSolve[D[expr,
{x,
2}]
==
0,
x]
co
={a,
b,
c,
d,
e,
f,
g,
h,
i,
j,
k};
expr
=
co.x^Range[10,
0,
-1]
NSolve[D[expr,
{x,
2}]
==
0,
x]
因为方程次数高于5,所以一般只有数值解。
第一句是给系数赋初值,对应x次数由高到低(一定要先赋值,否则没法计算)。
第二句是产生表达式,这里用了矢量的标量积来缩减表达式。
第三句是求导,求数值解。
作为例子,比如:
co
=
Range[11];
expr
=
co.x^Range[10,
0,
-1]
NSolve[D[expr,
{x,
2}]
==
0,
x]
已赞过
已踩过<
评论
收起
你对这个回答的评价是?
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询