mathjs 8.1.0 → 12.2.1

package/lib/esm/function/algebra/derivative.js

@@ -2,7 +2,7 @@ import { isConstantNode, typeOf } from '../../utils/is.js';
 import { factory } from '../../utils/factory.js';
 var name = 'derivative';
 var dependencies = ['typed', 'config', 'parse', 'simplify', 'equal', 'isZero', 'numeric', 'ConstantNode', 'FunctionNode', 'OperatorNode', 'ParenthesisNode', 'SymbolNode'];
-export var createDerivative = /* #__PURE__ */factory(name, dependencies, (_ref) => {
+export var createDerivative = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
     config,
@@ -17,7 +17,6 @@ export var createDerivative = /* #__PURE__ */factory(name, dependencies, (_ref)
     ParenthesisNode,
     SymbolNode
   } = _ref;
-
   /**
    * Takes the derivative of an expression expressed in parser Nodes.
    * The derivative will be taken over the supplied variable in the
@@ -30,14 +29,14 @@ export var createDerivative = /* #__PURE__ */factory(name, dependencies, (_ref)
    *
    * Syntax:
    *
-   *     derivative(expr, variable)
-   *     derivative(expr, variable, options)
+   *     math.derivative(expr, variable)
+   *     math.derivative(expr, variable, options)
    *
    * Examples:
    *
-   *     math.derivative('x^2', 'x')                     // Node {2 * x}
-   *     math.derivative('x^2', 'x', {simplify: false})  // Node {2 * 1 * x ^ (2 - 1)
-   *     math.derivative('sin(2x)', 'x'))                // Node {2 * cos(2 * x)}
+   *     math.derivative('x^2', 'x')                     // Node '2 * x'
+   *     math.derivative('x^2', 'x', {simplify: false})  // Node '2 * 1 * x ^ (2 - 1)'
+   *     math.derivative('sin(2x)', 'x'))                // Node '2 * cos(2 * x)'
    *     math.derivative('2*x', 'x').evaluate()          // number 2
    *     math.derivative('x^2', 'x').evaluate({x: 4})    // number 8
    *     const f = math.parse('x^2')
@@ -56,38 +55,23 @@ export var createDerivative = /* #__PURE__ */factory(name, dependencies, (_ref)
    *                         be simplified.
    * @return {ConstantNode | SymbolNode | ParenthesisNode | FunctionNode | OperatorNode}    The derivative of `expr`
    */
-  var derivative = typed('derivative', {
-    'Node, SymbolNode, Object': function NodeSymbolNodeObject(expr, variable, options) {
-      var constNodes = {};
-      constTag(constNodes, expr, variable.name);
-
-      var res = _derivative(expr, constNodes);
-
-      return options.simplify ? simplify(res) : res;
-    },
-    'Node, SymbolNode': function NodeSymbolNode(expr, variable) {
-      return this(expr, variable, {
-        simplify: true
-      });
-    },
-    'string, SymbolNode': function stringSymbolNode(expr, variable) {
-      return this(parse(expr), variable);
-    },
-    'string, SymbolNode, Object': function stringSymbolNodeObject(expr, variable, options) {
-      return this(parse(expr), variable, options);
-    },
-    'string, string': function stringString(expr, variable) {
-      return this(parse(expr), parse(variable));
-    },
-    'string, string, Object': function stringStringObject(expr, variable, options) {
-      return this(parse(expr), parse(variable), options);
-    },
-    'Node, string': function NodeString(expr, variable) {
-      return this(expr, parse(variable));
-    },
-    'Node, string, Object': function NodeStringObject(expr, variable, options) {
-      return this(expr, parse(variable), options);
-    } // TODO: replace the 8 signatures above with 4 as soon as typed-function supports optional arguments
+  function plainDerivative(expr, variable) {
+    var options = arguments.length > 2 && arguments[2] !== undefined ? arguments[2] : {
+      simplify: true
+    };
+    var constNodes = {};
+    constTag(constNodes, expr, variable.name);
+    var res = _derivative(expr, constNodes);
+    return options.simplify ? simplify(res) : res;
+  }
+  typed.addConversion({
+    from: 'identifier',
+    to: 'SymbolNode',
+    convert: parse
+  });
+  var derivative = typed(name, {
+    'Node, SymbolNode': plainDerivative,
+    'Node, SymbolNode, Object': plainDerivative
 
     /* TODO: implement and test syntax with order of derivatives -> implement as an option {order: number}
     'Node, SymbolNode, ConstantNode': function (expr, variable, {order}) {
@@ -100,22 +84,25 @@ export var createDerivative = /* #__PURE__ */factory(name, dependencies, (_ref)
       return res
     }
     */
-
+  });
+  typed.removeConversion({
+    from: 'identifier',
+    to: 'SymbolNode',
+    convert: parse
   });
   derivative._simplify = true;
-
   derivative.toTex = function (deriv) {
     return _derivTex.apply(null, deriv.args);
-  }; // FIXME: move the toTex method of derivative to latex.js. Difficulty is that it relies on parse.
-  // NOTE: the optional "order" parameter here is currently unused
-
+  };
 
+  // FIXME: move the toTex method of derivative to latex.js. Difficulty is that it relies on parse.
+  // NOTE: the optional "order" parameter here is currently unused
   var _derivTex = typed('_derivTex', {
     'Node, SymbolNode': function NodeSymbolNode(expr, x) {
       if (isConstantNode(expr) && typeOf(expr.value) === 'string') {
         return _derivTex(parse(expr.value).toString(), x.toString(), 1);
       } else {
-        return _derivTex(expr.toString(), x.toString(), 1);
+        return _derivTex(expr.toTex(), x.toString(), 1);
       }
     },
     'Node, ConstantNode': function NodeConstantNode(expr, x) {
@@ -130,16 +117,15 @@ export var createDerivative = /* #__PURE__ */factory(name, dependencies, (_ref)
     },
     'string, string, number': function stringStringNumber(expr, x, order) {
       var d;
-
       if (order === 1) {
         d = '{d\\over d' + x + '}';
       } else {
         d = '{d^{' + order + '}\\over d' + x + '^{' + order + '}}';
       }
-
       return d + "\\left[".concat(expr, "\\right]");
     }
   });
+
   /**
    * Does a depth-first search on the expression tree to identify what Nodes
    * are constants (e.g. 2 + 2), and stores the ones that are constants in
@@ -155,8 +141,6 @@ export var createDerivative = /* #__PURE__ */factory(name, dependencies, (_ref)
    * @return {boolean}  if node is constant
    */
   // TODO: can we rewrite constTag into a pure function?
-
-
   var constTag = typed('constTag', {
     'Object, ConstantNode, string': function ObjectConstantNodeString(constNodes, node) {
       constNodes[node] = true;
@@ -169,7 +153,6 @@ export var createDerivative = /* #__PURE__ */factory(name, dependencies, (_ref)
         constNodes[node] = true;
         return true;
       }
-
       return false;
     },
     'Object, ParenthesisNode, string': function ObjectParenthesisNodeString(constNodes, node, varName) {
@@ -180,26 +163,23 @@ export var createDerivative = /* #__PURE__ */factory(name, dependencies, (_ref)
         constNodes[node] = true;
         return true;
       }
-
       return constTag(constNodes, node.expr, varName);
     },
     'Object, FunctionNode | OperatorNode, string': function ObjectFunctionNodeOperatorNodeString(constNodes, node, varName) {
       if (node.args.length > 0) {
         var isConst = constTag(constNodes, node.args[0], varName);
-
         for (var i = 1; i < node.args.length; ++i) {
           isConst = constTag(constNodes, node.args[i], varName) && isConst;
         }
-
         if (isConst) {
           constNodes[node] = true;
           return true;
         }
       }
-
       return false;
     }
   });
+
   /**
    * Applies differentiation rules.
    *
@@ -207,7 +187,6 @@ export var createDerivative = /* #__PURE__ */factory(name, dependencies, (_ref)
    * @param  {Object} constNodes  Holds the nodes that are constant
    * @return {ConstantNode | SymbolNode | ParenthesisNode | FunctionNode | OperatorNode}    The derivative of `expr`
    */
-
   var _derivative = typed('_derivative', {
     'ConstantNode, Object': function ConstantNodeObject(node) {
       return createConstantNode(0);
@@ -216,7 +195,6 @@ export var createDerivative = /* #__PURE__ */factory(name, dependencies, (_ref)
       if (constNodes[node] !== undefined) {
         return createConstantNode(0);
       }
-
       return createConstantNode(1);
     },
     'ParenthesisNode, Object': function ParenthesisNodeObject(node, constNodes) {
@@ -226,33 +204,24 @@ export var createDerivative = /* #__PURE__ */factory(name, dependencies, (_ref)
       if (constNodes[node] !== undefined) {
         return createConstantNode(0);
       }
-
       return _derivative(node.expr, constNodes);
     },
     'FunctionNode, Object': function FunctionNodeObject(node, constNodes) {
-      if (node.args.length !== 1) {
-        funcArgsCheck(node);
-      }
-
       if (constNodes[node] !== undefined) {
         return createConstantNode(0);
       }
-
       var arg0 = node.args[0];
       var arg1;
       var div = false; // is output a fraction?
-
       var negative = false; // is output negative?
 
       var funcDerivative;
-
       switch (node.name) {
         case 'cbrt':
           // d/dx(cbrt(x)) = 1 / (3x^(2/3))
           div = true;
           funcDerivative = new OperatorNode('*', 'multiply', [createConstantNode(3), new OperatorNode('^', 'pow', [arg0, new OperatorNode('/', 'divide', [createConstantNode(2), createConstantNode(3)])])]);
           break;
-
         case 'sqrt':
         case 'nthRoot':
           // d/dx(sqrt(x)) = 1 / (2*sqrt(x))
@@ -261,19 +230,16 @@ export var createDerivative = /* #__PURE__ */factory(name, dependencies, (_ref)
             funcDerivative = new OperatorNode('*', 'multiply', [createConstantNode(2), new FunctionNode('sqrt', [arg0])]);
           } else if (node.args.length === 2) {
             // Rearrange from nthRoot(x, a) -> x^(1/a)
-            arg1 = new OperatorNode('/', 'divide', [createConstantNode(1), node.args[1]]); // Is a variable?
+            arg1 = new OperatorNode('/', 'divide', [createConstantNode(1), node.args[1]]);
 
+            // Is a variable?
             constNodes[arg1] = constNodes[node.args[1]];
             return _derivative(new OperatorNode('^', 'pow', [arg0, arg1]), constNodes);
           }
-
           break;
-
         case 'log10':
           arg1 = createConstantNode(10);
-
         /* fall through! */
-
         case 'log':
           if (!arg1 && node.args.length === 1) {
             // d/dx(log(x)) = 1 / x
@@ -287,175 +253,146 @@ export var createDerivative = /* #__PURE__ */factory(name, dependencies, (_ref)
             // d/dx(log(f(x), g(x))) = d/dx(log(f(x)) / log(g(x)))
             return _derivative(new OperatorNode('/', 'divide', [new FunctionNode('log', [arg0]), new FunctionNode('log', [node.args[1]])]), constNodes);
           }
-
           break;
-
         case 'pow':
-          constNodes[arg1] = constNodes[node.args[1]]; // Pass to pow operator node parser
-
-          return _derivative(new OperatorNode('^', 'pow', [arg0, node.args[1]]), constNodes);
-
+          if (node.args.length === 2) {
+            constNodes[arg1] = constNodes[node.args[1]];
+            // Pass to pow operator node parser
+            return _derivative(new OperatorNode('^', 'pow', [arg0, node.args[1]]), constNodes);
+          }
+          break;
         case 'exp':
           // d/dx(e^x) = e^x
           funcDerivative = new FunctionNode('exp', [arg0.clone()]);
           break;
-
         case 'sin':
           // d/dx(sin(x)) = cos(x)
           funcDerivative = new FunctionNode('cos', [arg0.clone()]);
           break;
-
         case 'cos':
           // d/dx(cos(x)) = -sin(x)
           funcDerivative = new OperatorNode('-', 'unaryMinus', [new FunctionNode('sin', [arg0.clone()])]);
           break;
-
         case 'tan':
           // d/dx(tan(x)) = sec(x)^2
           funcDerivative = new OperatorNode('^', 'pow', [new FunctionNode('sec', [arg0.clone()]), createConstantNode(2)]);
           break;
-
         case 'sec':
           // d/dx(sec(x)) = sec(x)tan(x)
           funcDerivative = new OperatorNode('*', 'multiply', [node, new FunctionNode('tan', [arg0.clone()])]);
           break;
-
         case 'csc':
           // d/dx(csc(x)) = -csc(x)cot(x)
           negative = true;
           funcDerivative = new OperatorNode('*', 'multiply', [node, new FunctionNode('cot', [arg0.clone()])]);
           break;
-
         case 'cot':
           // d/dx(cot(x)) = -csc(x)^2
           negative = true;
           funcDerivative = new OperatorNode('^', 'pow', [new FunctionNode('csc', [arg0.clone()]), createConstantNode(2)]);
           break;
-
         case 'asin':
           // d/dx(asin(x)) = 1 / sqrt(1 - x^2)
           div = true;
           funcDerivative = new FunctionNode('sqrt', [new OperatorNode('-', 'subtract', [createConstantNode(1), new OperatorNode('^', 'pow', [arg0.clone(), createConstantNode(2)])])]);
           break;
-
         case 'acos':
           // d/dx(acos(x)) = -1 / sqrt(1 - x^2)
           div = true;
           negative = true;
           funcDerivative = new FunctionNode('sqrt', [new OperatorNode('-', 'subtract', [createConstantNode(1), new OperatorNode('^', 'pow', [arg0.clone(), createConstantNode(2)])])]);
           break;
-
         case 'atan':
           // d/dx(atan(x)) = 1 / (x^2 + 1)
           div = true;
           funcDerivative = new OperatorNode('+', 'add', [new OperatorNode('^', 'pow', [arg0.clone(), createConstantNode(2)]), createConstantNode(1)]);
           break;
-
         case 'asec':
           // d/dx(asec(x)) = 1 / (|x|*sqrt(x^2 - 1))
           div = true;
           funcDerivative = new OperatorNode('*', 'multiply', [new FunctionNode('abs', [arg0.clone()]), new FunctionNode('sqrt', [new OperatorNode('-', 'subtract', [new OperatorNode('^', 'pow', [arg0.clone(), createConstantNode(2)]), createConstantNode(1)])])]);
           break;
-
         case 'acsc':
           // d/dx(acsc(x)) = -1 / (|x|*sqrt(x^2 - 1))
           div = true;
           negative = true;
           funcDerivative = new OperatorNode('*', 'multiply', [new FunctionNode('abs', [arg0.clone()]), new FunctionNode('sqrt', [new OperatorNode('-', 'subtract', [new OperatorNode('^', 'pow', [arg0.clone(), createConstantNode(2)]), createConstantNode(1)])])]);
           break;
-
         case 'acot':
           // d/dx(acot(x)) = -1 / (x^2 + 1)
           div = true;
           negative = true;
           funcDerivative = new OperatorNode('+', 'add', [new OperatorNode('^', 'pow', [arg0.clone(), createConstantNode(2)]), createConstantNode(1)]);
           break;
-
         case 'sinh':
           // d/dx(sinh(x)) = cosh(x)
           funcDerivative = new FunctionNode('cosh', [arg0.clone()]);
           break;
-
         case 'cosh':
           // d/dx(cosh(x)) = sinh(x)
           funcDerivative = new FunctionNode('sinh', [arg0.clone()]);
           break;
-
         case 'tanh':
           // d/dx(tanh(x)) = sech(x)^2
           funcDerivative = new OperatorNode('^', 'pow', [new FunctionNode('sech', [arg0.clone()]), createConstantNode(2)]);
           break;
-
         case 'sech':
           // d/dx(sech(x)) = -sech(x)tanh(x)
           negative = true;
           funcDerivative = new OperatorNode('*', 'multiply', [node, new FunctionNode('tanh', [arg0.clone()])]);
           break;
-
         case 'csch':
           // d/dx(csch(x)) = -csch(x)coth(x)
           negative = true;
           funcDerivative = new OperatorNode('*', 'multiply', [node, new FunctionNode('coth', [arg0.clone()])]);
           break;
-
         case 'coth':
           // d/dx(coth(x)) = -csch(x)^2
           negative = true;
           funcDerivative = new OperatorNode('^', 'pow', [new FunctionNode('csch', [arg0.clone()]), createConstantNode(2)]);
           break;
-
         case 'asinh':
           // d/dx(asinh(x)) = 1 / sqrt(x^2 + 1)
           div = true;
           funcDerivative = new FunctionNode('sqrt', [new OperatorNode('+', 'add', [new OperatorNode('^', 'pow', [arg0.clone(), createConstantNode(2)]), createConstantNode(1)])]);
           break;
-
         case 'acosh':
           // d/dx(acosh(x)) = 1 / sqrt(x^2 - 1); XXX potentially only for x >= 1 (the real spectrum)
           div = true;
           funcDerivative = new FunctionNode('sqrt', [new OperatorNode('-', 'subtract', [new OperatorNode('^', 'pow', [arg0.clone(), createConstantNode(2)]), createConstantNode(1)])]);
           break;
-
         case 'atanh':
           // d/dx(atanh(x)) = 1 / (1 - x^2)
           div = true;
           funcDerivative = new OperatorNode('-', 'subtract', [createConstantNode(1), new OperatorNode('^', 'pow', [arg0.clone(), createConstantNode(2)])]);
           break;
-
         case 'asech':
           // d/dx(asech(x)) = -1 / (x*sqrt(1 - x^2))
           div = true;
           negative = true;
           funcDerivative = new OperatorNode('*', 'multiply', [arg0.clone(), new FunctionNode('sqrt', [new OperatorNode('-', 'subtract', [createConstantNode(1), new OperatorNode('^', 'pow', [arg0.clone(), createConstantNode(2)])])])]);
           break;
-
         case 'acsch':
           // d/dx(acsch(x)) = -1 / (|x|*sqrt(x^2 + 1))
           div = true;
           negative = true;
           funcDerivative = new OperatorNode('*', 'multiply', [new FunctionNode('abs', [arg0.clone()]), new FunctionNode('sqrt', [new OperatorNode('+', 'add', [new OperatorNode('^', 'pow', [arg0.clone(), createConstantNode(2)]), createConstantNode(1)])])]);
           break;
-
         case 'acoth':
           // d/dx(acoth(x)) = -1 / (1 - x^2)
           div = true;
           negative = true;
           funcDerivative = new OperatorNode('-', 'subtract', [createConstantNode(1), new OperatorNode('^', 'pow', [arg0.clone(), createConstantNode(2)])]);
           break;
-
         case 'abs':
           // d/dx(abs(x)) = abs(x)/x
           funcDerivative = new OperatorNode('/', 'divide', [new FunctionNode(new SymbolNode('abs'), [arg0.clone()]), arg0.clone()]);
           break;
-
         case 'gamma': // Needs digamma function, d/dx(gamma(x)) = gamma(x)digamma(x)
-
         default:
-          throw new Error('Function "' + node.name + '" is not supported by derivative, or a wrong number of arguments is passed');
+          throw new Error('Cannot process function "' + node.name + '" in derivative: ' + 'the function is not supported, undefined, or the number of arguments passed to it are not supported');
       }
-
       var op, func;
-
       if (div) {
         op = '/';
         func = 'divide';
@@ -463,49 +400,42 @@ export var createDerivative = /* #__PURE__ */factory(name, dependencies, (_ref)
         op = '*';
         func = 'multiply';
       }
+
       /* Apply chain rule to all functions:
          F(x)  = f(g(x))
          F'(x) = g'(x)*f'(g(x)) */
-
-
       var chainDerivative = _derivative(arg0, constNodes);
-
       if (negative) {
         chainDerivative = new OperatorNode('-', 'unaryMinus', [chainDerivative]);
       }
-
       return new OperatorNode(op, func, [chainDerivative, funcDerivative]);
     },
     'OperatorNode, Object': function OperatorNodeObject(node, constNodes) {
       if (constNodes[node] !== undefined) {
         return createConstantNode(0);
       }
-
       if (node.op === '+') {
         // d/dx(sum(f(x)) = sum(f'(x))
         return new OperatorNode(node.op, node.fn, node.args.map(function (arg) {
           return _derivative(arg, constNodes);
         }));
       }
-
       if (node.op === '-') {
         // d/dx(+/-f(x)) = +/-f'(x)
         if (node.isUnary()) {
           return new OperatorNode(node.op, node.fn, [_derivative(node.args[0], constNodes)]);
-        } // Linearity of differentiation, d/dx(f(x) +/- g(x)) = f'(x) +/- g'(x)
-
+        }
 
+        // Linearity of differentiation, d/dx(f(x) +/- g(x)) = f'(x) +/- g'(x)
         if (node.isBinary()) {
           return new OperatorNode(node.op, node.fn, [_derivative(node.args[0], constNodes), _derivative(node.args[1], constNodes)]);
         }
       }
-
       if (node.op === '*') {
         // d/dx(c*f(x)) = c*f'(x)
         var constantTerms = node.args.filter(function (arg) {
           return constNodes[arg] !== undefined;
         });
-
         if (constantTerms.length > 0) {
           var nonConstantTerms = node.args.filter(function (arg) {
             return constNodes[arg] === undefined;
@@ -513,96 +443,68 @@ export var createDerivative = /* #__PURE__ */factory(name, dependencies, (_ref)
           var nonConstantNode = nonConstantTerms.length === 1 ? nonConstantTerms[0] : new OperatorNode('*', 'multiply', nonConstantTerms);
           var newArgs = constantTerms.concat(_derivative(nonConstantNode, constNodes));
           return new OperatorNode('*', 'multiply', newArgs);
-        } // Product Rule, d/dx(f(x)*g(x)) = f'(x)*g(x) + f(x)*g'(x)
-
+        }
 
+        // Product Rule, d/dx(f(x)*g(x)) = f'(x)*g(x) + f(x)*g'(x)
         return new OperatorNode('+', 'add', node.args.map(function (argOuter) {
           return new OperatorNode('*', 'multiply', node.args.map(function (argInner) {
             return argInner === argOuter ? _derivative(argInner, constNodes) : argInner.clone();
           }));
         }));
       }
-
       if (node.op === '/' && node.isBinary()) {
         var arg0 = node.args[0];
-        var arg1 = node.args[1]; // d/dx(f(x) / c) = f'(x) / c
+        var arg1 = node.args[1];
 
+        // d/dx(f(x) / c) = f'(x) / c
         if (constNodes[arg1] !== undefined) {
           return new OperatorNode('/', 'divide', [_derivative(arg0, constNodes), arg1]);
-        } // Reciprocal Rule, d/dx(c / f(x)) = -c(f'(x)/f(x)^2)
-
+        }
 
+        // Reciprocal Rule, d/dx(c / f(x)) = -c(f'(x)/f(x)^2)
         if (constNodes[arg0] !== undefined) {
           return new OperatorNode('*', 'multiply', [new OperatorNode('-', 'unaryMinus', [arg0]), new OperatorNode('/', 'divide', [_derivative(arg1, constNodes), new OperatorNode('^', 'pow', [arg1.clone(), createConstantNode(2)])])]);
-        } // Quotient rule, d/dx(f(x) / g(x)) = (f'(x)g(x) - f(x)g'(x)) / g(x)^2
-
+        }
 
+        // Quotient rule, d/dx(f(x) / g(x)) = (f'(x)g(x) - f(x)g'(x)) / g(x)^2
         return new OperatorNode('/', 'divide', [new OperatorNode('-', 'subtract', [new OperatorNode('*', 'multiply', [_derivative(arg0, constNodes), arg1.clone()]), new OperatorNode('*', 'multiply', [arg0.clone(), _derivative(arg1, constNodes)])]), new OperatorNode('^', 'pow', [arg1.clone(), createConstantNode(2)])]);
       }
-
       if (node.op === '^' && node.isBinary()) {
         var _arg = node.args[0];
         var _arg2 = node.args[1];
-
         if (constNodes[_arg] !== undefined) {
           // If is secretly constant; 0^f(x) = 1 (in JS), 1^f(x) = 1
           if (isConstantNode(_arg) && (isZero(_arg.value) || equal(_arg.value, 1))) {
             return createConstantNode(0);
-          } // d/dx(c^f(x)) = c^f(x)*ln(c)*f'(x)
-
+          }
 
+          // d/dx(c^f(x)) = c^f(x)*ln(c)*f'(x)
           return new OperatorNode('*', 'multiply', [node, new OperatorNode('*', 'multiply', [new FunctionNode('log', [_arg.clone()]), _derivative(_arg2.clone(), constNodes)])]);
         }
-
         if (constNodes[_arg2] !== undefined) {
           if (isConstantNode(_arg2)) {
             // If is secretly constant; f(x)^0 = 1 -> d/dx(1) = 0
             if (isZero(_arg2.value)) {
               return createConstantNode(0);
-            } // Ignore exponent; f(x)^1 = f(x)
-
-
+            }
+            // Ignore exponent; f(x)^1 = f(x)
             if (equal(_arg2.value, 1)) {
               return _derivative(_arg, constNodes);
             }
-          } // Elementary Power Rule, d/dx(f(x)^c) = c*f'(x)*f(x)^(c-1)
-
+          }
 
+          // Elementary Power Rule, d/dx(f(x)^c) = c*f'(x)*f(x)^(c-1)
           var powMinusOne = new OperatorNode('^', 'pow', [_arg.clone(), new OperatorNode('-', 'subtract', [_arg2, createConstantNode(1)])]);
           return new OperatorNode('*', 'multiply', [_arg2.clone(), new OperatorNode('*', 'multiply', [_derivative(_arg, constNodes), powMinusOne])]);
-        } // Functional Power Rule, d/dx(f^g) = f^g*[f'*(g/f) + g'ln(f)]
-
+        }
 
+        // Functional Power Rule, d/dx(f^g) = f^g*[f'*(g/f) + g'ln(f)]
         return new OperatorNode('*', 'multiply', [new OperatorNode('^', 'pow', [_arg.clone(), _arg2.clone()]), new OperatorNode('+', 'add', [new OperatorNode('*', 'multiply', [_derivative(_arg, constNodes), new OperatorNode('/', 'divide', [_arg2.clone(), _arg.clone()])]), new OperatorNode('*', 'multiply', [_derivative(_arg2, constNodes), new FunctionNode('log', [_arg.clone()])])])]);
       }
-
-      throw new Error('Operator "' + node.op + '" is not supported by derivative, or a wrong number of arguments is passed');
+      throw new Error('Cannot process operator "' + node.op + '" in derivative: ' + 'the operator is not supported, undefined, or the number of arguments passed to it are not supported');
     }
   });
-  /**
-   * Ensures the number of arguments for a function are correct,
-   * and will throw an error otherwise.
-   *
-   * @param {FunctionNode} node
-   */
-
-
-  function funcArgsCheck(node) {
-    // TODO add min, max etc
-    if ((node.name === 'log' || node.name === 'nthRoot' || node.name === 'pow') && node.args.length === 2) {
-      return;
-    } // There should be an incorrect number of arguments if we reach here
-    // Change all args to constants to avoid unidentified
-    // symbol error when compiling function
 
-
-    for (var i = 0; i < node.args.length; ++i) {
-      node.args[i] = createConstantNode(0);
-    }
-
-    node.compile().evaluate();
-    throw new Error('Expected TypeError, but none found');
-  }
   /**
    * Helper function to create a constant node with a specific type
    * (number, BigNumber, Fraction)
@@ -610,11 +512,8 @@ export var createDerivative = /* #__PURE__ */factory(name, dependencies, (_ref)
    * @param {string} [valueType]
    * @return {ConstantNode}
    */
-
-
   function createConstantNode(value, valueType) {
     return new ConstantNode(numeric(value, valueType || config.number));
   }
-
   return derivative;
 });

package/lib/esm/function/algebra/leafCount.js

@@ -0,0 +1,54 @@
+import { factory } from '../../utils/factory.js';
+var name = 'leafCount';
+var dependencies = ['parse', 'typed'];
+export var createLeafCount = /* #__PURE__ */factory(name, dependencies, _ref => {
+  var {
+    parse,
+    typed
+  } = _ref;
+  // This does the real work, but we don't have to recurse through
+  // a typed call if we separate it out
+  function countLeaves(node) {
+    var count = 0;
+    node.forEach(n => {
+      count += countLeaves(n);
+    });
+    return count || 1;
+  }
+
+  /**
+   * Gives the number of "leaf nodes" in the parse tree of the given expression
+   * A leaf node is one that has no subexpressions, essentially either a
+   * symbol or a constant. Note that `5!` has just one leaf, the '5'; the
+   * unary factorial operator does not add a leaf. On the other hand,
+   * function symbols do add leaves, so `sin(x)/cos(x)` has four leaves.
+   *
+   * The `simplify()` function should generally not increase the `leafCount()`
+   * of an expression, although currently there is no guarantee that it never
+   * does so. In many cases, `simplify()` reduces the leaf count.
+   *
+   * Syntax:
+   *
+   *     math.leafCount(expr)
+   *
+   * Examples:
+   *
+   *     math.leafCount('x') // 1
+   *     math.leafCount(math.parse('a*d-b*c')) // 4
+   *     math.leafCount('[a,b;c,d][0,1]') // 6
+   *
+   * See also:
+   *
+   *     simplify
+   *
+   * @param {Node|string} expr    The expression to count the leaves of
+   *
+   * @return {number}  The number of leaves of `expr`
+   *
+   */
+  return typed(name, {
+    Node: function Node(expr) {
+      return countLeaves(expr);
+    }
+  });
+});