原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
copputionalphysicshomework计算物理答案
Computational Physics姓名：王占雨导师：贾瑜 教授专业：凝聚态物理学号：201112131481第一题——1.1这实际是求解方程组的问题，将方程组写为行列式形式：0.125 0.015625 0.001953125 1.0 3.061467 0.0625 0 2-4 1.0 3.121445 0.03125 9.765625E-4 3.0517578125E-5 1.0 3.136548 0.015625 2-4 3.814697265625E-6 1.0 3.140331 通过逐项抵消，最终可得出PI的值，具体程序编写如下：publicclass Chapter_1 {/** * @param args */publicstaticvoid main(String[] args) {double [][] a=newdouble[4][5];int i,j,n;for(j=0,n=1;j=2;j++,n++) a[0][j]=1/Math.pow(8,n); a[0][3]=1; a[0][4]=3.061467;for(j=0,n=1;j=2;j++,n++) a[1][j]=1/Math.pow(16,n); a[1][3]=1; a[1][4]=3.121445;for(j=0,n=1;j=2;j++,n++) a[2][j]=1/Math.pow(32,n); a[2][3]=1; a[2][4]=3.136548;for(j=0,n=1;j=2;j++,n++) a[3][j]=1/Math.pow(64,n); a[3][3]=1; a[3][4]=3.140331;for(i=0;i=3;i++){for(j=0;j=4;j++){System.out.print(a[i][j]+ );if((j+1)%5==0)System.out.println(); }}for(i=0;i=2;i++) {double d=a[i][0]/a[3][0];for(j=0;j=4;j++){ a[i][j]=a[i][j]-a[3][j]*d; }}for(i=0;i=1;i++) {double d=a[i][1]/a[2][1];for(j=1;j=4;j++){ a[i][j]=a[i][j]-a[2][j]*d; } }for(j=2;j=4;j++) a[0][j]=a[0][j]-a[1][j]*a[0][2]/a[1][2];a[0][3]=(a[0][4]-a[1][4])/(a[0][3]-a[1][3]);for(i=0;i=3;i++){for(j=0;j=4;j++){System.out.print(a[i][j]+ );if((j+1)%5==0)System.out.println(); }}System.out.println(The value of PI is:+a[0][3]); }}运行程序后，可得出如下结果，求出PI的值为：3.14162255555555530.125 0.015625 0.001953125 1.0 3.061467 0.0625 0 2-4 1.0 3.121445 0.03125 9.765625E-4 3.0517578125E-5 1.0 3.136548 0.015625 2-4 3.814697265625E-6 1.0 3.140331 0.0 0.0 0.0 3.1416225555555553 65.974011 0.0 0.0 9.1552734375E-5 3.0 9.424805000000003 0.0 4.8828125E-4 2.288818359375E-5 -1.0 -3.1441140000000005 0.015625 2-4 3.814697265625E-6 1.0 3.140331 The value of PI is:3.1416225555555553第二题——1.2用欧拉方法给出在时刻的近似值为：用泰勒公式将y(t)在除展开可得：比较上面两式可知欧拉方法产生的误差为：我们
文档评论（0）