在R中,solve
将求解线性方程组,ode
可以求解初值问题微分方程。我有一个方程组,我无法弄清楚如何适应任何一个。有人能告诉我它是如何解决的吗?
其中有三点已知,但第四点有点棘手。我有最大y值的x值,但我不知道最大y值是多少。例如,假设在(0,0),(1,1)和(4,0)处有已知点。最大y值在x = 3处,因此第四个点是(3,ymax)。
有人向我展示了如何设置满足这些方程的方程组:
实际上有可能因为你需要x = 3的y值应该是最大的。因此,4次多项式有5个系数待确定,你有以下等式:
y(0)= 0
y(1)= 1
y(4)= 0
dy / dx(3)= 0(x = 3时的一阶导数应为0)
d2y / dx2(3)< 0(x = 3时的二阶导数应为负)
我一直试图弄清楚如何使用R解决这个问题,但我正在打砖墙。我不能使用solve
,因为方程式4和5是差分的。但我无法弄清楚如何使用ode
因为初始值是几个不同的x值。我错过了一些明显的东西吗?