Dear All, Could you please check what is wrong with my code?
hbar=1;​​m1=1.209;​​m2=1.209;​​μ=
m1*m2
m1+m2
;​​NSolvem1+m2+
6A
2
δ
+
3B
δ
-
2μ
2
hbar
2
C1+
8A
3
δ
+
3B
2
δ
3+
8μ
2
hbar
B
3
δ
+
3A
4
δ
-D1+
2
hbar
2μ
2
(0+0.5)
==3.097,m1+m2+
6A
2
δ
+
3B
δ
-
2μ
2
hbar
2
C1+
8A
3
δ
+
3B
2
δ
5+
8μ
2
hbar
B
3
δ
+
3A
4
δ
-D1+
2
hbar
2μ
2
(0+0.5)
==3.686,m1+m2+
6A
2
δ
+
3B
δ
-
2μ
2
hbar
2
C1+
8A
3
δ
+
3B
2
δ
7+
8μ
2
hbar
B
3
δ
+
3A
4
δ
-D1+
2
hbar
2μ
2
(0+0.5)
==4.039,m1+m2+
6A
2
δ
+
3B
δ
-
2μ
2
hbar
2
C1+
8A
3
δ
+
3B
2
δ
3+
8μ
2
hbar
B
3
δ
+
3A
4
δ
-D1+
2
hbar
2μ
2
(1+0.5)
3.511,m1+m2+
6A
2
δ
+
3B
δ
-
2μ
2
hbar
2
C1+
8A
3
δ
+
3B
2
δ
5+
8μ
2
hbar
B
3
δ
+
3A
4
δ
-D1+
2
hbar
2μ
2
(1+0.5)
3.927,m1+m2+
6A
2
δ
+
3B
δ
-
2μ
2
hbar
2
C1+
8A
3
δ
+
3B
2
δ
3+
8μ
2
hbar
B
3
δ
+
3A
4
δ
-D1+
2
hbar
2μ
2
(2+0.5)
3.77,{A,B,C1,D1,δ},Reals;​​v=Totalm1+m2+
6A
2
δ
+
3B
δ
-
2μ
2
hbar
2
C1+
8A
3
δ
+
3B
2
δ
3+
8μ
2
hbar
B
3
δ
+
3A
4
δ
-D1+
2
hbar
2μ
2
(0+0.5)
==3.097,m1+m2+
6A
2
δ
+
3B
δ
-
2μ
2
hbar
2
C1+
8A
3
δ
+
3B
2
δ
5+
8μ
2
hbar
B
3
δ
+
3A
4
δ
-D1+
2
hbar
2μ
2
(0+0.5)
==3.686,m1+m2+
6A
2
δ
+
3B
δ
-
2μ
2
hbar
2
C1+
8A
3
δ
+
3B
2
δ
7+
8μ
2
hbar
B
3
δ
+
3A
4
δ
-D1+
2
hbar
2μ
2
(0+0.5)
==4.039,m1+m2+
6A
2
δ
+
3B
δ
-
2μ
2
hbar
2
C1+
8A
3
δ
+
3B
2
δ
3+
8μ
2
hbar
B
3
δ
+
3A
4
δ
-D1+
2
hbar
2μ
2
(1+0.5)
3.511,m1+m2+
6A
2
δ
+
3B
δ
-
2μ
2
hbar
2
C1+
8A
3
δ
+
3B
2
δ
5+
8μ
2
hbar
B
3
δ
+
3A
4
δ
-D1+
2
hbar
2μ
2
(1+0.5)
3.927,m1+m2+
6A
2
δ
+
3B
δ
-
2μ
2
hbar
2
C1+
8A
3
δ
+
3B
2
δ
3+
8μ
2
hbar
B
3
δ
+
3A
4
δ
-D1+
2
hbar
2μ
2
(2+0.5)
3.77/.Equal[l_,r_]Norm[l-(r)];​​NMinimize[v,{A,B,C1,D1,δ},MaxIterations
10
10
]​​​​
​