lit
    JINRLIB

CUR2FT -

E201

: B.Escoubes, S.de Unamuno-Escoubes :

popaa apocp co:

y = p (A + B*), = p (A'+ B'* + C'*2))
:
y = A + B*, = A'+ B'* + C'*x2.

y'=log(y). :

y = A + B*,
= A'+ B'* + C'*x2

N (xi, yi) Δyi yi.
, ,
p, a c caccoo papoca. pooc c F-pp pa, oop o cooa pop aoc po apaa C'. popaa p a . o o ooao apoca co cacc aao oo ocooac popao E202(LSQFIT).

:

: - SUBROUTINE
: - CUR2FT
: - e201.dat

:

CALL CUR2FT
Bo a o e201.dat :
  1. pa cpoa c a oa 2 , :
    co cpa o N,
    = 1 - coao apoca,
                       = 0 - apoca ooo.
  2. Bopa pa N cpo, 3 c i , i , Δyi .

e201.dat .
B a e201.dat oo oc coo aopo 1,2, .. a o c oo o apoca pa apao.
o c a oa co caa cpoa oa c, .

:

Macao co cpa o N=100.
co o o, , opaop DIMENSION p popa.

:

  1. ..,..,..,..,
    , .IV, (A-G),
    P11-83-619, , 1983, c.228-235.
:
       . . .
       CALL CUR2FT
       . . .
'e201.dat':
        5         1
       1.0D0              10.0D0               1.0D0
       2.0D0              20.0D0               2.0D0
       3.0D0              60.0D0               6.0D0
       4.0D0             270.0D0              27.0D0
       5.0D0            1800.0D0             180.0D0
        0         0
:
    EXPONENTIAL CURVE FITTING
 
       2 PARAMETERS     EXTERNAL ERRORS     INTERNAL ERRORS
    .6007441502D+00     .1048808848D+00     .4553146335D+00
    .1298860339D+01     .3162277660D-01     .1372825273D+00
 
      CHISQ=     .5653947694D+02
 
       3 PARAMETERS     EXTERNAL ERRORS     INTERNAL ERRORS
    .2007449009D+01     .2144761059D+00     .6368975761D-02
    .9311288868D-01     .1634450541D+00     .4853583030D-02
    .2009579083D+00     .2672612419D-01     .7936456905D-03
 
      CHISQ=     .1763645750D-02
 
    F TEST
 
    F=     .6411459138D+05 FOR   2 DEGREES OF FREDOM
 
        X               Y              DY              Y2              Y3
    .10000000D+01 .100000000D+02 .100000000D+01 .6683250620D+01 .9989357793D+01
    .20000000D+01 .200000000D+02 .200000000D+01 .2449489743D+02 .2003552443D+02
    .30000000D+01 .600000000D+02 .600000000D+01 .8977667218D+02 .6006392127D+02
    .40000000D+01 .270000000D+03 .270000000D+02 .3290420338D+03 .2691388983D+03
    .50000000D+01 .180000000D+04 .180000000D+03 .1205977649D+04 .1802557305D+04


home up e-mail