fortran例题6

工程分析程序设计上机作业（六）

上机目的：练习派生类型的定义、构造、初始化，成员的调用和操作，操作符的重载，指针

1、构造一个由年、月、日组成的派生类型，并编写程序，实现从键盘上输入年、月、日，由程序计算出

是当年度的第几天。例如：

输入：2011, 4, 8

输出：4月8日是2011年的第98天。

注意：需考虑闰年的情形。

program main

type day !day 是年、月、日组成的派生类型

integer year

integer month

integer date

end type

type(day)::d

integer i,j,k,n,m ! n 为天数

n=0

print*,"please imput the year month date"

read*,d%year,d%month,d%date

!判断是否为闰年

if(mod(d%year,4)==0.and.mod(d%year,100)/=0.or.mod(d%year,400)==0) then

do i=1,d%month-1 !31天的月1 3 5 7 8 10 12，30天4 6 9 11 ，29天2

if(i==1.or.i==3.or.i==5.or.i==7.or.i==8.or.i==10.or.i==12) then

n=n+31

else if(i==2) then

n=n+29

else

n=n+30

end if

end do

else !不是闰年

do i=1,d%month-1 !31天的月1 3 5 7 8 10 12，30天4 6 9 11 ，28天2

if(i==1.or.i==3.or.i==5.or.i==7.or.i==8.or.i==10.or.i==12) then

n=n+31

else if(i==2) then

n=n+28

else

n=n+30

end if

end do

end if

n=n+d%date

write(*,"(I5,'年',I3,'月',I3,'日是第',I3,'天')",advance='no') d%year,d%month,d%date,n

end

2、设计一个“-”操作符的重载，把集合A中那些同时又出现在集合B中的元素去掉。

module chaji

implicit none

real,dimension:: A,B,C,D

interface operator(-)

MODULE PROCEDURE cha

end interface

contains

subroutine cha(A,B,C,D)

! program main !从A、B两个数列中，把同时出现在两个数列中的数据删去。例如：

do i=1,n

do j=1,n

if(A(i)==A(j).and.i/=j) then

A(i)=' '

end if

end do

end do

do i=1,m

do j=1,m

if(B(i)==B(j).and.i/=j) then

B(i)=' '

end if

end do

end do

do i=1,n

do j=1,m

if(A(i)==B(j)) then

A(i)=' '

B(j)=" "

end if

end do

end do

C=A

D=B

j=1

do i=1,n

if(a(i)/=' ') then

C(j)=A(i)

j=j+1

end if

end do

i=1

do J=1,m

if(B(j)/=" ") then

C(i)=D(j)

i=i+1

end if

end do

end subroutine

program main !从A 、B 两个数列中，把同时出现在两个数列中的数据删去。例如

use chaji

implicit none

integer::i,j,n,m

integer,allocatable::A(:),B(:),C(:),D(:)

print *,'how many in A,B: '

read*,n ,m

allocate (A(n))

allocate(C(n))

allocate(d(m))

allocate (B(m))

print*,"please input A"

do i=1,n

read*,A(i)

end do

print*,"input B:"

do j=1,m

read*,B(j)

end do

deallocate(A)

deallocate(B)

deallocate(C)

deallocate(D)

end

3、 编写一程序用四阶龙格-库塔法求解微分方程2*d y

y x y d x =-，当x=0时，y=1.0。试求出x=0.1, 0.2,

0.3, 0.4…., 1.0,…, 100时的y 值。算法如下：求解))(,('x y x f y =，定解条件：x =x 0，y =y 0。

()()()1123412132431226,1,221,22,n n n n n n n n n n y y k k k k k h f x y h k h f x y k h k h f x y k k h f x h y k +?=++++??

?=????=++? ????

???=++? ????

?=++?

已知n y y y x f y y x f y y x x ??????????==21110000),(,),(。

要求：1）把结果按 x 1, y 1

x 2, y 2

……

……

的排列顺序输出到一个无格式文档data1.txt中，一个有格式文档data2.txt。注意比较两个文档占用存储空间的大小。

2）编写另一个程序，从上述文档中提取出数据，然后找出这个函数的最大值y max，和最大值所在的x max，以及函数的最小值y min，和最小值所在的x min，从屏幕上输出。寻找最大值和最小值的算法要求用指针完成：设定一个最大值指针y max和x max，开始指向第一个y值和x值。然后提取下一个y值与指针指向的数值进行比较，如果更大，则令指针指向新的y值和对应的x值，如此扫描所有数据，指针指向的最后的数值，就是所求的数值。

选作：

4、有n个人围成一圈，从第一个人开始报数（1-3报数），凡报到3退出圈子，问最后留下的是原来第几

号位的，要求用单链表完成。