这种参数传递如何在数值积分中实现？

wlldhyz

我有如下函数：MYFUN(X,I,J,P) = sin(X)**(m+n+1) *cos(X)**k  *nfun(X,I)*nfun(X,P)/nfun(X,J) ，与I,J,P有关，其中nfun是另一个子函数，比如
            NFUN = sin(X)*cos(X)       when I=1
            NFUN = sin(X)                  when I=2
            NFUN = 1                         when I=3
在Fortran中如何才能实现MYFUN的数值积分呢
DO I=1,3
DO J=1,3
DO P=1,3
CALL DQDAG (MYFUN,A,B,ERRABS,ERRREL,IRULE...........)
END DO
END DO
END DO
有什么好办法吗？谢谢

penultimate

楼主不是已经写出来了吗
就是把一个子程序作为一个参数进行传递呀