lit
    JINRLIB

INSAUT, INSINT - -

D221

: .. :

popaa paaa ppoa cc n o pa pa poo opa

dyi/dx = fi(x,y1,y2,...,yn), i=1,2,...,n

aoac opo aa c aao ooc.

:

: - SUBROUTINE
: - INSAUT,INSINT
: - DRKSTP(D209)
EXTERN - /

:

CALL INSAUT(N,X,Y,H,EPS,T,EXTERN,W), :

N - (INTEGER) pa ;
X - (REAL*8) aca pa,
x;
Y - (REAL*8) acc N,
y1,y2,...,yn,
- y1,y2, ...yn;
H - (REAL*8) aa a ppoa;
EPS - (REAL*8) aaa ooc ppoa;
T - (REAL*8) o paa ppoa, c T 0;
c T 0, poco opa popa oc aoo aa ppoa H;
W - (REAL*8) pao acc papoc (N,3);
EXTERN - opopa , ca EXTERNAL a popa.
:
SUBROUTINE EXTERN(X,Y,F), :
   X - (REAL*8) aca pa,
   Y,F - (REAL*8) oop acc, Y(*), F(*).
c a pa ac p pa:
F(I) = fI(X,Y(1),...,Y(N)), I=1,2,...,N.

oo o INSINT paa ppoa c [X,T] c oco ao H pop ooc.

CALL INSINT(N,X,Y,H,EPS,T,EXTERN,W), :
N,X,Y,H,EXTERN - o cc, o opopa INSAUT;
T - (REAL*8) o ppoa;
EPS - (REAL*8) apap.

p opa opopaa INSAUT INSINT
X Y(I) (I=1,2,...,N) o copa aa a
ppo cc pa pa.

:
       IMPLICIT REAL*8 (A-H,O-Z)
       DIMENSION Y(2),W(2,3)
       EXTERNAL EXTERN
       X=0.D0
       Y(1)=-1.D0
       Y(2)=1.D0
       CALL INSAUT(2,X,Y,0.1D0,0.01D0,1.D0,EXTERN,W)
       WRITE(*,12) X,Y(1),Y(2)
       X=0.D0
       Y(1)=-1.D0
       Y(2)=1.D0
       CALL INSINT(2,X,Y,0.1D0,0.01D0,1.D0,EXTERN,W)
       WRITE(*,14) X,Y(1),Y(2)
    12 FORMAT(/'  INSAUT'//3X,' X=',D10.5,5X,' Y1=',D18.10,5X,
      *' Y2=',D18.10)
    14 FORMAT(/'  INSINT'//3X,' X=',D10.5,5X,' Y1=',D18.10,5X,
      *' Y2=',D18.10)
       END
 
       SUBROUTINE EXTERN(X,Y,F)
       IMPLICIT REAL*8 (A-H,O-Z)
       DIMENSION Y(*),F(*)
       F(1)=DSIN(X)
       F(2)=DEXP(X)
       RETURN
       END
:
    INSAUT:
       X=.10000D+01      Y1=  -.5403023064D+00      Y2=   .2718281830D+01
 
    INSINT:
       X=.10000D+01      Y1=  -.5403022899D+00      Y2=   .2718281888D+01


home up e-mail