8.3 Arithmetic Operators

The following arithmetic operators are available, and work on scalars and matrices. The element-by-element operators and functions broadcast (see Broadcasting).

x + y

Addition. If both operands are matrices, the number of rows and columns must both agree, or they must be broadcastable to the same shape.

x .+ y

Element-by-element addition. This operator is equivalent to +.

x - y

Subtraction. If both operands are matrices, the number of rows and columns of both must agree, or they must be broadcastable to the same shape.

x .- y

Element-by-element subtraction. This operator is equivalent to -.

x * y

Matrix multiplication. The number of columns of x must agree with the number of rows of y.

x .* y

Element-by-element multiplication. If both operands are matrices, the number of rows and columns must both agree, or they must be broadcastable to the same shape.

x / y

Right division. This is conceptually equivalent to the expression

(inv (y') * x')'

but it is computed without forming the inverse of y’.

If the system is not square, or if the coefficient matrix is singular, a minimum norm solution is computed.

x ./ y

Element-by-element right division.

x \ y

Left division. This is conceptually equivalent to the expression

inv (x) * y

but it is computed without forming the inverse of x.

If the system is not square, or if the coefficient matrix is singular, a minimum norm solution is computed.

x .\ y

Element-by-element left division. Each element of y is divided by each corresponding element of x.

x ^ y

x ** y

Power operator. If x and y are both scalars, this operator returns x raised to the power y. If x is a scalar and y is a square matrix, the result is computed using an eigenvalue expansion. If x is a square matrix, the result is computed by repeated multiplication if y is an integer, and by an eigenvalue expansion if y is not an integer. An error results if both x and y are matrices.

The implementation of this operator needs to be improved.

x .^ y

x .** y

Element-by-element power operator. If both operands are matrices, the number of rows and columns must both agree, or they must be broadcastable to the same shape. If several complex results are possible, the one with smallest non-negative argument (angle) is taken. This rule may return a complex root even when a real root is also possible. Use realpow, realsqrt, cbrt, or nthroot if a real result is preferred.

-x

Negation.

+x

Unary plus. This operator has no effect on the operand.

x’

Complex conjugate transpose. For real arguments, this operator is the same as the transpose operator. For complex arguments, this operator is equivalent to the expression

conj (x.')

x.’

Transpose.

Note that because Octave’s element-by-element operators begin with a ‘.’, there is a possible ambiguity for statements like

1./m

because the period could be interpreted either as part of the constant or as part of the operator. To resolve this conflict, Octave treats the expression as if you had typed

(1) ./ m

and not

(1.) / m

Although this is inconsistent with the normal behavior of Octave’s lexer, which usually prefers to break the input into tokens by preferring the longest possible match at any given point, it is more useful in this case.

: ctranspose (x)

Return the complex conjugate transpose of x.

This function and x' are equivalent.

See also: transpose.

: ldivide (x, y)

Return the element-by-element left division of x and y.

This function and x .\ y are equivalent.

See also: rdivide, mldivide, times, plus.

: minus (x, y)

This function and x - y are equivalent.

See also: plus, uminus.

: mldivide (x, y)

Return the matrix left division of x and y.

This function and x \ y are equivalent.

If the system is not square, or if the coefficient matrix is singular, a minimum norm solution is computed.

See also: mrdivide, ldivide, rdivide, linsolve.

: mpower (x, y)

Return the matrix power operation of x raised to the y power.

This function and x ^ y are equivalent.

See also: power, mtimes, plus, minus.

: mrdivide (x, y)

Return the matrix right division of x and y.

This function and x / y are equivalent.

If the system is not square, or if the coefficient matrix is singular, a minimum norm solution is computed.

See also: mldivide, rdivide, plus, minus.

: mtimes (x, y)

: mtimes (x1, x2, …)

Return the matrix multiplication product of inputs.

This function and x * y are equivalent. If more arguments are given, the multiplication is applied cumulatively from left to right:

(…((x1 * x2) * x3) * …)

See also: times, plus, minus, rdivide, mrdivide, mldivide, mpower.

: plus (x, y)

: plus (x1, x2, …)

This function and x + y are equivalent.

If more arguments are given, the summation is applied cumulatively from left to right:

(…((x1 + x2) + x3) + …)

See also: minus, uplus.

: power (x, y)

Return the element-by-element operation of x raised to the y power.

This function and x .^ y are equivalent.

If several complex results are possible, returns the one with smallest non-negative argument (angle). Use realpow, realsqrt, cbrt, or nthroot if a real result is preferred.

See also: mpower, realpow, realsqrt, cbrt, nthroot.

: rdivide (x, y)

Return the element-by-element right division of x and y.

This function and x ./ y are equivalent.

See also: ldivide, mrdivide, times, plus.

: times (x, y)

: times (x1, x2, …)

Return the element-by-element multiplication product of inputs.

This function and x .* y are equivalent. If more arguments are given, the multiplication is applied cumulatively from left to right:

(…((x1 .* x2) .* x3) .* …)

See also: mtimes, rdivide.

: transpose (x)

Return the transpose of x.

This function and x.' are equivalent.

See also: ctranspose.

: uminus (x)

This function and - x are equivalent.

See also: uplus, minus.

: uplus (x)

This function and + x are equivalent.

See also: uminus, plus.