mumus 設定値

http://www-m2.ma.tum.de/Software/mumus/src/

mumus.f より

SUBROUTINES SUPPLIED BY THE USER

C   RIGHT HAND SIDE OF ODE

C     SUBROUTINE FCN   (X,Y,DY,ISIDE)

C   2-POINT AND INNER BOUNDARY-CONDITIONS

C     SUBROUTINE R     (X,YA,YB,W,NSW,ISIDE)

C   INTEGRATION-ROUTINE

C     SUBROUTINE METHOD (NDGL,FCN,XA,YM,XS,TOL,HMAX,H,ISIDE)

SUBROUTINES ESPECIALLY DESIGNED FOR THE USE IN MUMUS

C   COMPUTING OF JACOBI MATRICES

C     SUBROUTINE JACOBM

C    *    (NDGL,FCN,XA,YM,XS,TOLJAC,HMAX,H,ISIDE,

C    *     GJ(1,1),WKJ,NDGL,NDGL,NORD,ETA)

C   DECOMPOSITION OF THE LINEAR SYSTEM

C     SUBROUTINE DECGIV (...)

C   BACKWARD-SUBSTITUTION

C     SUBROUTINE BACGIV (...)

INPUT VARIABLES

C     X(M)         NODES OR INITIAL ESTIMATES FOR SWITCHING POINTS

C     Y(NY,M)      INITIAL ESTIMATES FOR Y(X(I)+)-VALUES

C     EPS          REQUIRED ACCURACY FOR NEWTON METHOD

C     TOL          REQUIRED ACCURACY FOR INTEGRATION METHOD

C     TOLJAC       REQUIRED ACCURACY FOR JACOBI MATRICES

C     ITMAX        MAXIMUM ALLOWED NUMBER OF ITERATION STEPS

C     JPRINT       PRINT PARAMETER

C     MAXBR        MAX. NUMBER OF SUCCESSIVE BROYDEN-UPDATES

C                  MAXBR = 0 IS RECOMMENDED FOR LINEAR PROBLEMS

C     COND         MAXIMUM LIMIT FOR ESTIMATE OF CONDITION NUMBER OF

C                  THE LINEAR SYSTEM

C     FC0          INITIAL RELAXATION FACTOR FOR MOD. NEWTON METHOD

C     ETA          RELATIVE DEVIATION FOR THE NUMERICAL DIFFERENTIATION

C     NCD(M)       NCD(I) = NUMBER OF BOUNDARY CONDITIONS IN

C                  SWITCHING-POINT I;

C                  MUST BE GIVEN BY THE USER IF ISW(I) .NE. 0

C     ISW(M)       INITIAL IDENTIFICATION OF NODES WITH SWITCHING POINTS

C                  ISW(I) = K MEANS: X(I) IS SW.P. K

C                  ISW(I) = 0 MEANS: X(I) IS SIMPLE NODE

C     NDGL         DIMENSION OF THE ODE SYSTEM

C     M            NUMBER OF NODES

C     NY           FIRST DIMENSION OF ARRAY Y      (NY   .GE. NDGL+1)

C     NYMM         FIRST DIMENSION OF ARRAY YMM    (NYMM .GE. NDGL)

C     LWORK(NLW)   LOGICAL WORK ARRAY

C     IWORK(NIW)   INTEGER WORK ARRAY

C     RWORK(NRW)   double precision    WORK ARRAY

C     NLW          LENGTH OF THE LOGICAL WORK ARRAY

C     NIW          LENGTH OF THE INTEGER WORK ARRAY

C     NRW          LENGTH OF THE double precision    WORK ARRAY

C                  THE MINIMUM LENGTHS OF THE WORK ARRAYS WILL BE

C                  RETURNED

OUTPUT VARIABLES

C     X(M)         SOLUTION-OR FINAL VALUES

C     Y(NY,M)      SOLUTION-OR FINAL VALUES

C     YMM(NY,M)    SOLUTION-OR FINAL VALUES FOR Y(X(I-)) VALUES

C     ISW(1:M)     FINAL IDENTIFICATION OF NODES WITH SWITCHING POINTS

C                  ISW(I) = K MEANS SW.P. K IS NODE X(I)

C     IST(1:M)     FINAL IDENTIFICATION OF SWITCHING POINTS WITH NODES

C                  IST(K) = I MEANS SW.P. K IS NODE X(I)

C     JPRINT       ERROR INDICATOR

FORMAL PARAMETERS OF THE USER SUPPLIED ROUTINES

C  IN FCN (X,Y,DY,ISIDE)

C     X

C     Y (1:NDGL)

C     DY(1:NDGL)   DY = DY/DX = F(X,Y)

C     ISIDE(1:M)   ISIDE(K) = -1 MEANS: (X - SW.P. K)*(X(M)-X(1) .LE. 0

C                  ISIDE(K) =  1 MEANS: (X - SW.P. K)*(X(M)-X(1) .GT. 0

C

C  IN R (X,YA,YB,W,NSW,ISIDE)

C     NSW          NUMBER OF SW.P.

C                  IF NSW .EQ. 0 THEN X=X(M) IS ASSUMED

C     YA (1:NDGL)  Y+  IF X .NE. X(M)

C     YB (1:NDGL)  Y-  IF X .NE. X(M)

C     YA           Y(X(1)) IF X .EQ. X(M)

C     YB           Y(X(M)) IF X .EQ. X(M)

C                  X, YA, YB MAY NOT BE ALTERED BY THE USER!

C     W  (1:2N+1)  W = S(X,YA,YB)

C

C IN METHOD (NDGL,FCN,XA,YM,XS,TOL,HMAX,H,ISIDE)

C     METHOD INTEGRATES FROM XA TO XS

C     H            INITIAL STEPSIZE

C     TOL          REQUIRED ACCURACY