diff --git a/src/expand.js b/src/expand.js
new file mode 100644
index 0000000..2319adf
--- /dev/null
+++ b/src/expand.js
@@ -0,0 +1,19 @@
+'use strict';
+
+import {broadcast, getBroadcastShape} from './lib/broadcast.js';
+import {Scalar} from '../src/lib/tensor.js';
+
+/**
+ * Expand any dimension of size 1 of the input tensor to a
+ * larger size according to the new shape.
+ * @param {Tensor} input
+ * @param {Array} newShape
+ * @return {Tensor}
+ */
+
+
+export function expand(input, newShape) {
+  const inputReshape = input.shape.length === 0 ? new Scalar(input.data) : input;
+  const outputShape = getBroadcastShape(inputReshape.shape, newShape);
+  return broadcast(inputReshape, outputShape);
+}
diff --git a/test/expand_test.js b/test/expand_test.js
new file mode 100644
index 0000000..9e66610
--- /dev/null
+++ b/test/expand_test.js
@@ -0,0 +1,124 @@
+'use strict';
+
+import {expand} from '../src/expand.js';
+import {Tensor} from '../src/lib/tensor.js';
+import * as utils from './utils.js';
+
+describe('test expand', function() {
+  function testExpand(input, newShape, expected) {
+    const tensor = new Tensor(input.shape, input.data);
+    const outputTensor = expand(tensor, newShape);
+    utils.checkShape(outputTensor, expected.shape);
+    utils.checkValue(outputTensor, expected.data);
+  }
+
+  it('expand a 3D input with a 4D newShape to a 4D output.', function() {
+    const input = {
+      shape: [2, 1, 4],
+      data: [
+        1, 2, 3, 4, 5, 6, 7, 8,
+      ],
+    };
+
+    const newShape = [5, 1, 3, 4];
+
+    const expected = {
+      shape: [5, 2, 3, 4],
+      data: [
+        1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4,
+        5, 6, 7, 8, 5, 6, 7, 8, 5, 6, 7, 8,
+        1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4,
+        5, 6, 7, 8, 5, 6, 7, 8, 5, 6, 7, 8,
+        1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4,
+        5, 6, 7, 8, 5, 6, 7, 8, 5, 6, 7, 8,
+        1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4,
+        5, 6, 7, 8, 5, 6, 7, 8, 5, 6, 7, 8,
+        1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4,
+        5, 6, 7, 8, 5, 6, 7, 8, 5, 6, 7, 8,
+      ],
+    };
+    testExpand(input, newShape, expected);
+  });
+
+  it('expand a 3D input with a 2D newShape to a 3D output', function() {
+    const input = {
+      shape: [2, 1, 4],
+      data: [
+        1, 2, 3, 4, 5, 6, 7, 8,
+      ],
+    };
+
+    const newShape = [3, 1];
+
+    const expected = {
+      shape: [2, 3, 4],
+      data: [
+        1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4,
+        5, 6, 7, 8, 5, 6, 7, 8, 5, 6, 7, 8,
+      ],
+    };
+    testExpand(input, newShape, expected);
+  });
+
+  it('expand a 2D input with a 2D newShape to a 2D output', function() {
+    const input = {
+      shape: [3, 1],
+      data: [
+        1, 2, 3,
+      ],
+    };
+
+    const newShape = [3, 4];
+
+    const expected = {
+      shape: [3, 4],
+      data: [
+        1, 1, 1, 1,
+        2, 2, 2, 2,
+        3, 3, 3, 3,
+      ],
+    };
+    testExpand(input, newShape, expected);
+  });
+
+  it('expand a 0D input with a 2D newShape to a 2D output.', function() {
+    const input = {
+      shape: [],
+      data: [
+        6,
+      ],
+    };
+
+    const newShape = [2, 3];
+
+    const expected = {
+      shape: [2, 3],
+      data: [
+        6, 6, 6,
+        6, 6, 6,
+      ],
+    };
+    testExpand(input, newShape, expected);
+  });
+
+  it('expand a 2D input with a 0D newShape to a 2D output.', function() {
+    const input = {
+      shape: [2, 3],
+      data: [
+        1, 2, 3,
+        4, 5, 6,
+      ],
+    };
+
+    const newShape = [];
+
+    const expected = {
+      shape: [2, 3],
+      data: [
+        1, 2, 3,
+        4, 5, 6,
+      ],
+    };
+    testExpand(input, newShape, expected);
+  });
+});