Ideone.com

fork download

copy

!! PROGRAM coupled_oscillators !!
!! n 個の質点からなる連成振動子の運動を         !!
!! 0 < t < TEND の範囲で解く。                  !!
!! 長さ 2*N の配列 w[0~(N-1)][0~1] を定義して， !!
!! j 番目の質点の変位を w[j][0]，               !!
!! 速度を w[j][1] に格納する。                  !!
!! 質点の質量の逆数に比例する無次元パラメータは !!
!! 配列 c[0~(N-1)] に格納する。                 !!
!! 初期条件は変位ゼロ，                         !!
!! 速度は 1 番目の質点のみ VZERO とする。       !!
program coupled_oscillators
implicit none
integer(4), parameter :: n = 10
integer(4) :: j
real(8) :: t, dt
real(8), dimension(n) :: c
real(8), dimension(n, 2) :: w, p, k1, k2, k3, k4
real(8), parameter :: vzero = 0.1d0, tend = 100.0d0
!! 時間ステップの刻み幅 dt を標準入力から読み込む。!!
read(*, *) dt
!! 質量はみんな同じ。!!
c(:) = 1.0d0
!! 未知変数 w と時刻 t の初期化。!!
w(:, 1) = 0.0d0
w(1, 2) = vzero
w(2:n, 2) = 0.0d0
t = 0.0d0
!! 時間積分の開始(無限 do loop)。!!
do
  !! 時刻 t での質点の位置を書き出す。!!
  write(*, "(f12.7, 99e16.8)") t, &
    (w(j, 1) + dble(j)/dble(n+1), j = 1, n)
  !! 時刻か゛ tend を超えたら無限 do loop から脱出する。!!
  if ( t >= tend ) stop
  !! 4次ルンゲ・クッタ法。微分方程式を !!
  !! dw/dt = f(w) だと思って，機械的に !!
  !! アルゴリズムを適用する。!!
  !! まず補助変数 k1 = f(p) を計算する。!!
  p(:, :) = w(:, :)
  call rhs (k1, p, c, n)
  !! 補助変数 k2 = f(p) を計算する。!!
  p(:, :) = w(:, :) + 0.5d0 * dt * k1(:, :)
  call rhs (k2, p, c, n)
  !! 補助変数 k3 = f(p) を計算する。!!
  p(:, :) = w(:, :) + 0.5d0 * dt * k2(:, :)
  call rhs (k3, p, c, n)
  !! 補助変数 k4 = f(p) を計算する。!!
  p(:, :) = w(:, :) + dt * k3(:, :)
  call rhs (k4, p, c, n)
  !! 未知変数 w と時刻 t を次の時間ステップでの値で上書きする。!!
  w(:, :) = w(:, :) + (dt/6.0d0) * &
    (k1(:, :) + k4(:, :) + 2.d0*(k2(:, :) + k3(:, :)))
  t = t + dt
end do
end
 
!! 微分方程式 dw/dt = f(w) の右辺を計算するサブルーチン。!!
subroutine rhs (f, w, c, n)
  implicit none
  integer(4) :: n, j
  real(8), dimension(n) :: c
  real(8), dimension(n, 2) :: f, w
    do j = 1, n
      !! du(j)/dt = v(j) !!
      f(j, 1) = w(j, 2)
      !! dv(j)/dt = -c(j) * (2u(j) - u(j-1) - u(j+1)) !!
      if ( j == 1 ) then
        f(j, 2) = - c(j) * (2.d0*w(j, 1) - w(j+1, 1))
      else if ( j == n ) then
        f(j, 2) = - c(j) * (2.d0*w(j, 1) - w(j-1, 1))
      else
        f(j, 2) = - c(j) * (2.d0*w(j, 1) - w(j-1, 1) - w(j+1, 1))
    end if
  end do
  return
end subroutine rhs

