六面体の頂点をたどる（7）

固有値が求まらなかったので、数値的に求める。
ニュートン法でやってみた。



Array.prototype.value = function(x) {

    var v = this[this.length-1];

    for(var i = this.length - 2; i >= 0; i--)

        v = v * x + this[i];

    return v;

}
Array.prototype.differential = function() {

    var f = [ ];

    for(var i = 0; i < this.length - 1; i++)

        f[i] = this[i+1] * (i + 1);

    return f;

}
var a = [ 35, 35, -47, -12, 13, 1, -1 ];        // 関数

var x = [ -2.7, -2.2, -0.6, 1.6, 2.2, 2.6 ];    // 仮の解
var eps = 1e-13;

var df = a.differential();

for(var i = 0; i < x.length; i++) {

    var v;

    while(Math.abs(v = a.value(x[i])) > eps) {

        x[i] -= v / df.value(x[i]);

    }

    WScript.Echo(x[i]);

}