第六篇 函数.pdfVIP

1998年诺贝尔化学奖得主库恩（Kohn） 计算物理和量子化学里著名的密度泛函理论 核心方程Kohn-Sham方程 18岁时 89岁时 人都会老去，要趁年轻多干点事情，不是吗？ 问题 1. 例子ex0714.f90 a(1,1)=2 a(1,2)=3 a(1,3)=3 a(1,4)=3 a(1,5)=3 a(2,1)=1 a(2,2)=2 a(2,3)=3 a(2,4)=3 a(2,5)=3 a(3,1)=1 a(3,2)=1 a(3,3)=2 a(3,4)=3 a(3,5)=3 a(4,1)=1 a(4,2)=1 a(4,3)=1 a(4,4)=2 a(4,5)=3 a(5,1)=1 a(5,2)=1 a(5,3)=1 a(5,4)=1 a(1,5)=2 write(*,(5(5I5,/))) a a(1,1)=2 a(2,1)=1 a(3,1)=1 a(4,1)=1 a(5,1)=1 a(1,2)=3 a(2,2)=2 a(3,2)=1 a(4,2)=1 a(5,2)=1 问题 2. 例子ex0719.f90 real :: A(L,M) = (/ 1,2,3,4,5,6,7,8,9,10,11,12 /) do i=1,L write(*,*) A(i,:) end do write(*,*) A(i,:) A(1,1)=1 a(1,2)=4 a(1,3)=7 a(1,4)=10 第六章 函数 一、子程序（subroutine ）的使用 写程序时，可以把某一段常常被使用、具有特定功能的程序代码独 立出来，封装成子程序，以后只要经过调用的CALL命令就可以执行这一 段程序代码。先来看一个使用子程序的实例。 ex0801.f90 程序的第3、4行出现了一个新的命令“CALL”，CALL在程序中的意 思就是“调用”的意思。 子程序中也包含可执行的程序代码，这段程序代码的内容可以是任 何的Fortran命令描述。在这里面也可以声明变量、做流程控制、执行循 环、甚至调用其他子程序。 一、子程序（subroutine ）的使用 子程序和第4章一开始所提到过的 “主程序” （以program开头，end 来结束之间的这一段程序代码）最大不同处在于： “主程序的程序代码， 在程序一开始就自动会被执行，而子程序则不会自动执行，它需要被别人 ‘调用’才会执行。” 子程序的程序代码以subroutine开头，它同样要取一个名字，以end或 end subroutine来结束。严格地讲，end subroutine后还可以再接上这个子 程序的名字。 子程序的名字比起主程序来重要多了，因为主程序的名字可以随时更 改，而且在改名时只需更动程序中program后的那几个字就行了，而子程 序的名字是用来给别人调用时所要使用的，在取名时最好能配合这个子程 序的功能来给予适当的称呼（和变量一样）。改变子程序的名字时，还要 把在程序中所有CALL它旧名字的部分全部改成CALL它新的名字。 一、子程序（subroutine ）的使用 一个包含子程序的Fortran程序在结构上的模样大致如下： 一、子程序（subroutine ）的使用 主程序并不一定要放在程序代码的最开头，它可以安排在程序中的任意 位置，可以先写子程序再写主程序，以下的结构也是合乎标准的写法。 一、子程序（subroutine