24小时热门版块排行榜    

查看: 532  |  回复: 9
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 开心@豆豆 的 5 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

开心@豆豆

[交流] 【求助】fortran小程序,帮忙一下

帮忙编一下
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cxz314

木虫 (正式写手)

★ ★ ★ ★ ★
开心@豆豆(金币+5,VIP+0):谢谢,fm是单位 9-5 12:35
PROGRAM LZ
        IMPLICIT NONE

        REAL*8 :: SQRC,SQRLD,DELTSQRC,DELTSQRLD
        REAL*8 :: SQRC_2008,SQRLD_2008
        REAL*8 :: A,SQR0,F,M,R0,FINAL

        F=1.D0
        M=2.D0
        R0=1.2D0*F*M

        PRINT*,'ENTER THE VALUE OF A'
        READ*, A

        CALL FSQRC(A,SQRC)
        CALL FSQRLD(A,R0,SQRLD)
        CALL FSQRC(2008,SQRC_2008)
        CALL FSQRLD(2008,R0,SQRLD_2008)
       
        DELTSQRC=SQRC-SQRC_2008
        DELTSQRLD=SQRLD-SQRLD_2008

        FINAL=DELTSQRC-DELTSQRLD

        PRINT*,'RESULT=',FINAL

        END


        SUBROUTINE FSQRC(X,Y)
        IMPLICIT NONE
        REAL*8 :: X,Y

        Y=X*X

        END SUBROUTINE

        SUBROUTINE FSQRLD(X,Y,Z)
        IMPLICIT NONE
        REAL*8 :: X,Y,Z

        Z=(3.D0/5.D0)*Y*Y*(X**(2.D0/3.D0))

        END SUBROUTINE
lz不知道你的fm是什么,所以这里取F=1,M=2
函数Rc2(A)的表达式也不知道,这里取Rc2(A)=A*A

[ Last edited by cxz314 on 2009-9-5 at 12:28 ]
9楼2009-09-05 12:26:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 开心@豆豆 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见