אני מנסה לחשב את אלמנט המטריצה בין מצבי s1,s2 ל pz
בפתרון Y12 , עמוד 2, בשורה שבה מחושב הצימוד בין s1 לpz , למעשה נטען שהפתרון של האינטגרל הראדיאלי הוא c1, ולגבי הצימוד השני התוצאה היא c2. ע"פ חישוב אנליטי שלי + בדיקה בwolfram התוצאה היא לא c1:
integral_0^infinity sqrt(2) exp(-x) x^3 (1-x\/2) exp(-x\/2) dx = -(32 sqrt(2))/81 ~~ -0.55870
אנליטית זה מחושב ע"י
ואת זה צריך להכפיל ב
לגבי הצימוד השני:
integral_0^infinity (exp(-x) x^3 (1-x\/2) x)\/(sqrt(12) 2) dx = -3 sqrt(3) ~~ -5.1962
ושוב יש להכפיל ב
לסיכומו של דבר, אם אני צודק, אני ממליץ לשנות את ההגדרה של c1 כך שתתאם להיות האינטגרל שרשום ב(*).