数组的使用及Fortran函数.pptVIP

苏州科技大学计算中心7.5语句函数/示例7.5语句函数[例7.3]计算空间一点P(x,y,z)的方向弦,其计算公式如下：其中:α、β、γ为点P到原点直线与坐标轴x、y、z的夹角。PROGRAMF751F(a,b,c,r)=r／SQRT(a*a+b*b+c*c)WRITE(*,*)输入x,y,z:READ(*,*)x,y,zDA＝F(x,y,z,x)DB＝F(x,y,z,y)DG＝F(x,y,z,z)WRITE(*,*)DA,DB,DGEND概述定义调用说明示例第七讲数组的使用及Fortran函数7.5语句函数苏州科技大学计算中心7.5语句函数/概述/示例在求解实际问题时,特别是在求解有关科学与工程数值计算问题时,常常会遇到计算一些数学表达式,如：这些数学表达式在程序中多次出现,需要重复书写它们,不但浪费了宝贵的时间和精力,而且也容易出错。任何数学表达式实际上都可表示成一元、二元或多元函数关系,即F(x1,x2,…,xn),这些函数关系如果没有对应的标准函数,则需要用户设法定义新的函数来实现它们。FORTRAN90提供了语句函数定义功能来定义一些简单的函数,用于实现有关的函数关系。用户自定义语句函数其使用方法与标准函数相同,可以作为操作数在表达式中调用。（示例）[例7.2]已知角度A＝36?43?27?,B＝15?24?,C＝8?16?54?,计算：S＝cos2A-4sinB?tgC角度Dn＝x1?x2?x3?的弧度转换公式:PROGRAMF91p1＝3.14159/180.0A=p1*(36.0+43.0/60+27.0/3600.0)B=p1*(15.0+24.0/60)C=p1*(8.0+16.0/60+54.0/3600)S=COS(A)**2-4.0*SIN(B)*TAN(C)WRITE(*,10)A,B,C,S10FORMAT(A=,F10.7,B=,F10.7,C=,F10.7,/S=,F10.7)END概述定义调用说明示例第七讲数组的使用及Fortran函数7.6外部函数子程序苏州科技大学计算中心7.6外部函数子程序/概述无法用一条语句来定义的函数，可以用一个单独的程序单位定义，它不属于调用函数语句所在的程序单位，因此称之为外部函数。是用户自定义的函数子程序。所在程序的其他程序单位都可以调用它。PROGRAMmainINTEGERaddPRINT*,add(5,3)ENDPROGRAMmainFUNCTIONadd(m,n)INTEGERm,n,addadd=m+nENDFUNCTIONaddFUNCTIONadd(m,n)INTEGERm,n,addadd=m+nENDFUNCTIONaddPROGRAMmainINTEGERaddINTEGER*,add(5,3)ENDPROGRAMmain概述外部函数定义外部函数说明外部函数调用第七讲数组的使用及Fortran函数7.6外部函数子程序苏州科技大学计算中心7.6外部函数子程序/定义外部函数只产生一个结果，由函数名返回。外部函数定义一般格式：外部函数定义→[类型说明]FUNCTION函数名（[形式参数表]）说明语句部分执行语句部分ENDFUNCTION[函数名]例如：REALFUNCTIONadd(m,n)INTEGERm,nadd=m+nENDFUNCTIONaddFUNCTIONadd(m,n)INTEGERm,n,addadd=m+nENDFUNCTIONadd概述外部函数定义外部函数说明外部函数调用第七讲数组的使用及Fortran函数7.6外部函数子程序苏州科技大学计算中心7.6外部函数子程序/外部函数说明关键字FUNCTION前面可指定内部数据类型,如：INTEGER、REAL、COMPLEX、CHARACTER或LOGICAL等。如果没指定,则按I-N规则或类型声明语句声明类型。形式参数表称为哑元表;无形参时圆括号也不能省略。形参可以是变量名、数组名或子程序名函数体位于FUNCTION语句和ENDFUNCTION语句之间。函数体有说明部分、执行部分。当返回值只有一个时