From 6e77c81997a622b11687a50df230e1cd9a88caaa Mon Sep 17 00:00:00 2001
From: Michael McGuffin <47359290+MJMcGuffin@users.noreply.github.com>
Date: Wed, 10 Jul 2024 14:21:41 -0400
Subject: [PATCH] new version of mathjs, and added comments about concurrence
 and von neumann entropy

---
 index.html | 7 ++++++-
 1 file changed, 6 insertions(+), 1 deletion(-)

diff --git a/index.html b/index.html
index e258581..b0e4d68 100644
--- a/index.html
+++ b/index.html
@@ -7,7 +7,7 @@
 Or go to the <a href="https://github.com/MJMcGuffin/muqcs">github repository</a> page.
 </p>
 
-<script src="https://cdnjs.cloudflare.com/ajax/libs/mathjs/12.3.2/math.js"></script>
+<script src="https://cdnjs.cloudflare.com/ajax/libs/mathjs/13.0.2/math.js"></script>
 
 <script>
 
@@ -1402,6 +1402,8 @@
     //    https://docs.quantum.ibm.com/api/qiskit/0.24/qiskit.quantum_info.concurrence
     //    https://github.com/Qiskit/qiskit/blob/stable/0.16/qiskit/quantum_info/states/measures.py
     //       (lines 208-213 at the end of the concurrence() routine)
+    // More discussion is at
+    //    https://physics.stackexchange.com/questions/46443/what-is-the-motivation-for-the-definition-of-concurrence-in-quantum-information
     static computePairwiseQubitConcurrences(
         n, // number of qubits in the circuit
         densityMatrix // a matrix of size (2**n)x(2**n)
@@ -1475,6 +1477,9 @@
     // (In general, when considering n qubits, the von Neumann entropy computed from their 2^n x 2^n reduced density matrix
     // is in the range [0,n])
     //
+    // See
+    //     https://quantumcomputing.stackexchange.com/questions/28708/how-is-the-von-neumann-entropy-of-a-state-defined-from-its-eigendecomposition
+    //
     static computePairwiseQubitVonNeumannEntropy(
         n, // number of qubits in the circuit
         densityMatrix // a matrix of size (2**n)x(2**n)