double a = 1, Lspan = 30, Lcable = 33, fn, fdn, dfn, j;
fn = (2 * a * (Math.Sinh(Lspan / 2 * a))) - Lcable;
fdn = (2 * (Math.Sinh(Lspan / 2 * a)) - ((Lspan / 2 * a) * Math.Cosh(Lspan / 2 * a)));
dfn = -fn / fdn;
do
j = a + dfn;
while (dfn > 0.00000000001);
这就是我的代码。我知道它写得很糟糕!我是初学者。所以这是我的问题,从上面我想做的是牛顿拉夫森方法。
x n+1 = x n - f(x n )/f'(x n )
现在我知道如何计算 f(x n )/f'(x n ),如上所示,但我遇到的问题是 x n部分。
对于那些不知道的人,用 newton raphson 给出 x n的猜测估计值,假设为 1 并填写它。然后您得到的新答案将成为方程的 x n并填写一个等等,直到你达到了极限。现在,如果您能遵循这一点,说实话,它写得非常糟糕!那么请我需要你的帮助!