Решить уравнение методом Ньютона (касательных)

using System;
using System.Collections.Generic;
using System.ComponentModel;
using System.Data;
using System.Drawing;
using System.Linq;
using System.Text;
using System.Threading.Tasks;
using System.Windows.Forms;
 
namespace WindowsFormsApplication4
{
class SOSAT
{
private double a, b, eps, x, fx, fa, fb, dfx, dfa, dfb, d2fx, d2fb, d2fa, s;
public SOSAT(double a, double b, double eps, double k)
{
this.a = a;
this.b = b;
this.eps = eps;
s = Math.Log10(1 + Math.Sqrt(k));
Getdix2();
}
 
public double FFx(double z)
{
return (Math.Pow(z, 2.0) * Math.Sin(z) - s);
}
 
public double dFx(double z)
{
return (2 * z * Math.Sin(z) + Math.Pow(z, 2.0) * Math.Cos(z));
}
 
public double d2Fx(double z)
{
return (4 * z * Math.Cos(z) + Math.Sin(z)*(2 - Math.Pow(z,2.0)));
}
 
public double Getdix2()
{
    if ((FFx(a) * dFx(a)) > 0)
    {
        do
        {
            x = this.a - FFx(this.a) / dFx(this.a);
            fx = FFx(x);
            fa = FFx(this.a);
            fb = FFx(this.x);
 
            dfa = dFx(this.a);
            dfx = dFx(x);
            dfb = dFx(this.x);
 
            d2fa = d2Fx(this.a);
            d2fx = d2Fx(x);
            d2fb = d2Fx(this.x);
 
            if (((fx * fa) < 0) && ((d2fa * d2fx) > 0) && ((dfa * dfx) > 0))
            {
 
                this.a = x;
                fa = fx;
                dfa = dfx;
                d2fa = d2fx;
 
 
            }
 
         /*   else
            {
 
                this.b = x;
                fb = fx;
                dfb = dfx;
                d2fb = d2fx;
 
 
            }*/
 
 
        }
        while ((Math.Abs(fx)) > this.eps);
        return x;
    }
 
    else
    {
         
            x = this.b - FFx(this.b) / dFx(this.b);
            fx = FFx(x);
            fa = FFx(this.x);
            fb = FFx(this.b);
 
            dfa = dFx(this.x);
            dfx = dFx(x);
            dfb = dFx(this.b);
 
            d2fa = d2Fx(this.x);
            d2fx = d2Fx(x);
            d2fb = d2Fx(this.b);
 
            if (((fx * fa) < 0) && ((d2fa * d2fx) > 0) && ((dfa * dfx) > 0))
            {
 
                this.b = x;
                fa = fx;
                dfa = dfx;
                d2fa = d2fx;
 
 
            }
           }
        while ((Math.Abs(fx)) > this.eps);
        return x;
    }
    }
 
}