MATLAB_Comparison_4

fobos123deimos · Jul 23, 2024 · 56cce4e · 56cce4e
1 parent 7581d95
commit 56cce4e
Show file tree

Hide file tree

Showing 5 changed files with 36 additions and 0 deletions.
diff --git a/exemples/MATLAB_Comparison/MATLAB_Wavefunction_Comparison.mlx b/exemples/MATLAB_Comparison/MATLAB_Wavefunction_Comparison.mlx
diff --git a/exemples/MATLAB_Comparison/MATLAB_Wavefunction_Comparison.pdf b/exemples/MATLAB_Comparison/MATLAB_Wavefunction_Comparison.pdf
diff --git a/exemples/MATLAB_Comparison/MATLAB_Wavefunction_Precision_Tests.mlx b/exemples/MATLAB_Comparison/MATLAB_Wavefunction_Precision_Tests.mlx
diff --git a/exemples/MATLAB_Comparison/wavefunction_MATLAB_3.m b/exemples/MATLAB_Comparison/wavefunction_MATLAB_3.m
@@ -0,0 +1,18 @@
+function wavefunction = wavefunction_MATLAB_3(n,x,prec)
+
+    digits(prec);
+    x = vpa(x);
+    x_size = numel(x);
+
+    wavefunction = vpa(zeros(x_size, n+1));
+
+    for index = 1:n+1
+
+        norm = (2^(-0.5*vpa(index-1))) * (gamma(vpa(index-1)+1)^(-0.5)) * (pi^(-0.25));
+
+        H = hermiteH(vpa(index-1),x);
+
+        wavefunction(:,index) = norm * exp(-0.5 * x^2) * H;
+    end
+
+end
diff --git a/exemples/MATLAB_Comparison/wavefunction_MATLAB_4.m b/exemples/MATLAB_Comparison/wavefunction_MATLAB_4.m
@@ -0,0 +1,18 @@
+function wavefunction = wavefunction_MATLAB_4(n,x,prec)
+
+    digits(prec);
+    x = vpa(x);
+    x_size = numel(x);
+
+    wavefunction = vpa(zeros(n+1, x_size));
+
+    for index = 1:n+1
+
+        norm = (2^(-0.5*vpa(index-1))) * (gamma(vpa(index-1)+1)^(-0.5)) * (pi^(-0.25));
+
+        H = hermiteH(vpa(index-1),x);
+
+        wavefunction(index,:) = norm * exp(-0.5 * x.^2) .* H;
+    end
+
+end