Make Jacobian types offered by dlsode.f accessible by lsode (bug #31626). * liboctave/numeric/LSODE-opts.in: Add options "jacobian type", "lower jacobian subdiagonals", and "upper jacobian subdiagonals". * liboctave/numeric/LSODE.cc (file scope, lsode_j, LSODE::do_integrate (double)): Handle new configurable Jacobian types. * build-aux/mk-opts.pl: Don't implicitly convert to integer in condition.

function r = nested_test_2 ()
  a = 3;
  b = 5;
  c = 7;
  f1 = @(f, x) f(x) + a;
  f2 = @(y) f1 (@(z) z^2 + b*y, y) + c;
  r = f2 (2);
endfunction