Дискретный метод Ньютона для поиска минимума(перевод программного кода из Java в MATLAB)

package org.mtuci.mm.zero_max_min;
 
// ПРИМЕНЕНИЕ ДИСКРЕТНОГО МЕТОДА НЬЮТОНА ДЛЯ ПОИСКА МИНИМУМА
public class Bond 
{
  static final double e2 = 14.4;
  static final double v0 = 1090;
  static final double r0 = 0.33;
  
  public static void main(String argv[]) 
  {
    double del = 1e-12;
    double r   = 3;
    double dr  = 0.1;
    
    int n = 200;
    r = secant2(n, del, r, dr);
    
    System.out.println("The bond length is " + r + " angstroms");
    System.out.println("The value is: " + g(r));
  }
 
// Method to carry out the secant search for the
// minimum of g(x) via the root of f(x)=-dg(x)/dx=0.
  public static double secant2(int n, double del, double x, double dx) 
  {
    int     k = 0;
    double x1 = x+dx;
    double g0 = g(x);
    double g1 = g(x1);
    
    if (g1 > g0)
        x1 = x-dx;
    
    while ((Math.abs(dx) > del) && (k < n)) 
    {
      double d = f(x1) - f(x);
      dx = -(x1 - x) * f(x1) / d;
      x  = x1;
      x1 += dx;
      k++;
    }
    
    if (k == n)
        System.out.println("Convergence not found after " + n + " iterations");
    
    return x1;
  }
 
  // Method to provide function g(x)=-e2/x+v0*exp(-x/r0).
  public static double g(double x) 
  {
      return -e2 / x + v0 * Math.exp(-x / r0);
  }
 
  // Method to provide function f(x)=-dg(x)/dx.
  public static double f(double x) 
  {
    return -e2 / (x * x) + v0 * Math.exp(-x / r0) / r0;
  }
}