ISEgUFJPR1JBTSBjb3VwbGVkX29zY2lsbGF0b3JzICEhCiEhIG4g5YCL44Gu6LOq54K544GL44KJ44Gq44KL6YCj5oiQ5oyv5YuV5a2Q44Gu6YGL5YuV44KSICAgICAgICAgISEKISEgMCA8IHQgPCBURU5EIOOBruevhOWbsuOBp+ino+OBj+OAgiAgICAgICAgICAgICAgICAgICEhCiEhIOmVt+OBlSAyKk4g44Gu6YWN5YiXIHdbMH4oTi0xKV1bMH4xXSDjgpLlrprnvqnjgZfjgabvvIwgISEKISEgaiDnlarnm67jga7os6rngrnjga7lpInkvY3jgpIgd1tqXVswXe+8jCAgICAgICAgICAgICAgICEhCiEhIOmAn+W6puOCkiB3W2pdWzFdIOOBq+agvOe0jeOBmeOCi+OAgiAgICAgICAgICAgICAgICAgICEhCiEhIOizqueCueOBruizqumHj+OBrumAhuaVsOOBq+avlOS+i+OBmeOCi+eEoeasoeWFg+ODkeODqeODoeODvOOCv+OBryAhIQohISDphY3liJcgY1swfihOLTEpXSDjgavmoLzntI3jgZnjgovjgIIgICAgICAgICAgICAgICAgICEhCiEhIOWIneacn+adoeS7tuOBr+WkieS9jeOCvOODre+8jCAgICAgICAgICAgICAgICAgICAgICAgICAhIQohISDpgJ/luqbjga8gMSDnlarnm67jga7os6rngrnjga7jgb8gVlpFUk8g44Go44GZ44KL44CCICAgICAgICEhCnByb2dyYW0gY291cGxlZF9vc2NpbGxhdG9ycwppbXBsaWNpdCBub25lCmludGVnZXIoNCksIHBhcmFtZXRlciA6OiBuID0gMTAKaW50ZWdlcig0KSA6OiBqCnJlYWwoOCkgOjogdCwgZHQKcmVhbCg4KSwgZGltZW5zaW9uKG4pIDo6IGMKcmVhbCg4KSwgZGltZW5zaW9uKG4sIDIpIDo6IHcsIHAsIGsxLCBrMiwgazMsIGs0CnJlYWwoOCksIHBhcmFtZXRlciA6OiB2emVybyA9IDAuMWQwLCB0ZW5kID0gMTAwLjBkMAohISDmmYLplpPjgrnjg4bjg4Pjg5fjga7liLvjgb/luYUgZHQg44KS5qiZ5rqW5YWl5Yqb44GL44KJ6Kqt44G/6L6844KA44CCISEKcmVhZCgqLCAqKSBkdAohISDos6rph4/jga/jgb/jgpPjgarlkIzjgZjjgIIhIQpjKDopID0gMS4wZDAKISEg5pyq55+l5aSJ5pWwIHcg44Go5pmC5Yi7IHQg44Gu5Yid5pyf5YyW44CCISEKdyg6LCAxKSA9IDAuMGQwCncoMSwgMikgPSB2emVybwp3KDI6biwgMikgPSAwLjBkMAp0ID0gMC4wZDAKISEg5pmC6ZaT56mN5YiG44Gu6ZaL5aeLKOeEoemZkCBkbyBsb29wKeOAgiEhCmRvCiAgISEg5pmC5Yi7IHQg44Gn44Gu6LOq54K544Gu5L2N572u44KS5pu444GN5Ye644GZ44CCISEKICB3cml0ZSgqLCAiKGYxMi43LCA5OWUxNi44KSIpIHQsICYKICAgICh3KGosIDEpICsgZGJsZShqKS9kYmxlKG4rMSksIGogPSAxLCBuKQogICEhIOaZguWIu+OBi+OCmyB0ZW5kIOOCkui2heOBiOOBn+OCieeEoemZkCBkbyBsb29wIOOBi+OCieiEseWHuuOBmeOCi+OAgiEhCiAgaWYgKCB0ID49IHRlbmQgKSBzdG9wCiAgISEgNOasoeODq+ODs+OCsuODu+OCr+ODg+OCv+azleOAguW+ruWIhuaWueeoi+W8j+OCkiAhIQogICEhIGR3L2R0ID0gZih3KSDjgaDjgajmgJ3jgaPjgabvvIzmqZ/morDnmoTjgasgISEKICAhISDjgqLjg6vjgrTjg6rjgrrjg6DjgpLpgannlKjjgZnjgovjgIIhIQogICEhIOOBvuOBmuijnOWKqeWkieaVsCBrMSA9IGYocCkg44KS6KiI566X44GZ44KL44CCISEKICBwKDosIDopID0gdyg6LCA6KQogIGNhbGwgcmhzIChrMSwgcCwgYywgbikKICAhISDoo5zliqnlpInmlbAgazIgPSBmKHApIOOCkuioiOeul+OBmeOCi+OAgiEhCiAgcCg6LCA6KSA9IHcoOiwgOikgKyAwLjVkMCAqIGR0ICogazEoOiwgOikKICBjYWxsIHJocyAoazIsIHAsIGMsIG4pCiAgISEg6KOc5Yqp5aSJ5pWwIGszID0gZihwKSDjgpLoqIjnrpfjgZnjgovjgIIhIQogIHAoOiwgOikgPSB3KDosIDopICsgMC41ZDAgKiBkdCAqIGsyKDosIDopCiAgY2FsbCByaHMgKGszLCBwLCBjLCBuKQogICEhIOijnOWKqeWkieaVsCBrNCA9IGYocCkg44KS6KiI566X44GZ44KL44CCISEKICBwKDosIDopID0gdyg6LCA6KSArIGR0ICogazMoOiwgOikKICBjYWxsIHJocyAoazQsIHAsIGMsIG4pCiAgISEg5pyq55+l5aSJ5pWwIHcg44Go5pmC5Yi7IHQg44KS5qyh44Gu5pmC6ZaT44K544OG44OD44OX44Gn44Gu5YCk44Gn5LiK5pu444GN44GZ44KL44CCISEKICB3KDosIDopID0gdyg6LCA6KSArIChkdC82LjBkMCkgKiAmCiAgICAoazEoOiwgOikgKyBrNCg6LCA6KSArIDIuZDAqKGsyKDosIDopICsgazMoOiwgOikpKQogIHQgPSB0ICsgZHQKZW5kIGRvCmVuZAoKISEg5b6u5YiG5pa556iL5byPIGR3L2R0ID0gZih3KSDjga7lj7PovrrjgpLoqIjnrpfjgZnjgovjgrXjg5bjg6vjg7zjg4Hjg7PjgIIhIQpzdWJyb3V0aW5lIHJocyAoZiwgdywgYywgbikKICBpbXBsaWNpdCBub25lCiAgaW50ZWdlcig0KSA6OiBuLCBqCiAgcmVhbCg4KSwgZGltZW5zaW9uKG4pIDo6IGMKICByZWFsKDgpLCBkaW1lbnNpb24obiwgMikgOjogZiwgdwogICAgZG8gaiA9IDEsIG4KICAgICAgISEgZHUoaikvZHQgPSB2KGopICEhCiAgICAgIGYoaiwgMSkgPSB3KGosIDIpCiAgICAgICEhIGR2KGopL2R0ID0gLWMoaikgKiAoMnUoaikgLSB1KGotMSkgLSB1KGorMSkpICEhCiAgICAgIGlmICggaiA9PSAxICkgdGhlbgogICAgICAgIGYoaiwgMikgPSAtIGMoaikgKiAoMi5kMCp3KGosIDEpIC0gdyhqKzEsIDEpKQogICAgICBlbHNlIGlmICggaiA9PSBuICkgdGhlbgogICAgICAgIGYoaiwgMikgPSAtIGMoaikgKiAoMi5kMCp3KGosIDEpIC0gdyhqLTEsIDEpKQogICAgICBlbHNlCiAgICAgICAgZihqLCAyKSA9IC0gYyhqKSAqICgyLmQwKncoaiwgMSkgLSB3KGotMSwgMSkgLSB3KGorMSwgMSkpCiAgICBlbmQgaWYKICBlbmQgZG8KICByZXR1cm4KZW5kIHN1YnJvdXRpbmUgcmhzCg==

