@pawells/math-extended 1.0.5 → 2.0.0

package/build/vectors/core.js

@@ -13,9 +13,11 @@ import { AssertVector, AssertVectors, AssertVectorValue, VectorError, AssertVect
  * @returns A new vector with identical components
  *
  * @example
- * const original = [1, 2, 3];
- * const copy = VectorClone(original);
- * copy[0] = 10; // original remains unchanged
+     * ```typescript
+     * const original = [1, 2, 3];
+     * const copy = VectorClone(original);
+     * copy[0] = 10; // original remains unchanged
+     * ```
  */
 export function VectorClone(vector) {
     AssertVector(vector);
@@ -32,10 +34,12 @@ export function VectorClone(vector) {
  * @returns True if vectors are equal within tolerance, false otherwise
  *
  * @example
- * const a = [1.0001, 2.0001];
- * const b = [1.0002, 2.0002];
- * const exactlyEqual = VectorEquals(a, b); // false
- * const approximatelyEqual = VectorEquals(a, b, 0.001); // true
+     * ```typescript
+     * const a = [1.0001, 2.0001];
+     * const b = [1.0002, 2.0002];
+     * const exactlyEqual = VectorEquals(a, b); // false
+     * const approximatelyEqual = VectorEquals(a, b, 0.001); // true
+     * ```
  */
 export function VectorEquals(a, b, tolerance = 0) {
     AssertVectors([a, b]);
@@ -64,9 +68,11 @@ export function VectorEquals(a, b, tolerance = 0) {
  * @returns String representation of the vector
  *
  * @example
- * const vec = [1, 2, 3];
- * const parens = VectorToString(vec, 'parens'); // "(1, 2, 3)"
- * const brackets = VectorToString(vec, 'brackets'); // "[1, 2, 3]"
+     * ```typescript
+     * const vec = [1, 2, 3];
+     * const parens = VectorToString(vec, 'parens'); // "(1, 2, 3)"
+     * const brackets = VectorToString(vec, 'brackets'); // "[1, 2, 3]"
+     * ```
  */
 export function VectorToString(vector, style = 'parens') {
     AssertVector(vector);
@@ -74,12 +80,7 @@ export function VectorToString(vector, style = 'parens') {
     if (style === 'parens') {
         return `(${components})`;
     }
-    else if (style === 'brackets') {
-        return `[${components}]`;
-    }
-    else {
-        throw new Error(`Invalid style: ${style}. Use 'parens' or 'brackets'.`);
-    }
+    return `[${components}]`;
 }
 /**
  * Performs component-wise addition of two vectors.
@@ -91,9 +92,11 @@ export function VectorToString(vector, style = 'parens') {
  * @returns New vector where each component is the sum of corresponding components
  *
  * @example
- * const position = [10, 20, 30];
- * const velocity = [1, -2, 0.5];
- * const newPosition = VectorAdd(position, velocity); // [11, 18, 30.5]
+     * ```typescript
+     * const position = [10, 20, 30];
+     * const velocity = [1, -2, 0.5];
+     * const newPosition = VectorAdd(position, velocity); // [11, 18, 30.5]
+     * ```
  */
 export function VectorAdd(a, b) {
     AssertVectors([a, b]);
@@ -117,9 +120,11 @@ export function VectorAdd(a, b) {
  * @returns New vector where each component is the difference of corresponding components
  *
  * @example
- * const target = [100, 50, 0];
- * const current = [80, 30, 0];
- * const direction = VectorSubtract(target, current); // [20, 20, 0]
+     * ```typescript
+     * const target = [100, 50, 0];
+     * const current = [80, 30, 0];
+     * const direction = VectorSubtract(target, current); // [20, 20, 0]
+     * ```
  */
 export function VectorSubtract(a, b) {
     AssertVectors([a, b]);
@@ -143,10 +148,12 @@ export function VectorSubtract(a, b) {
  * @returns New vector with multiplied components
  *
  * @example
- * const velocity = [10, 5, 0];
- * const scaled = VectorMultiply(velocity, 2); // [20, 10, 0] - scalar multiplication
- * const factors = [1, -1, 0.5];
- * const componentWise = VectorMultiply(velocity, factors); // [10, -5, 0] - component-wise
+     * ```typescript
+     * const velocity = [10, 5, 0];
+     * const scaled = VectorMultiply(velocity, 2); // [20, 10, 0] - scalar multiplication
+     * const factors = [1, -1, 0.5];
+     * const componentWise = VectorMultiply(velocity, factors); // [10, -5, 0] - component-wise
+     * ```
  */
 export function VectorMultiply(a, b) {
     const result = [];
@@ -180,9 +187,11 @@ export function VectorMultiply(a, b) {
  * @returns The straight-line distance between the two points represented by the vectors
  *
  * @example
- * const pointA = [0, 0, 0];
- * const pointB = [3, 4, 0];
- * const distance = VectorDistance(pointA, pointB); // 5.0 (3-4-5 triangle)
+     * ```typescript
+     * const pointA = [0, 0, 0];
+     * const pointB = [3, 4, 0];
+     * const distance = VectorDistance(pointA, pointB); // 5.0 (3-4-5 triangle)
+     * ```
  */
 export function VectorDistance(a, b) {
     return Math.sqrt(VectorDistanceSquared(a, b));
@@ -197,9 +206,11 @@ export function VectorDistance(a, b) {
  * @returns The squared distance between vectors
  *
  * @example
- * const pointA = [1, 1];
- * const pointB = [4, 5];
- * const distSq = VectorDistanceSquared(pointA, pointB); // 25 (faster than distance comparison)
+     * ```typescript
+     * const pointA = [1, 1];
+     * const pointB = [4, 5];
+     * const distSq = VectorDistanceSquared(pointA, pointB); // 25 (faster than distance comparison)
+     * ```
  */
 export function VectorDistanceSquared(a, b) {
     AssertVectors([a, b]);
@@ -209,7 +220,8 @@ export function VectorDistanceSquared(a, b) {
         AssertVectorValue(av, {});
         const bv = b[i];
         AssertVectorValue(bv, {});
-        sum += Math.pow(bv - av, 2);
+        const diff = bv - av;
+        sum += diff * diff;
     }
     return sum;
 }
@@ -223,12 +235,13 @@ export function VectorDistanceSquared(a, b) {
  * @returns The dot product (scalar value)
  *
  * @example
- * const forward = [0, 0, 1];
- * const direction = [0, 0, 2];
- * const dot = VectorDot(forward, direction); // 2 (same direction)
- *
- * const perpendicular = [1, 0, 0];
- * const dotPerp = VectorDot(forward, perpendicular); // 0 (perpendicular)
+     * ```typescript
+     * const forward = [0, 0, 1];
+     * const direction = [0, 0, 2];
+     * const dot = VectorDot(forward, direction); // 2 (same direction)
+     * const perpendicular = [1, 0, 0];
+     * const dotPerp = VectorDot(forward, perpendicular); // 0 (perpendicular)
+     * ```
  */
 export function VectorDot(a, b) {
     AssertVectors([a, b]);
@@ -253,11 +266,12 @@ export function VectorDot(a, b) {
  * @throws {VectorError} If the vector is zero or has infinite magnitude
  *
  * @example
- * const vector = [3, 4, 0];
- * const normalized = VectorNormalize(vector); // [0.6, 0.8, 0] (magnitude = 1)
- *
- * const direction = [10, 0, 0];
- * const unitDirection = VectorNormalize(direction); // [1, 0, 0]
+     * ```typescript
+     * const vector = [3, 4, 0];
+     * const normalized = VectorNormalize(vector); // [0.6, 0.8, 0] (magnitude = 1)
+     * const direction = [10, 0, 0];
+     * const unitDirection = VectorNormalize(direction); // [1, 0, 0]
+     * ```
  */
 export function VectorNormalize(a) {
     AssertVector(a);
@@ -282,18 +296,19 @@ export function VectorNormalize(a) {
  * @returns The magnitude (length) of the vector
  *
  * @example
- * const velocity = [3, 4, 0];
- * const speed = VectorMagnitude(velocity); // 5.0
- *
- * const unitVector = [1, 0, 0];
- * const unitLength = VectorMagnitude(unitVector); // 1.0
+     * ```typescript
+     * const velocity = [3, 4, 0];
+     * const speed = VectorMagnitude(velocity); // 5.0
+     * const unitVector = [1, 0, 0];
+     * const unitLength = VectorMagnitude(unitVector); // 1.0
+     * ```
  */
 export function VectorMagnitude(a) {
     AssertVector(a);
     let sum = 0;
     for (const av of a) {
         AssertVectorValue(av, {});
-        sum += Math.pow(av, 2);
+        sum += av * av;
     }
     return Math.sqrt(sum);
 }
@@ -306,11 +321,12 @@ export function VectorMagnitude(a) {
  * @returns New vector with absolute values of all components
  *
  * @example
- * const vector = [-3, 4, -2];
- * const absolute = VectorAbs(vector); // [3, 4, 2]
- *
- * const mixed = [1.5, -2.7, 0];
- * const absValues = VectorAbs(mixed); // [1.5, 2.7, 0]
+     * ```typescript
+     * const vector = [-3, 4, -2];
+     * const absolute = VectorAbs(vector); // [3, 4, 2]
+     * const mixed = [1.5, -2.7, 0];
+     * const absValues = VectorAbs(mixed); // [1.5, 2.7, 0]
+     * ```
  */
 export function VectorAbs(a) {
     AssertVector(a);
@@ -329,11 +345,12 @@ export function VectorAbs(a) {
  * @returns True if all components are zero, false otherwise
  *
  * @example
- * const zero = [0, 0, 0];
- * const isZero = VectorIsZero(zero); // true
- *
- * const notZero = [0, 0.001, 0];
- * const isNotZero = VectorIsZero(notZero); // false
+     * ```typescript
+     * const zero = [0, 0, 0];
+     * const isZero = VectorIsZero(zero); // true
+     * const notZero = [0, 0.001, 0];
+     * const isNotZero = VectorIsZero(notZero); // false
+     * ```
  */
 export function VectorIsZero(vector) {
     return vector.every((v) => v === 0);
@@ -349,13 +366,14 @@ export function VectorIsZero(vector) {
  * @throws {VectorError} If either vector is zero
  *
  * @example
- * const right = [1, 0, 0];
- * const up = [0, 1, 0];
- * const angle = VectorAngle(right, up); // π/2 (90 degrees)
- *
- * const forward = [0, 0, 1];
- * const backward = [0, 0, -1];
- * const oppositeAngle = VectorAngle(forward, backward); // π (180 degrees)
+     * ```typescript
+     * const right = [1, 0, 0];
+     * const up = [0, 1, 0];
+     * const angle = VectorAngle(right, up); // π/2 (90 degrees)
+     * const forward = [0, 0, 1];
+     * const backward = [0, 0, -1];
+     * const oppositeAngle = VectorAngle(forward, backward); // π (180 degrees)
+     * ```
  */
 export function VectorAngle(a, b) {
     AssertVectors([a, b]);
@@ -375,9 +393,11 @@ export function VectorAngle(a, b) {
  * @returns New rotated 2D vector
  *
  * @example
- * const right = [1, 0];
- * const rotated90 = Vector2Rotate(right, Math.PI / 2); // [0, 1] (up)
- * const rotated180 = Vector2Rotate(right, Math.PI); // [-1, 0] (left)
+     * ```typescript
+     * const right = [1, 0];
+     * const rotated90 = Vector2Rotate(right, Math.PI / 2); // [0, 1] (up)
+     * const rotated180 = Vector2Rotate(right, Math.PI); // [-1, 0] (left)
+     * ```
  */
 export function Vector2Rotate(vector, radians) {
     AssertVector2(vector);
@@ -397,9 +417,11 @@ export function Vector2Rotate(vector, radians) {
  * @returns Unit vector pointing in the specified direction
  *
  * @example
- * const right = Vector2FromAngle(0); // [1, 0]
- * const up = Vector2FromAngle(Math.PI / 2); // [0, 1]
- * const diagonal = Vector2FromAngle(Math.PI / 4); // [0.707, 0.707]
+     * ```typescript
+     * const right = Vector2FromAngle(0); // [1, 0]
+     * const up = Vector2FromAngle(Math.PI / 2); // [0, 1]
+     * const diagonal = Vector2FromAngle(Math.PI / 4); // [0.707, 0.707]
+     * ```
  */
 export function Vector2FromAngle(radians) {
     return [Math.cos(radians), Math.sin(radians)];
@@ -414,10 +436,12 @@ export function Vector2FromAngle(radians) {
  * @returns Scalar cross product (positive = counterclockwise, negative = clockwise)
  *
  * @example
- * const right = [1, 0];
- * const up = [0, 1];
- * const cross = Vector2Cross(right, up); // 1 (counterclockwise)
- * const crossReverse = Vector2Cross(up, right); // -1 (clockwise)
+     * ```typescript
+     * const right = [1, 0];
+     * const up = [0, 1];
+     * const cross = Vector2Cross(right, up); // 1 (counterclockwise)
+     * const crossReverse = Vector2Cross(up, right); // -1 (clockwise)
+     * ```
  */
 export function Vector2Cross(a, b) {
     AssertVector2(a);
@@ -435,9 +459,11 @@ export function Vector2Cross(a, b) {
  * @throws {VectorError} If vector b is zero
  *
  * @example
- * const force = [5, 3, 0];
- * const surface = [1, 0, 0];
- * const perpendicular = Vector3Reject(force, surface); // [0, 3, 0]
+     * ```typescript
+     * const force = [5, 3, 0];
+     * const surface = [1, 0, 0];
+     * const perpendicular = Vector3Reject(force, surface); // [0, 3, 0]
+     * ```
  */
 export function Vector3Reject(a, b) {
     AssertVector3(a);
@@ -459,16 +485,18 @@ export function Vector3Reject(a, b) {
  * @throws {VectorError} If vector b is zero
  *
  * @example
- * const force = [5, 3, 0];
- * const surface = [1, 0, 0];
- * const parallel = VectorProject(force, surface); // [5, 0, 0]
+     * ```typescript
+     * const force = [5, 3, 0];
+     * const surface = [1, 0, 0];
+     * const parallel = VectorProject(force, surface); // [5, 0, 0]
+     * ```
  */
 export function VectorProject(a, b) {
     AssertVectors([a, b]);
     if (VectorIsZero(b))
         throw new VectorError('Cannot project onto a zero vector');
     const dot = VectorDot(a, b);
-    const magSquared = Math.pow(VectorMagnitude(b), 2);
+    const magSquared = VectorDot(b, b);
     const scalar = dot / magSquared;
     const result = [];
     for (const bv of b) {
@@ -488,9 +516,11 @@ export function VectorProject(a, b) {
  * @throws {VectorError} If the normal is a zero vector
  *
  * @example
- * const incoming = [1, -1, 0];
- * const normal = [0, 1, 0]; // surface normal (upward)
- * const reflected = Vector3Reflect(incoming, normal); // [1, 1, 0]
+     * ```typescript
+     * const incoming = [1, -1, 0];
+     * const normal = [0, 1, 0]; // surface normal (upward)
+     * const reflected = Vector3Reflect(incoming, normal); // [1, 1, 0]
+     * ```
  */
 export function Vector3Reflect(incident, normal) {
     AssertVector3(incident);
@@ -510,11 +540,12 @@ export function Vector3Reflect(incident, normal) {
  * @returns Vector perpendicular to both a and b (following right-hand rule)
  *
  * @example
- * const right = [1, 0, 0];
- * const forward = [0, 0, 1];
- * const up = Vector3Cross(right, forward); // [0, 1, 0]
- *
- * const normal = Vector3Cross([1, 0, 0], [0, 1, 0]); // [0, 0, 1]
+     * ```typescript
+     * const right = [1, 0, 0];
+     * const forward = [0, 0, 1];
+     * const up = Vector3Cross(right, forward); // [0, 1, 0]
+     * const normal = Vector3Cross([1, 0, 0], [0, 1, 0]); // [0, 0, 1]
+     * ```
  */
 export function Vector3Cross(a, b) {
     AssertVector3(a);
@@ -535,9 +566,11 @@ export function Vector3Cross(a, b) {
  * @returns Magnitude of the cross product
  *
  * @example
- * const side1 = [3, 0, 0];
- * const side2 = [0, 4, 0];
- * const area = VectorCrossMagnitude(side1, side2); // 12 (area of rectangle)
+     * ```typescript
+     * const side1 = [3, 0, 0];
+     * const side2 = [0, 4, 0];
+     * const area = VectorCrossMagnitude(side1, side2); // 12 (area of rectangle)
+     * ```
  */
 export function VectorCrossMagnitude(a, b) {
     AssertVector3(a);
@@ -556,10 +589,12 @@ export function VectorCrossMagnitude(a, b) {
  * @returns Signed volume (positive = right-handed orientation)
  *
  * @example
- * const x = [1, 0, 0];
- * const y = [0, 1, 0];
- * const z = [0, 0, 1];
- * const volume = Vector3ScalarTripleProduct(x, y, z); // 1 (unit cube)
+     * ```typescript
+     * const x = [1, 0, 0];
+     * const y = [0, 1, 0];
+     * const z = [0, 0, 1];
+     * const volume = Vector3ScalarTripleProduct(x, y, z); // 1 (unit cube)
+     * ```
  */
 export function Vector3ScalarTripleProduct(a, b, c) {
     AssertVector3(a);
@@ -579,10 +614,12 @@ export function Vector3ScalarTripleProduct(a, b, c) {
  * @returns Vector result of a × (b × c)
  *
  * @example
- * const a = [1, 0, 0];
- * const b = [0, 1, 0];
- * const c = [0, 0, 1];
- * const result = Vector3TripleProduct(a, b, c); // [0, 0, 0]
+     * ```typescript
+     * const a = [1, 0, 0];
+     * const b = [0, 1, 0];
+     * const c = [0, 0, 1];
+     * const result = Vector3TripleProduct(a, b, c); // [0, 0, 0]
+     * ```
  */
 export function Vector3TripleProduct(a, b, c) {
     AssertVector3(a);
@@ -602,9 +639,11 @@ export function Vector3TripleProduct(a, b, c) {
  * @returns The reflected vector
  *
  * @example
- * const incoming = [1, -1, 0];
- * const wall = [0, 1, 0]; // vertical wall normal
- * const bounced = VectorReflect(incoming, wall); // [1, 1, 0]
+     * ```typescript
+     * const incoming = [1, -1, 0];
+     * const wall = [0, 1, 0]; // vertical wall normal
+     * const bounced = VectorReflect(incoming, wall); // [1, 1, 0]
+     * ```
  */
 export function VectorReflect(incident, normal) {
     AssertVectors([incident, normal]);
@@ -630,11 +669,12 @@ export function VectorReflect(incident, normal) {
  * @returns Vector with all components negated
  *
  * @example
- * const forward = [0, 0, 1];
- * const backward = VectorNegate(forward); // [0, 0, -1]
- *
- * const velocity = [5, -3, 2];
- * const opposite = VectorNegate(velocity); // [-5, 3, -2]
+     * ```typescript
+     * const forward = [0, 0, 1];
+     * const backward = VectorNegate(forward); // [0, 0, -1]
+     * const velocity = [5, -3, 2];
+     * const opposite = VectorNegate(velocity); // [-5, 3, -2]
+     * ```
  */
 export function VectorNegate(a) {
     const result = [];
@@ -661,10 +701,12 @@ export function VectorNegate(a) {
  * @throws {VectorError} If any divisor component is zero
  *
  * @example
- * const velocity = [20, 10, 0];
- * const halved = VectorDivide(velocity, 2); // [10, 5, 0] - scalar division
- * const factors = [2, 5, 1];
- * const componentWise = VectorDivide(velocity, factors); // [10, 2, 0] - component-wise
+     * ```typescript
+     * const velocity = [20, 10, 0];
+     * const halved = VectorDivide(velocity, 2); // [10, 5, 0] - scalar division
+     * const factors = [2, 5, 1];
+     * const componentWise = VectorDivide(velocity, factors); // [10, 2, 0] - component-wise
+     * ```
  */
 export function VectorDivide(a, b) {
     const result = [];
@@ -705,12 +747,13 @@ export function VectorDivide(a, b) {
  * @returns New vector with each component clamped between min and max
  *
  * @example
- * const v = [5, -3, 12, 0];
- * VectorClamp(v, 0, 10); // [5, 0, 10, 0] - scalar bounds
- *
- * const mins = [0, -5, 0, -1];
- * const maxs = [10, 5, 8, 1];
- * VectorClamp(v, mins, maxs); // [5, -3, 8, 0] - per-component bounds
+     * ```typescript
+     * const v = [5, -3, 12, 0];
+     * VectorClamp(v, 0, 10); // [5, 0, 10, 0] - scalar bounds
+     * const mins = [0, -5, 0, -1];
+     * const maxs = [10, 5, 8, 1];
+     * VectorClamp(v, mins, maxs); // [5, -3, 8, 0] - per-component bounds
+     * ```
  */
 export function VectorClamp(a, min, max) {
     AssertVector(a);
@@ -736,44 +779,58 @@ export function VectorClamp(a, min, max) {
  * @throws {VectorError} If max is negative
  *
  * @example
- * const velocity = [15, 20, 0]; // magnitude ≈ 25
- * const limited = VectorLimit(velocity, 10); // magnitude = 10, same direction
- *
- * const small = [1, 1, 0]; // magnitude ≈ 1.414
- * const unchanged = VectorLimit(small, 5); // unchanged since already under limit
- */
-export function VectorLimit(a, max) {
-    if (max < 0)
-        throw new VectorError('Maximum magnitude cannot be negative');
-    const magnitude = VectorMagnitude(a);
-    if (magnitude <= max || magnitude === 0)
-        return VectorClone(a);
-    const scaleFactor = max / magnitude;
-    return VectorMultiply(a, scaleFactor);
-}
-/**
- * Validates if the input is a properly formatted vector.
- * Performs comprehensive validation without throwing errors.
- * Useful for input validation and defensive programming.
- *
- * @param vector - Input to validate
- * @returns True if input is a valid vector, false otherwise
- *
- * @example
- * const valid = VectorIsValid([1, 2, 3]); // true
- * const invalid = VectorIsValid("not a vector"); // false
- * const nullVector = VectorIsValid(null); // false
- * const emptyArray = VectorIsValid([]); // depends on implementation
- */
-export function VectorIsValid(vector) {
-    try {
+     * ```typescript
+     * const velocity = [15, 20, 0]; // magnitude ≈ 25
+     * const limited = VectorLimit(velocity, 10); // magnitude = 10, same direction
+     * const small = [1, 1, 0]; // magnitude ≈ 1.414
+     * const unchanged = VectorLimit(small, 5); // unchanged since already under limit
+     * ```
+export function VectorGramSchmidt<T extends TAnyVector>(vectors: T[], normalize: boolean = false): TVectorResult<T>[] {
+    if (vectors.length === 0) throw new VectorError('GramSchmidt: Empty Vector Set');
+    const [firstVector] = vectors;
+    if (!firstVector) throw new VectorError('GramSchmidt: Undefined First Vector');
+    const dimension = firstVector.length;
+    for (const [i, vector] of vectors.entries()) {
         AssertVector(vector);
-        return true;
+        if (vector.length !== dimension) throw new VectorError(`GramSchmidt: Vector at index ${i} has different dimension than first vector. Expected ${dimension}, got ${vector.length}`);
+        if (VectorIsZero(vector)) throw new VectorError(`GramSchmidt: Vector at index ${i} is a zero vector. Cannot orthogonalize zero vectors.`);
     }
-    catch {
-        return false;
+
+    const result: TVectorResult<T>[] = [];
+
+    for (const [i, currentVector] of vectors.entries()) {
+        AssertVector(currentVector);
+
+        // Start with a clone of the current vector
+        let orthogonalVector = VectorClone(currentVector);
+
+        // Subtract projections of all previous vectors
+        for (let j = 0; j < i; j++) {
+            const previousVector = result[j];
+            AssertVector(previousVector);
+
+            // Both vectors have the same dimension, projection result is compatible
+            const vectorT = currentVector as unknown as T;
+            const prevT = previousVector as unknown as T;
+            const projection = VectorProject(vectorT, prevT);
+            orthogonalVector = VectorSubtract(orthogonalVector as unknown as T, projection as unknown as T);
+        }
+        
+        if (VectorIsZero(orthogonalVector as unknown as T)) {
+            throw new VectorError(`GramSchmidt: Vector at index ${i} is linearly dependent on previous vectors. Cannot orthogonalize linearly dependent vectors.`);
+        }
+        
+        if (normalize) {
+            orthogonalVector = VectorNormalize(orthogonalVector as unknown as T);
+        }
+        
+        result.push(orthogonalVector as TVectorResult<T>);
     }
+
+    return result;
+}}
 }
+
 /**
  * Performs Gram-Schmidt orthogonalization on a set of vectors.
  * Converts a set of linearly independent vectors into an orthogonal (or orthonormal) set.
@@ -786,9 +843,11 @@ export function VectorIsValid(vector) {
  * @throws {VectorError} If vectors are linearly dependent or invalid
  *
  * @example
- * const vectors = [[1, 1, 0], [1, 0, 1], [0, 1, 1]];
- * const orthogonal = VectorGramSchmidt(vectors); // Orthogonal set
- * const orthonormal = VectorGramSchmidt(vectors, true); // Orthonormal set
+     * ```typescript
+     * const vectors = [[1, 1, 0], [1, 0, 1], [0, 1, 1]];
+     * const orthogonal = VectorGramSchmidt(vectors); // Orthogonal set
+     * const orthonormal = VectorGramSchmidt(vectors, true); // Orthonormal set
+     * ```
  */
 export function VectorGramSchmidt(vectors, normalize = false) {
     if (vectors.length === 0)