《Fortran程序设计》实验.docx

实验一

实验内容：将以下的代码输入并编译验证

实验目的：掌握linux下Fortran编程的方法

注意事项：源代码命名为pa1.f90，存放目录在自己的主目录下的P1文件夹

在修改时间与作者后面填写相关信息，作者填写学号与姓名

PROGRAMroots

!目的：

!本程序是对形式为a*x**2+b*x+c=0的二次方程的根进行求解。

!最后修改时间：作者：

IMPLICITNONE

!数据字典：声明变量类型、定义和单位

REAL::a!二次项系数

REAL::b!一次项系数

REAL::c!常数项

REAL::discriminant!方程的判别式

REAL::imag_part!复根的虚部

REAL::real_part !复根的实部

REAL::x1!方程的第一个实根

REAL::x2!方程的第二个实根

!提示用户输入方程式的系数

WRITE(*,*)Thisprogramsolvesfortherootsofaquadratic

WRITE(*,*)equationoftheformA*X**2+B*X+C=0

WRITE(*,*)EnterthecoefficientsA,B,andC:

READ(*,*)a,b,c

!回显系数

WRITE(*,*)ThecoefficientsA,B,andCare:,a,b,c

!计算判别式

discriminant=b**2-4.*a*c

IF(discriminant0)THEN!如果判别式大于0，有两个不等的实根

x1=(-b+sqrt(discriminant))/(2.*a)

x2=(-b-sqrt(discriminant))/(2.*a)

WRITE(*,*)Thisequationhastworealroots:

WRITE(*,*)X1=,x1

WRITE(*,*)X2=,x2

ELSEIF(discriminant0)THEN!如果判别式小于0，有两个虚根

real_part=(-b)/(2.*a)

imag_part=sqrt(abs(discriminant))/(2.*a)

WRITE(*,*)Thisequationhascomplexroots:

WRITE(*,*)X1=,real_part,+i,imag_part

WRITE(*,*)X2=,real_part,-i,imag_part

ELSE!如果判别式等于0，有两个相等的实根

x1=(-b)/(2.*a)

WRITE(*,*)Thisequationhastwoidenticalrealroots:

WRITE(*,*)X1=X2=,x1

ENDIF

ENDPROGRAMroots

实验二

实验内容：

复合梯形求积分法的公式如下：

将积分区域分为n个相等的子空间，对每个子空间近似进行梯形求积

式中h=(b-a)/n，也就是每个小梯形的宽度

编程利用这个公式，计算积分

分别把区间分为40,80和200个小区间，比较计算精度的情况

实验目的：掌握控制语句、循环语句的用法

注意事项：源代码命名为pa.f90，存放目录在自己的主目录下的P2文件夹

在修改时间与作者后面填写相关信息，作者填写学号与姓名

实验文档提交截止时间：2015年3月29日24:00

实验三

实验内容：

矩阵的乘法：

矩阵AN,L和BL,M相乘，结果放在CN,M中

C的元素由下式计算：

编程计算C

设

注意事项：源代码命名为pa.f90，存放目录在自己的主目录下的P3文件夹

实验文档提交截止时间：2015年4月5日24:00，不要在这个时间之后再修改你的文档

实验三

实验内容：

在三维空间中的任意一点P的位置可以表示为三个数值的集合（x,y,z），两点间的距离可以用下列公式计算：

d=

编程，从键盘读入两点的坐标值，计算它们之间的距离并输出结果，显示结果用ES格式保留4位有效数值。

注意事项：源代码命名为pa.f90，存放目录在自己的主目录下的P4文件夹

实验文档提交截止时间：2015年4月12日24:00，不要在这个时间之后再修改你的文档

实验三

实验内容：

金属盘温度分布：

在稳定条件下，一个金属