Success #stdin #stdout 0s 4396KB

stdin

copy

1.0

stdout

copy

   0.0000000  0.90909091E-01  0.18181818E+00  0.27272727E+00  0.36363636E+00  0.45454545E+00  0.54545455E+00  0.63636364E+00  0.72727273E+00  0.81818182E+00  0.90909091E+00
   1.0000000  0.15757576E+00  0.19848485E+00  0.27272727E+00  0.36363636E+00  0.45454545E+00  0.54545455E+00  0.63636364E+00  0.72727273E+00  0.81818182E+00  0.90909091E+00
   2.0000000  0.12979798E+00  0.23459596E+00  0.28939394E+00  0.36502525E+00  0.45454545E+00  0.54545455E+00  0.63636364E+00  0.72727273E+00  0.81818182E+00  0.90909091E+00
   3.0000000  0.83270202E-01  0.22391888E+00  0.32099116E+00  0.38041877E+00  0.45674453E+00  0.54554135E+00  0.63636364E+00  0.72727273E+00  0.81818182E+00  0.90909091E+00
   4.0000000  0.80395974E-01  0.18078616E+00  0.31699811E+00  0.40799400E+00  0.47126999E+00  0.54831913E+00  0.63657583E+00  0.72727755E+00  0.81818182E+00  0.90909091E+00
   5.0000000  0.93318347E-01  0.16812310E+00  0.27679669E+00  0.40832917E+00  0.49605539E+00  0.56206691E+00  0.63975930E+00  0.72762919E+00  0.81819890E+00  0.90909116E+00
   6.0000000  0.94057165E-01  0.18166749E+00  0.25711664E+00  0.37247170E+00  0.49899718E+00  0.58469520E+00  0.65284060E+00  0.73110532E+00  0.81868894E+00  0.90912731E+00
   7.0000000  0.90300954E-01  0.18650534E+00  0.26932388E+00  0.34786568E+00  0.46731348E+00  0.58933146E+00  0.67374851E+00  0.74360411E+00  0.82237961E+00  0.90974557E+00
   8.0000000  0.90101789E-01  0.18192345E+00  0.27848602E+00  0.35673236E+00  0.43967455E+00  0.56143589E+00  0.67950857E+00  0.76309880E+00  0.83435820E+00  0.91350644E+00
   9.0000000  0.90989975E-01  0.18047175E+00  0.27422701E+00  0.36937404E+00  0.44451462E+00  0.53212125E+00  0.65494487E+00  0.76960791E+00  0.85254714E+00  0.92417819E+00
  10.0000000  0.91103419E-01  0.18177060E+00  0.27085127E+00  0.36692858E+00  0.45934980E+00  0.53286481E+00  0.62489213E+00  0.74779256E+00  0.85862643E+00  0.93741754E+00
  11.0000000  0.90952645E-01  0.18212808E+00  0.27223798E+00  0.36170511E+00  0.45958191E+00  0.54864451E+00  0.62162736E+00  0.71663333E+00  0.83535872E+00  0.93406123E+00
  12.0000000  0.90850903E-01  0.18187813E+00  0.27324124E+00  0.36241624E+00  0.45307940E+00  0.55167652E+00  0.63625268E+00  0.70598384E+00  0.79520722E+00  0.90145183E+00
  13.0000000  0.90868364E-01  0.18176632E+00  0.27287882E+00  0.36417932E+00  0.45223265E+00  0.54354430E+00  0.63826043E+00  0.71085012E+00  0.77082403E+00  0.86457978E+00
  14.0000000  0.90935918E-01  0.18175710E+00  0.27257674E+00  0.36369385E+00  0.45338343E+00  0.53705581E+00  0.62186911E+00  0.70535286E+00  0.77659634E+00  0.86608069E+00
  15.0000000  0.90918081E-01  0.18176348E+00  0.27228748E+00  0.36182748E+00  0.44943994E+00  0.53017959E+00  0.60405331E+00  0.69185977E+00  0.79961068E+00  0.90538802E+00
  16.0000000  0.90819467E-01  0.18143465E+00  0.27095348E+00  0.35759923E+00  0.43934520E+00  0.51791307E+00  0.59855794E+00  0.69692998E+00  0.82185523E+00  0.93383551E+00
  17.0000000  0.90459252E-01  0.17993868E+00  0.26674306E+00  0.34864628E+00  0.42636235E+00  0.50907451E+00  0.60998323E+00  0.72398436E+00  0.83200629E+00  0.92337331E+00
  18.0000000  0.89053041E-01  0.17570807E+00  0.25784732E+00  0.33623771E+00  0.41868576E+00  0.51718783E+00  0.63346426E+00  0.74476694E+00  0.82646768E+00  0.90078016E+00
  19.0000000  0.85257094E-01  0.16716270E+00  0.24577954E+00  0.32836738E+00  0.42636405E+00  0.54075850E+00  0.65161236E+00  0.74001591E+00  0.81447643E+00  0.89847039E+00
  20.0000000  0.78367276E-01  0.15592029E+00  0.23812137E+00  0.33516782E+00  0.44864964E+00  0.56044532E+00  0.64940723E+00  0.72346734E+00  0.81097934E+00  0.91032436E+00
  21.0000000  0.71194214E-01  0.14981089E+00  0.24471682E+00  0.35668164E+00  0.46864416E+00  0.55968661E+00  0.63443020E+00  0.71809245E+00  0.81741546E+00  0.91470037E+00
  22.0000000  0.71665902E-01  0.15944859E+00  0.26683418E+00  0.37755164E+00  0.46970875E+00  0.54481571E+00  0.62723273E+00  0.72571287E+00  0.82207047E+00  0.90970421E+00
  23.0000000  0.87316591E-01  0.18702268E+00  0.29154159E+00  0.38166932E+00  0.45601753E+00  0.53631193E+00  0.63350164E+00  0.73175243E+00  0.81968463E+00  0.90660688E+00
  24.0000000  0.11359868E+00  0.21818376E+00  0.30320181E+00  0.37228169E+00  0.44777354E+00  0.54216820E+00  0.64069219E+00  0.72951218E+00  0.81651268E+00  0.90840448E+00
  25.0000000  0.13040737E+00  0.23073261E+00  0.30076687E+00  0.36890351E+00  0.45596006E+00  0.55173502E+00  0.64013011E+00  0.72576629E+00  0.81715425E+00  0.91007313E+00
  26.0000000  0.11946297E+00  0.21551801E+00  0.29614800E+00  0.38172222E+00  0.47189198E+00  0.55572256E+00  0.63757992E+00  0.72661150E+00  0.81902736E+00  0.90956541E+00
  27.0000000  0.88111569E-01  0.18630614E+00  0.29459050E+00  0.39772371E+00  0.48273534E+00  0.55858494E+00  0.64114865E+00  0.73042334E+00  0.81963325E+00  0.90898270E+00
  28.0000000  0.66552914E-01  0.16607271E+00  0.28761781E+00  0.39725514E+00  0.48535743E+00  0.56692569E+00  0.65057761E+00  0.73476210E+00  0.82065547E+00  0.90972328E+00
  29.0000000  0.74211624E-01  0.16570005E+00  0.27076134E+00  0.37767486E+00  0.48099943E+00  0.57617583E+00  0.66069634E+00  0.74102598E+00  0.82457468E+00  0.91165347E+00
  30.0000000  0.96677674E-01  0.17792640E+00  0.25677830E+00  0.35493980E+00  0.46840605E+00  0.57532339E+00  0.66683593E+00  0.74994358E+00  0.83173796E+00  0.91497205E+00
  31.0000000  0.10573063E+00  0.18844017E+00  0.26015369E+00  0.34597789E+00  0.45030155E+00  0.56053323E+00  0.66458258E+00  0.75705712E+00  0.84013713E+00  0.91995468E+00
  32.0000000  0.95210937E-01  0.18921054E+00  0.27545862E+00  0.35383625E+00  0.43843081E+00  0.54034179E+00  0.65128869E+00  0.75502234E+00  0.84505693E+00  0.92483842E+00
  33.0000000  0.83650248E-01  0.18292158E+00  0.28368263E+00  0.36747381E+00  0.44236441E+00  0.52860505E+00  0.63165640E+00  0.74014534E+00  0.83984526E+00  0.92495829E+00
  34.0000000  0.85060252E-01  0.17781667E+00  0.27741320E+00  0.37321265E+00  0.45611357E+00  0.53219259E+00  0.61740962E+00  0.71728674E+00  0.82086524E+00  0.91541804E+00
  35.0000000  0.92869549E-01  0.17866287E+00  0.26795803E+00  0.36748450E+00  0.46359887E+00  0.54384281E+00  0.61646613E+00  0.69809627E+00  0.79392824E+00  0.89683646E+00
  36.0000000  0.94975743E-01  0.18254276E+00  0.26702330E+00  0.35866232E+00  0.45710403E+00  0.54841187E+00  0.62333778E+00  0.69208201E+00  0.77426251E+00  0.87918705E+00
  37.0000000  0.91050034E-01  0.18377974E+00  0.27196185E+00  0.35552149E+00  0.44421436E+00  0.53832006E+00  0.62460498E+00  0.69848228E+00  0.77605897E+00  0.87681798E+00
  38.0000000  0.88346389E-01  0.18110038E+00  0.27293777E+00  0.35648545E+00  0.43565412E+00  0.52140919E+00  0.61530901E+00  0.70871303E+00  0.79918826E+00  0.89502426E+00
  39.0000000  0.89028570E-01  0.17744358E+00  0.26679523E+00  0.35349366E+00  0.43261452E+00  0.51190401E+00  0.60636892E+00  0.71698661E+00  0.82704208E+00  0.92101471E+00
  40.0000000  0.89086296E-01  0.17423649E+00  0.25808263E+00  0.34411089E+00  0.43017675E+00  0.51621942E+00  0.61223554E+00  0.72501200E+00  0.84012390E+00  0.93294147E+00
  41.0000000  0.85834598E-01  0.16967953E+00  0.25099110E+00  0.33515147E+00  0.42872243E+00  0.53004561E+00  0.63283537E+00  0.73455088E+00  0.83288958E+00  0.92169234E+00
  42.0000000  0.80843006E-01  0.16308773E+00  0.24669965E+00  0.33492714E+00  0.43486375E+00  0.54568628E+00  0.65205149E+00  0.74016956E+00  0.81682484E+00  0.90133145E+00
  43.0000000  0.77117858E-01  0.15823562E+00  0.24738614E+00  0.34561664E+00  0.45062927E+00  0.55711555E+00  0.65473891E+00  0.73523849E+00  0.80756491E+00  0.89361579E+00
  44.0000000  0.77338457E-01  0.16122148E+00  0.25700050E+00  0.36276567E+00  0.46719915E+00  0.55983179E+00  0.64204401E+00  0.72341098E+00  0.81054906E+00  0.90366585E+00
  45.0000000  0.84031362E-01  0.17547462E+00  0.27580925E+00  0.37866119E+00  0.47285078E+00  0.55274469E+00  0.62848706E+00  0.71707511E+00  0.81906639E+00  0.91646445E+00
  46.0000000  0.97685770E-01  0.19791076E+00  0.29659754E+00  0.38635270E+00  0.46568591E+00  0.54216732E+00  0.62632662E+00  0.72251276E+00  0.82371902E+00  0.91735122E+00
  47.0000000  0.11327760E+00  0.21846095E+00  0.30886164E+00  0.38439914E+00  0.45614056E+00  0.53889894E+00  0.63500002E+00  0.73218430E+00  0.82171126E+00  0.90886733E+00
  48.0000000  0.12077875E+00  0.22471079E+00  0.30726045E+00  0.37907389E+00  0.45661241E+00  0.54770280E+00  0.64491336E+00  0.73527124E+00  0.81753906E+00  0.90372048E+00
  49.0000000  0.11242072E+00  0.21086445E+00  0.29554852E+00  0.37884237E+00  0.46919761E+00  0.56194328E+00  0.64903678E+00  0.73162650E+00  0.81678213E+00  0.90708596E+00
  50.0000000  0.91388138E-01  0.18442381E+00  0.28218049E+00  0.38447762E+00  0.48369534E+00  0.57121777E+00  0.64933673E+00  0.73040040E+00  0.82064643E+00  0.91276210E+00
  51.0000000  0.72315402E-01  0.16259540E+00  0.27265290E+00  0.38673620E+00  0.48737205E+00  0.57209925E+00  0.65215706E+00  0.73712816E+00  0.82635997E+00  0.91457121E+00
  52.0000000  0.69924032E-01  0.15906777E+00  0.26688660E+00  0.37674514E+00  0.47643179E+00  0.56833982E+00  0.65893600E+00  0.74749652E+00  0.83138001E+00  0.91421817E+00
  53.0000000  0.84921855E-01  0.17266745E+00  0.26352075E+00  0.35801007E+00  0.45821824E+00  0.56259848E+00  0.66335554E+00  0.75367188E+00  0.83546395E+00  0.91630146E+00
  54.0000000  0.10225839E+00  0.18919702E+00  0.26403483E+00  0.34492635E+00  0.44367093E+00  0.55313938E+00  0.65802758E+00  0.75193834E+00  0.83829460E+00  0.92054363E+00
  55.0000000  0.10576716E+00  0.19488779E+00  0.27013877E+00  0.34808577E+00  0.43905832E+00  0.53983848E+00  0.64266347E+00  0.74272578E+00  0.83676706E+00  0.92200506E+00
  56.0000000  0.94509367E-01  0.18806553E+00  0.27831845E+00  0.36287268E+00  0.44379624E+00  0.52900471E+00  0.62479214E+00  0.72746087E+00  0.82673361E+00  0.91673999E+00
  57.0000000  0.82467972E-01  0.17810758E+00  0.28103562E+00  0.37432681E+00  0.45268673E+00  0.52801206E+00  0.61325723E+00  0.70911957E+00  0.80816293E+00  0.90515829E+00
  58.0000000  0.81976200E-01  0.17448566E+00  0.27512761E+00  0.37238147E+00  0.45863508E+00  0.53578351E+00  0.61157450E+00  0.69421716E+00  0.78806834E+00  0.89162601E+00
  59.0000000  0.90632443E-01  0.17806625E+00  0.26634768E+00  0.36047894E+00  0.45594400E+00  0.54203863E+00  0.61634027E+00  0.69011161E+00  0.77746994E+00  0.88291993E+00
  60.0000000  0.96503338E-01  0.18246197E+00  0.26270185E+00  0.34950683E+00  0.44460491E+00  0.53762898E+00  0.62077995E+00  0.69873963E+00  0.78403739E+00  0.88544834E+00
  61.0000000  0.93258661E-01  0.18180356E+00  0.26458560E+00  0.34569563E+00  0.43150910E+00  0.52457400E+00  0.62070825E+00  0.71426754E+00  0.80564637E+00  0.90026424E+00
  62.0000000  0.85865681E-01  0.17592766E+00  0.26480552E+00  0.34607979E+00  0.42515048E+00  0.51463838E+00  0.61865640E+00  0.72799162E+00  0.83007844E+00  0.91950888E+00
  63.0000000  0.81901260E-01  0.16883562E+00  0.25844267E+00  0.34482013E+00  0.42840552E+00  0.51809973E+00  0.62174029E+00  0.73514835E+00  0.84255630E+00  0.93013164E+00
  64.0000000  0.82157619E-01  0.16407497E+00  0.24957222E+00  0.34160099E+00  0.43756076E+00  0.53433212E+00  0.63360951E+00  0.73656447E+00  0.83660966E+00  0.92436979E+00
  65.0000000  0.82384471E-01  0.16275297E+00  0.24686118E+00  0.34240690E+00  0.44794407E+00  0.55288815E+00  0.64835268E+00  0.73483555E+00  0.81934770E+00  0.90765930E+00
  66.0000000  0.81336205E-01  0.16532871E+00  0.25472350E+00  0.35244126E+00  0.45798972E+00  0.56276378E+00  0.65452034E+00  0.73095193E+00  0.80553862E+00  0.89479035E+00
  67.0000000  0.83102516E-01  0.17338041E+00  0.27049916E+00  0.36950364E+00  0.46699014E+00  0.56057379E+00  0.64674184E+00  0.72543870E+00  0.80505719E+00  0.89640105E+00
  68.0000000  0.91424802E-01  0.18792811E+00  0.28829186E+00  0.38502351E+00  0.47189642E+00  0.55135325E+00  0.63241297E+00  0.72110935E+00  0.81525212E+00  0.90897241E+00
  69.0000000  0.10418686E+00  0.20577278E+00  0.30262490E+00  0.39142385E+00  0.46986662E+00  0.54349006E+00  0.62516114E+00  0.72199057E+00  0.82520347E+00  0.91908765E+00
  70.0000000  0.11430811E+00  0.21869572E+00  0.30899432E+00  0.38823345E+00  0.46371512E+00  0.54314578E+00  0.63162969E+00  0.72864197E+00  0.82689969E+00  0.91769181E+00
  71.0000000  0.11500361E+00  0.21808497E+00  0.30461822E+00  0.38140949E+00  0.46143436E+00  0.55134356E+00  0.64584212E+00  0.73630702E+00  0.82194128E+00  0.90884943E+00
  72.0000000  0.10443734E+00  0.20193343E+00  0.29103017E+00  0.37728476E+00  0.46803789E+00  0.56382242E+00  0.65659590E+00  0.73967855E+00  0.81808464E+00  0.90377908E+00
  73.0000000  0.87435988E-01  0.17808308E+00  0.27482830E+00  0.37700748E+00  0.47879355E+00  0.57337361E+00  0.65877059E+00  0.73905828E+00  0.82072323E+00  0.90796025E+00
  74.0000000  0.73671638E-01  0.16011808E+00  0.26366504E+00  0.37616937E+00  0.48261173E+00  0.57424771E+00  0.65621238E+00  0.73974589E+00  0.82839572E+00  0.91639370E+00
  75.0000000  0.72005439E-01  0.15819269E+00  0.26082018E+00  0.36977302E+00  0.47283371E+00  0.56601966E+00  0.65459796E+00  0.74474794E+00  0.83553894E+00  0.92107453E+00
  76.0000000  0.83464609E-01  0.17155548E+00  0.26403965E+00  0.35824211E+00  0.45413204E+00  0.55332829E+00  0.65382385E+00  0.74983009E+00  0.83786202E+00  0.91996635E+00
  77.0000000  0.99006979E-01  0.18897636E+00  0.26919464E+00  0.34864938E+00  0.43854759E+00  0.54161640E+00  0.64863316E+00  0.74731511E+00  0.83447872E+00  0.91676783E+00
  78.0000000  0.10612500E+00  0.19742468E+00  0.27383718E+00  0.34883973E+00  0.43508119E+00  0.53352851E+00  0.63589692E+00  0.73417787E+00  0.82615348E+00  0.91389269E+00
  79.0000000  0.99735735E-01  0.19222400E+00  0.27703708E+00  0.35920432E+00  0.44290270E+00  0.52934906E+00  0.61968894E+00  0.71522273E+00  0.81358503E+00  0.90916417E+00
  80.0000000  0.86856676E-01  0.17996146E+00  0.27748326E+00  0.37069658E+00  0.45337622E+00  0.52911770E+00  0.60845337E+00  0.69885728E+00  0.79851395E+00  0.90022123E+00
  81.0000000  0.79646277E-01  0.17159984E+00  0.27378546E+00  0.37241883E+00  0.45751574E+00  0.53226508E+00  0.60746352E+00  0.69122753E+00  0.78580781E+00  0.88995038E+00
  82.0000000  0.83433867E-01  0.17239652E+00  0.26687715E+00  0.36182462E+00  0.45202693E+00  0.53559376E+00  0.61443016E+00  0.69406566E+00  0.78253232E+00  0.88549872E+00
  83.0000000  0.92088661E-01  0.17813849E+00  0.26048985E+00  0.34693652E+00  0.44028172E+00  0.53441864E+00  0.62233651E+00  0.70554312E+00  0.79304558E+00  0.89184753E+00
  84.0000000  0.95332380E-01  0.18058275E+00  0.25779415E+00  0.33823933E+00  0.42926648E+00  0.52775188E+00  0.62619963E+00  0.72115961E+00  0.81403398E+00  0.90683046E+00
  85.0000000  0.89635324E-01  0.17584379E+00  0.25792957E+00  0.33911344E+00  0.42532781E+00  0.52150566E+00  0.62708140E+00  0.73471432E+00  0.83471162E+00  0.92204208E+00
  86.0000000  0.80913272E-01  0.16743000E+00  0.25731609E+00  0.34478753E+00  0.43086152E+00  0.52417677E+00  0.62988658E+00  0.74100391E+00  0.84338228E+00  0.92838171E+00
  87.0000000  0.77161795E-01  0.16212168E+00  0.25488635E+00  0.34970626E+00  0.44330394E+00  0.53816128E+00  0.63747160E+00  0.73905407E+00  0.83580119E+00  0.92217231E+00
  88.0000000  0.80422292E-01  0.16407079E+00  0.25482731E+00  0.35401924E+00  0.45703680E+00  0.55599128E+00  0.64682410E+00  0.73246179E+00  0.81858697E+00  0.90806059E+00
  89.0000000  0.86850519E-01  0.17284206E+00  0.26276411E+00  0.36198181E+00  0.46703969E+00  0.56620496E+00  0.65113802E+00  0.72612975E+00  0.80448830E+00  0.89635697E+00
  90.0000000  0.92931042E-01  0.18558983E+00  0.27919116E+00  0.37504480E+00  0.47147902E+00  0.56332509E+00  0.64628084E+00  0.72284928E+00  0.80287064E+00  0.89559182E+00
  91.0000000  0.98982417E-01  0.19933444E+00  0.29750671E+00  0.38826801E+00  0.47143075E+00  0.55227496E+00  0.63575334E+00  0.72293271E+00  0.81296762E+00  0.90552043E+00
  92.0000000  0.10602640E+00  0.21077725E+00  0.30879938E+00  0.39429029E+00  0.46902118E+00  0.54360345E+00  0.62886803E+00  0.72593071E+00  0.82544877E+00  0.91713294E+00
  93.0000000  0.11135568E+00  0.21542267E+00  0.30826235E+00  0.39036107E+00  0.46645918E+00  0.54474237E+00  0.63290054E+00  0.73133813E+00  0.83087620E+00  0.92072040E+00
  94.0000000  0.10958493E+00  0.20916453E+00  0.29741825E+00  0.38083518E+00  0.46605584E+00  0.55503789E+00  0.64628979E+00  0.73770283E+00  0.82755496E+00  0.91482038E+00
  95.0000000  0.98558447E-01  0.19224936E+00  0.28176423E+00  0.37276147E+00  0.46915377E+00  0.56749216E+00  0.65975147E+00  0.74252703E+00  0.82189497E+00  0.90711891E+00
  96.0000000  0.83239034E-01  0.17165440E+00  0.26741892E+00  0.36939846E+00  0.47379242E+00  0.57426986E+00  0.66452897E+00  0.74420635E+00  0.82147910E+00  0.90630485E+00
  97.0000000  0.72920362E-01  0.15820555E+00  0.25905105E+00  0.36813847E+00  0.47444314E+00  0.57137750E+00  0.65956773E+00  0.74367914E+00  0.82779366E+00  0.91337704E+00
  98.0000000  0.74172639E-01  0.15942721E+00  0.25861565E+00  0.36443530E+00  0.46629133E+00  0.56012249E+00  0.65058868E+00  0.74297601E+00  0.83531268E+00  0.92130034E+00
  99.0000000  0.85725350E-01  0.17363570E+00  0.26451836E+00  0.35737690E+00  0.45091524E+00  0.54557561E+00  0.64289670E+00  0.74182751E+00  0.83691982E+00  0.92262727E+00
 100.0000000  0.99265650E-01  0.19056024E+00  0.27233912E+00  0.35123543E+00  0.43687255E+00  0.53347052E+00  0.63632325E+00  0.73675504E+00  0.82977037E+00  0.91635393E+00

https://ideone.com/ZXQY49

language:

Fortran (gfortran 8.3)

created:

visibility:

public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!

Discover > Sphere Engine API

Discover > IDE Widget

Choose your language