PARI / GP中某个范围内的函数零

时间:2019-07-11 08:20:41

标签: equation-solving pari pari-gp

我正在尝试查找[0,1]范围内的函数零。我的函数通常看起来像这样:

F(t)=2*sin(8*Pi*t)+2*sin(2*Pi*t)-1

到目前为止,我尝试过的幼稚方法是将函数的粗略图绘制到Pari plot(t=0,1,F(t))中,然后从中“猜测”较小的范围[a1,a2],以表示第一个零将是。然后,我一直在使用solve(t=a1,a2,F(t))来找到那个零。然后找到第二个小范围[a2,a3],然后重复进行,直到找到所有零为止。

显然不是一个很好的方法,但是它已经完成了工作,到目前为止,我看过的许多示例只有两个零。但是,我现在正在查看的示例(例如上面定义的F(t))具有更多的零,因此我正在寻找一种更快的方法。

我知道solvestep,但是我尝试使用它,但是它从不返回任何解决方案。我猜这种类型的函数实际上不适合它使用的拆分方法(或者可能我使用的是错误的方法。)。

有没有什么好方法可以使我在返回的[0,1]中得到该函数的零,而不必通过这种漫长而幼稚的方法?

0 个答案:

没有答案