关于一道习题的数值解法及其讨论.pdfVIP

第 28卷第3期 大 学 物 理 Vo1．28N0．3 2009年 3月 C0LLEGE PHYSICS Mar．2009 关于一道习题 的数值解法及其讨论 冯育俊 (云南大学 材料科学与工程系 ，云南 昆明 650091) 摘要：利用 自编 c语言程序的模拟运算 ，指出了 《近代物理》(杨桂林，江兴方 ，柯善哲主编——科学出版社，2004年，21世 纪高等院校教材)中第三章 “薛定谔方程”习题 13的答案并不正确 ，通过数值计算给 出了正确 的结果，并对该 问题进行 了较为 深入的探讨 ． 关键词 ：二维无限深势阱；简并度 ；量子态 ；C语言编程 ；不定方程 中图分类号 ：0413．1 文献标识码 ：A 文章编号 ：1000—0712(2009)03—0052—03 在二维无 限深势阱 中，粒子 能量 E ：En= d=0： E(n+n：)，E为基态能量，量子数n，n分别是不 f0r(X=1；X=sqrt(n)；X++) 一 定相同的整数组合． f0r(Y=1；Y=sqrt(n)；Y++) (原问题)对于二维盒 中的粒子 ，能量在 E= if(x x+Y Y=n) 50E．之上 ，简并度 比2大 的最邻近 的能级是 哪 { 个?… if(13==X X+Y y) 教材给出的能级n是200，简并度为3，量子态 ： { (n，n)：(2，14)、(14，2)和 (10，10)．其实 ，这个答 d++ ： 案是不符合题意的． a[d]=X； b[d]=y； 1 数值解法给 出正确结果 } 根据题意 ，假设各量子态 ：n，n在 c程序 中的 } 代码为 ，y；能级仍然为 7／,，即n= +Y；该能级的 elsebreak； 简并度为D，输 出时为 d，冒号后边为量子态 ，编写 if(d2) 了如下 C程序 ： { # ineludestdio．h printf(n=％d，d=％d： ，n，d)； # includemath．h for(i=1；i=d；i++) #defineZ50 { voidmain() printf(”(％d，”，a[i])； { printf(”％d)”，b[i])； ints，m，n，x，y，d，i，a[z]，b[z]； } printf(”Input[s，m]=”)； printf(＼n)； S~Janf(”[％d，％d])”，s，＆m)；