子程序上机实验（补充） .pdfVIP

program ex0808 !implicit none ! Variables logical f do i=1,1000 if (f(i)) print*,i enddo ! Body of ex0808 end program ex0808 logical function f(n) f= false i=n m=0 do while (i0) m=m+mod (i,10)**3 i=i/10 end do if (m==n) f= true end ==== program ex0815 !implicit none ! Variables complex z1,z2,z3 integer line z1=(1,1) z2=(2,1) z3=(3,1) print*,line (z1,z2,z3) ! Body of ex0815 end program ex0815 integer function line (z1,z2,z3) complex z1,z2,z3 line=0 if ((z1-z2)/abs (z1-z2)==(z3-z2)/abs (z3-z2) or.(z1-z2)/abs (z1-z2)==(z2-z3)/abs (z3-z2)) line=1 end ======== program ex0816 !implicit none ! Variables dimension a (3,3) a=reshape ((/1,2,3,4,5,6,7,8,9/),(/3,3/)) print*,calhls (a,3) ! Body of ex0816 end program ex0816 function calhls (a,n) real a (n,n) - calhls=a (1,1)*(a (2,2)*a (3,3)-a (2,3)*a (3,2))+a (1,2)*(a (3,3)*a (3,1)-a (2,1)*a (3,3)) + a (3,3)*(a (2,1)*a (3,2)-a (2,2)*a (3,1)) end ==== program ex0819 use jfs !implicit none intrinsic cos !external xinm !real jfm a=0 0 b=3 1415926/2 n=20 print*,use module print*,jfm (cos) b=1 0 print*,jfm (xinm) end program ex0819 real function jf (f,a,b,n) jf=0 0 h=(b-a)/n do i=0,n-1 jf=jf+h*f(a+i*h) enddo end function xin (x) xin=x*x end real function jf1 (f) common a,b,n jf1=0 0 h=(b-a)/n do i=0,n-1 jf1=jf1+h*f(a+i*h) enddo end !================================= ! mod0819.f90 ! module file module jfs real a,b integer n contains function xinm (x) xinm=x*x end function xinm - real func