《Fortran程序设计》实验

实验一实验内容将以下的代码输入并编译验证实验目的掌握下编程的方法linux Fortran注意事项源代码命名为存放目录在自己的主目录下的文件夹paLf90,Pl在修改时间与作者后面填写相关信息，作者填写学号与姓名PROGRAM roots!目的!最后修改时间:作者:!本程序是对形式为a*x**2+b*x+c=0的二次方程的根进行求解!数据字典:声明变量类型、定义和单位REAL!二次项系数REAL!一次项系数REAL!常数项REAL:discriminant!方程的判别式REAL:imag_part!复根的虚部REAL:real_part!复根的实部REAL:xl!方程的第一个实根REAL:x2!方程的第二个实根IMPLICIT NONE!提示用户输入方程式的系数WRITE*,*This programWRITE*,**equation ofsolves forthe rootsof aquadraticthe formA*X**2+B*X+C=OWRITE*,*J Enterthe coefficientsA,B,and C:READ*,*a,b,c!回显系数WRITE*,*The coefficientsA,B,and Care/,a,b,c!计算判别式discriminant=b**2—

4.*a*cIF discriminant0THEN!如果判别式大于3有两个不等的实根xl=-b+sqrt discriminant/

2.*ax2二-b-sqrtdiscriminant/

2.*aWRITE*,*This equation has tworeal roots/WRITE*,*Xl=,xlELSE IFdiscriminant0THEN!如果判别式小于3有两个虚根WRITE*,*X2=,x2realpart~b/

2.*aimagpart sqrtabsdiscriminant/

2.*aWRITE*,*J This equationhascomplex rootsWRITE*,*Xl-,real_part,J+i,imag partWRITE*,*X2-,real part,*-i,imag partELSE!如果判别式等于0,有两个相等的实根xl=-b/

2.*aWRITE*,*Thisequationhas twoidentical realroots:WRITE*,*，X1=X2=,xl ENDIF实验二实验内容复合梯形求END PROGRAMroots积分法的公式如下Ja271-1将积分区域分为n个相等的/a+/S+2Z/〃+//=1子空间，对每个子空间近似进行梯形求积式中h=b-a/n,也就是每个小梯形的宽度编程利用这个公式，计算积分，+sinx公£分别把区间分为40,80和200个小区间，比较计算精度的情况实验目的掌握控制语句、循环语句的用法注意事项源代码命名为存放目录在自己的主目录下的文件夹在修改时pa.f90,P2间与作者后面填写相关信息，作者填写学号与姓名实验文档提交截止时间2015年3月29日24:00实验三实验内容矩阵的乘法矩阵和相乘，结果放在中AN,L BL,M CN,MC的元素由下式计算刀附葭cik=£j=T编程计算c设117-24A=—7143B=2-22-35注意事项源代码命名为pa.f90,存放目录在自己的主目录下的P3文件夹实验文档提交截止时间2015年4月5日24:00,不要在这个时间之后再修改你的文档实验三实验内容在三维空间中的任意一点P的位置可以表示为三个数值的集合x,y,z,两点间的距离可以用下列公式计算d=7xl—x22+yl—y22+zl—z22编程，从键盘读入两点的坐标值，计算它们之间的距离并输出结果，显示结果用ES格式保留4位有效数值注意事项源代码命名为pa.f90,存放目录在自己的主目录下的P4文件夹实验文档提交截止时间2015年4月12日24:00,不要在这个时间之后再修改你的文档实验三实验内容:金属盘温度分布:在稳定条件下,一个金属盘表面的任一点的温度是它周围所有点的温度的平均值这一特性可用于重复地计算，求出分布在一个盘子上的所有点的温度现有一个方形的盘子，被网格分成100个小方块或节点10行10列节点的温度形成了一个二维数组T在盘子的边缘的所有节点的温度被一个冷却系统控制为20度，节点3,8的温度通过接触沸腾的水使其固定在100度任意给定节点的温度值可以通过围绕它的所有方块的温度的平均值计算得到:+工+虱Tij,new=I Ti+i,j+TJI jj+i j_i为了求得分布在盘子表面的温度,必须对每个节点的温度做一个初始的设定，然后用以上的公式估算每一个温度不固定的节点的温度值,这些更新过的温度值又可以用于计算更新的估算值，重复这个过程,直到每个节点的估算值和旧的值都仅有一点点不同在那点上，将会找到一个比较稳定的结果编写程序，计算整个盘子上的稳定温度值分布，假设所有内部节点的初始值都被设为50度记住所有外部节点都稳定在20度，节点3,8的温度固定在100度程序应该反复应用计算公式，直到任意节点的前后变化小于

0.01度最后按照10x10的格式输出整个注盘的温度值意事项源代码命名为pa.f90,存放目录在自己的主目录下的P5文件夹实验文档提交截止时间2015年4月19日24:00,不要在这个时间之后再修改你的文档实验三实验内容:编写程序，打开一个文件，将1-10的整数写入到前10个记录中然后，在文件中向回移动6个记录，读取存储在该记录中的数值将这个数值保存在变量x中然后再向回移动3个记录，读取存储在该记录中的数值将该数值保存在变量y中，将x与y相乘，输出结果注意事项源代码命名为pa.f90,存放目录在自己的主目录下的P6文件夹实验文档提交截止时间2015年5月2日24:00,不要在这个时间之后再修改你的文档。