*! xtfpanic 1.0.2  23jul2026   (xtflexur library)
*! Fourier-PANIC panel unit root test: PANIC with smooth (Fourier) structural
*! breaks and a common factor structure.
*!
*! Author : Dr Merwan Roudane  merwanroudane920@gmail.com
*!          https://github.com/merwanroudane
*! Part of the -xtflexur- library : factor-augmented panel time-series tests.
*!
*! Faithful port of the GAUSS routine in appl3_PANIC_fourier.gss by
*!   Saban Nazlioglu et al.  Non-commercial use.
*!
*! References:
*!   Nazlioglu, S., et al. (2023). Smooth structural changes and common factors
*!     in nonstationary panel data: an analysis of healthcare expenditures.
*!     Econometric Reviews 42(1): 78-97. <doi:10.1080/07474938.2022.2156740>
*!   Bai, J., Ng, S. (2004). A PANIC attack on unit roots and cointegration.
*!     Econometrica 72: 1127-1177.
*!   Enders, W., Lee, J. (2012). The flexible Fourier form and Dickey-Fuller type
*!     unit root tests. Economics Letters 117: 196-199.

program xtfpanic, rclass
    version 14.0

    syntax varname(numeric ts) [if] [in] [ ,   ///
        Freq(integer 1)                         ///
        Kmax(integer 5)                         ///
        ICfactor(integer 2)                     ///
        Pmax(integer 3)                         ///
        IClag(integer 3)                        ///
    ]

    if `freq' < 1 | `freq' > 3 {
        di as error "xtfpanic: freq() must be 1, 2 or 3 (max cumulative frequency)."
        exit 198
    }

    quietly xtset
    local pv "`r(panelvar)'"
    local tv "`r(timevar)'"
    if "`pv'" == "" | "`tv'" == "" {
        di as error "xtfpanic: data must be {help xtset}."
        exit 459
    }
    if "`r(balanced)'" != "strongly balanced" {
        di as error "xtfpanic: requires a strongly balanced panel."
        exit 459
    }

    marksample touse
    markout `touse' `pv' `tv'
    quietly count if `touse'
    local NT = r(N)
    quietly levelsof `pv' if `touse', local(ids)
    local N : word count `ids'
    local T = `NT'/`N'
    if `T' != int(`T') | `N'*`T' != `NT' {
        di as error "xtfpanic: unbalanced over the selected sample."
        exit 459
    }
    if `T' < 12 {
        di as error "xtfpanic: need at least 12 time periods (found `T')."
        exit 2001
    }

    tempname R OUT
    mata: _xtfp_run("`varlist'", "`touse'", `N', `T', `freq', `kmax', `icfactor', `pmax', `iclag')
    matrix `R'   = r(pooled)
    matrix `OUT' = r(units)

    * ----- display -----------------------------------------------------------
    local icfd = cond(`icfactor'==1,"PCp",cond(`icfactor'==2,"ICp","AIC/BIC"))
    local mstr = cond(`freq'==1,"1",cond(`freq'==2,"1,2","1,2,3"))
    di ""
    di as text "Fourier-PANIC panel unit root test (PANIC with smooth breaks)"
    di as text "Nazlioglu et al. (2023, Econometric Reviews)"
    di as text "{hline 62}"
    di as text "Panel var" _col(14) ": " as result "`pv'" ///
        as text _col(42) "N" _col(50) "= " as result %5.0f `N'
    di as text "Time var" _col(14) ": " as result "`tv'" ///
        as text _col(42) "T" _col(50) "= " as result %5.0f `T'
    di as text "Frequencies" _col(14) ": " as result "m = `mstr'" ///
        as text _col(42) "Factors" _col(50) "= " as result %3.0f `R'[3,1] ///
        as text " (`icfd')"
    di as text "H0: idiosyncratic component has a unit root (non-stationary)"
    di as text "{hline 62}"
    di as text "  Pooled test" _col(26) "Statistic" _col(42) "p-value"
    di as text "{hline 62}"
    local P = `R'[1,1]
    local Pp = `R'[1,2]
    local Pm = `R'[2,1]
    local Pmp = `R'[2,2]
    local s1 = cond(`Pp'<0.01,"***",cond(`Pp'<0.05,"**",cond(`Pp'<0.10,"*","")))
    local s2 = cond(`Pmp'<0.01,"***",cond(`Pmp'<0.05,"**",cond(`Pmp'<0.10,"*","")))
    di as text "  P  (Fisher, chi2)" _col(26) as result %9.3f `P' _col(42) %7.3f `Pp' "  `s1'"
    di as text "  Pm (standardized)" _col(26) as result %9.3f `Pm' _col(42) %7.3f `Pmp' "  `s2'"
    di as text "{hline 62}"
    di as text "  Per-unit Fourier-LM (idiosyncratic):"
    di as text "   id" _col(12) "LM" _col(26) "p-value" _col(38) "lags"
    forvalues i = 1/`N' {
        di as text "  " %4.0f `OUT'[`i',1] _col(9) as result %9.3f `OUT'[`i',2] ///
            _col(26) %7.3f `OUT'[`i',3] _col(39) %3.0f `OUT'[`i',4]
    }
    di as text "{hline 62}"
    di as text "Note: reject H0 (=> stationary idiosyncratic component) for large"
    di as text "      P/Pm.  *** p<.01, ** p<.05, * p<.10."

    * ----- returns -----------------------------------------------------------
    return scalar P    = `P'
    return scalar P_p  = `Pp'
    return scalar Pm   = `Pm'
    return scalar Pm_p = `Pmp'
    return scalar nf   = `R'[3,1]
    return scalar N    = `N'
    return scalar T    = `T'
    return matrix units = `OUT'
    return local  cmd  "xtfpanic"
end


* =========================================================================
* Mata engine
* =========================================================================
version 14.0
mata:

real scalar _xtfp_argmin(real colvector v)
{
    real scalar i, mn, mi
    mn = v[1]
    mi = 1
    i = 2
    while (i <= rows(v)) {
        if (v[i] < mn) {
            mn = v[i]
            mi = i
        }
        i = i + 1
    }
    return(mi)
}

real colvector _xtfp_lagn_k(real colvector x, real scalar k)
{
    real scalar t
    t = rows(x)
    return((J(k,1,0)) \ x[1..t-k])
}

void _xtfp_ols(real colvector y, real matrix X, real colvector b, real colvector se, real colvector e)
{
    real matrix m
    real scalar sig2
    m = luinv(cross(X,X))
    b = m*cross(X,y)
    e = y - X*b
    sig2 = cross(e,e)/(rows(y)-cols(X))
    se = sqrt(diagonal(m)*sig2)
}

// ---- principal components (shared factor engine) ----------------------
void _xtfp_pca(real matrix X, real scalar nf, real matrix Fhat, real matrix lambda, real matrix ehat)
{
    real scalar t, n
    real matrix evec, A
    real rowvector eval
    real colvector ord
    t = rows(X)
    n = cols(X)
    if (t < n) {
        A = X*X'
        symeigensystem(A, evec, eval)
        ord = order(eval', -1)
        Fhat = sqrt(t) * evec[., ord[1..nf]]
        lambda = luinv(cross(Fhat,Fhat)) * cross(Fhat,X)
        ehat = X - Fhat*lambda
    }
    else {
        A = X'X
        symeigensystem(A, evec, eval)
        ord = order(eval', -1)
        lambda = sqrt(n) * evec[., ord[1..nf]]
        Fhat = X*lambda/n
        ehat = X - Fhat*lambda'
    }
}

real scalar _xtfp_fnumber2(real matrix X, real scalar kmax, real scalar ic)
{
    real scalar t, n, sigma2, k, mn
    real matrix Fhat, lambda, ehat, ICp
    real colvector V
    t = rows(X)
    n = cols(X)
    _xtfp_pca(X, kmax, Fhat, lambda, ehat)
    sigma2 = sum(ehat:^2)/(n*t)
    mn = min((n,t))
    V = J(kmax,1,0)
    ICp = J(kmax,1,0)
    k = 1
    while (k <= kmax) {
        _xtfp_pca(X, k, Fhat, lambda, ehat)
        V[k] = sum(ehat:^2)/(n*t)
        if (ic == 1)      ICp[k] = V[k] + k*sigma2*(n+t)/(n*t)*ln(mn)
        else if (ic == 2) ICp[k] = ln(V[k]) + k*(n+t)/(n*t)*ln(mn)
        else              ICp[k] = V[k] + k*sigma2*(ln(t)/t)
        k = k + 1
    }
    return(_xtfp_argmin(ICp))
}

// ---- 221-percentile grid used by the response surface -----------------
real colvector _xtfp_pgrid()
{
    real colvector p
    p = (0.0001 \ 0.0002 \ 0.0005) \ ((1::9)/1000) \ (0.005:*(2::198)) \ ((991::999)/1000) \ (0.9995 \ 0.9998 \ 0.9999)
    return(p)
}

real matrix _xtfp_rsblock(real scalar m)
{
    if (m == 0) return(_xtfp_rs0())
    if (m == 1) return(_xtfp_rs1())
    if (m == 2) return(_xtfp_rs2())
    return(_xtfp_rs3())
}

// ---- response-surface p-value (MacKinnon interpolation) ---------------
real scalar _xtfp_pval(real scalar test, real scalar T, real scalar m, real scalar plags)
{
    real matrix params, Xw
    real colvector p, reg, cv, yv, bw, idx
    real scalar it, pt, mn, mx, bands, lower, lo, hi, pval
    params = _xtfp_rsblock(m)
    it = 1/T
    pt = plags/T
    reg = (1 \ it \ it^2 \ pt \ pt^2)
    cv = params*reg
    p = _xtfp_pgrid()
    mn = min(cv)
    mx = max(cv)
    if (test < mn) return(0.0001)
    if (test > mx) return(1)
    bands = 11
    idx = selectindex(cv :<= test)
    lower = max(idx)
    lo = lower - bands + 1
    hi = lower + bands
    if (lower < bands) {
        lo = 1
        hi = bands
    }
    if (lower > rows(cv)-bands) {
        lo = rows(cv) - bands - 1
        hi = rows(cv)
    }
    yv = invnormal(p[lo..hi])
    Xw = J(hi-lo+1,1,1) , cv[lo..hi] , cv[lo..hi]:^2 , cv[lo..hi]:^3
    bw = qrsolve(Xw, yv)
    pval = normal((1, test, test^2, test^3) * bw)
    if (pval > 1) pval = 1
    if (pval < 0) pval = 0.0001
    return(pval)
}

// ---- Fourier-LM statistic on one unit's idiosyncratic component -------
real rowvector _xtfp_lmtau(real colvector dyi, real colvector wwi, real matrix dz, real matrix Fhat, real scalar pmax, real scalar icp)
{
    real colvector sti, dsti, s1, aic, sic, lmt, lmp, b, se, e, dyitr
    real matrix ind, indk, indtr, lmat
    real scalar k, j, laglm, tau, r
    sti = runningsum(wwi)
    dsti = (0 \ (sti[2..rows(sti)] - sti[1..rows(sti)-1]))
    s1 = (0 \ sti[1..rows(sti)-1])
    ind = s1 , dz , Fhat
    aic = J(pmax+1,1,0)
    sic = J(pmax+1,1,0)
    lmt = J(pmax+1,1,0)
    lmp = J(pmax+1,1,0)
    k = 0
    while (k <= pmax) {
        if (k > 0) {
            lmat = J(rows(dsti), k, 0)
            j = 1
            while (j <= k) {
                lmat[.,j] = _xtfp_lagn_k(dsti, j)
                j = j + 1
            }
            indk = ind , lmat[., 1..k]
        }
        else indk = ind
        dyitr = dyi[(k+1)..rows(dyi)]
        indtr = indk[(k+1)..rows(indk), .]
        _xtfp_ols(dyitr, indtr, b, se, e)
        r = rows(indtr)
        aic[k+1] = ln(cross(e,e)/r) + 2*(cols(indtr)+2)/r
        sic[k+1] = ln(cross(e,e)/r) + (cols(indtr)+2)*ln(r)/r
        lmt[k+1] = b[1]/se[1]
        lmp[k+1] = b[cols(indtr)]/se[cols(indtr)]
        k = k + 1
    }
    if (icp == 1) laglm = _xtfp_argmin(aic)
    else if (icp == 2) laglm = _xtfp_argmin(sic)
    else {
        laglm = 1
        j = pmax + 1
        while (j >= 1) {
            if (abs(lmp[j]) > 1.645 | j == 1) {
                laglm = j
                break
            }
            j = j - 1
        }
    }
    tau = lmt[laglm]
    return((tau, laglm-1))
}

// ---- driver ------------------------------------------------------------
void _xtfp_run(string scalar vname, string scalar touse, real scalar N, real scalar T, real scalar m, real scalar kmax, real scalar ick, real scalar pmax, real scalar icp)
{
    real colvector yv, ttt, dyi, tau, lags, pval
    real matrix Y, dy, z, dz, dyh, rr, Fhat, lambda, ehat, ww, ccc, sss, units, pooled
    real scalar nf, i, kf, P, Pm, Tm1
    real colvector tdel
    real rowvector lr

    yv = st_data(., vname, touse)
    Y = rowshape(yv, N)'              // T x N
    dy = Y[2..T,.] - Y[1..T-1,.]      // (T-1) x N
    Tm1 = rows(dy)
    ttt = (1::T)

    // Fourier terms z = [t, cos(2 pi k t/T), sin(2 pi k t/T)], k=1..m
    ccc = J(T, m, 0)
    sss = J(T, m, 0)
    kf = 1
    while (kf <= m) {
        ccc[.,kf] = cos(2*pi()*kf*ttt/T)
        sss[.,kf] = sin(2*pi()*kf*ttt/T)
        kf = kf + 1
    }
    z = ttt , ccc , sss
    dz = z[2..T,.] - z[1..T-1,.]      // (T-1) x (1+2m)

    // Fourier-detrend each differenced series -> residuals rr
    dyh = J(Tm1, N, 0)
    rr  = J(Tm1, N, 0)
    i = 1
    while (i <= N) {
        tdel = qrsolve(dz, dy[.,i])
        dyh[.,i] = dz*tdel
        rr[.,i]  = dy[.,i] - dyh[.,i]
        i = i + 1
    }

    // factor extraction on the Fourier-detrended differences
    nf = _xtfp_fnumber2(rr, kmax, ick)
    _xtfp_pca(rr, nf, Fhat, lambda, ehat)
    if (T > N) ww = dy - dyh - Fhat*lambda'
    else       ww = dy - dyh - Fhat*lambda

    tau  = J(N,1,.)
    lags = J(N,1,.)
    pval = J(N,1,.)
    i = 1
    while (i <= N) {
        lr = _xtfp_lmtau(dy[.,i], ww[.,i], dz, Fhat, pmax, icp)
        tau[i]  = lr[1]
        lags[i] = lr[2]
        pval[i] = _xtfp_pval(tau[i], Tm1, m, lags[i])
        i = i + 1
    }

    P  = -2*sum(ln(pval))
    Pm = (P - 2*N) / sqrt(4*N)

    units = (1::N) , tau , pval , lags
    pooled = (P , chi2tail(2*N, P)) \ (Pm , 1-normal(Pm)) \ (nf , .)
    st_matrix("r(units)", units)
    st_matrix("r(pooled)", pooled)
}


// ---- Fourier-PANIC LM response-surface tables (221x5, max_m=0..3) ----
// CV(p,T,lags) = c1 + c2/T + c3/T^2 + c4*(lags/T) + c5*(lags/T)^2 (appl3 rsf)
real matrix _xtfp_rs0()
{
    string scalar s
    s = "-4.65260 -65.17050 1036.22190 4.46840 -5.35230 -4.52400 -49.49360 643.06510 4.79310 -11.16730 -4.32830 -15.67680 -138.33440 2.15190 -3.74360 -4.19070 -16.56060 -65.27290 2.74530 -4.40940 -4.06460 -6.80310 -213.83640 3.32120 -10.67230 -3.94300 -1.80570 -249.84290 2.45660 -7.95110 -3.86970 -2.25530 -191.39410 2.61580 -8.63060 -3.80520 -5.92910 -77.97180 3.04830 -9.85260 -3.74740 -3.73120 -128.50440 2.59140 -7.36700 -3.69770 -2.05790 -165.78010 2.29920 -5.58080 -3.66050 -1.36770 -163.58120 2.25420 -5.46160 -3.62570 -0.43210 -187.37960 2.24960 -5.65500 -3.59390 -1.54460 -146.13400 2.35080 -5.68100 -3.45110 0.21610 -157.66790 1.76250 -2.68260 -3.35260 -0.94620 -103.59880 1.70550 -1.83760 -3.26800 -0.88130 -79.94900 1.45080 -0.58330 -3.21040 1.45950 -128.64150 1.35750 -0.66250 -3.15140 1.09680 -110.38610 1.26130 -0.00990 -3.10530 1.08000 -103.62300 1.27060 0.17260 -3.06030 1.28970 -105.07730 1.16730 0.76030 -3.01740 1.80290 -109.56360 1.01940 1.32180 -2.98190 1.19680 -90.93000 1.10910 1.00570 -2.94840 2.04410 -107.82200 0.97550 1.55530 -2.91640 2.32210 -105.74160 0.88810 2.00590 -2.88890 2.77920 -113.37390 0.92010 1.64020 -2.86080 2.64100 -101.73560 0.92700 1.53770 -2.83800 2.76820 -101.29740 0.97290 1.32410 -2.81430 2.48080 -89.59800 1.01190 1.16570 -2.79260 2.83930 -97.99960 1.02050 1.09850 -2.77160 2.73310 -93.18390 1.03760 1.11330 "
    s = s + "-2.75010 3.15930 -97.34440 0.98830 1.14510 -2.72800 3.12830 -93.90600 0.95340 1.27300 -2.70810 2.78140 -79.97610 0.97140 1.18200 -2.68970 2.21650 -62.49010 1.01180 1.20210 -2.67210 2.24360 -61.35430 1.02790 1.10530 -2.65470 2.10310 -53.34660 1.04830 0.98970 -2.63630 2.54250 -61.63600 0.99660 1.05320 -2.61920 2.61890 -62.33470 0.95660 1.26340 -2.60490 2.40680 -55.50160 1.03960 0.84240 -2.58860 2.71300 -59.35360 0.99840 0.91840 -2.57340 2.71000 -58.37660 1.01060 0.83470 -2.55980 2.60210 -54.69880 1.04570 0.70210 -2.54540 2.68350 -56.17830 1.04040 0.72300 -2.52800 2.79080 -56.79390 0.94570 1.15960 -2.51470 2.51670 -48.37870 0.97880 1.10630 -2.50220 2.31400 -42.33560 1.02960 0.89880 -2.48950 2.37940 -44.33330 1.04520 0.80830 -2.47590 2.20170 -39.52080 1.05720 0.77030 -2.46280 2.25810 -39.10990 1.03530 0.87540 -2.44990 1.96430 -31.63720 1.05300 0.86690 -2.43740 1.84050 -28.59720 1.06220 0.84750 -2.42540 2.13570 -34.87240 1.02930 0.93690 -2.41240 2.30630 -38.60260 0.99520 1.04450 -2.40040 2.31390 -39.13260 0.98180 1.12630 -2.38950 2.35590 -38.93400 0.99470 1.02950 -2.37770 2.64080 -46.36480 0.96360 1.10060 -2.36610 2.86510 -49.76660 0.93140 1.17890 -2.35500 2.70410 -44.01690 0.95230 1.07650 -2.34380 2.86670 -45.66140 0.93120 1.08380 -2.33370 2.64850 -38.86880 0.95930 0.98000 "
    s = s + "-2.32250 2.60380 -38.02160 0.95300 0.98840 -2.31250 2.65030 -38.38260 0.95650 0.93770 -2.30190 2.97200 -45.52360 0.92180 1.01530 -2.29050 2.88820 -42.82950 0.90830 1.07080 -2.28130 2.80420 -40.64820 0.93510 0.94270 -2.27090 2.74520 -38.31900 0.92150 1.02460 -2.26050 2.78960 -38.33230 0.90510 1.06450 -2.25080 3.09630 -45.42330 0.88540 1.03940 -2.24170 2.95990 -40.85120 0.92060 0.82990 -2.23200 3.03110 -42.56160 0.90620 0.87230 -2.22310 2.82510 -36.92680 0.93800 0.74930 -2.21340 2.67880 -32.69790 0.94880 0.71210 -2.20480 2.76840 -33.90610 0.95310 0.65250 -2.19500 2.84270 -34.67880 0.93750 0.66070 -2.18610 2.65210 -29.31460 0.96470 0.53870 -2.17660 2.93400 -34.92780 0.92930 0.60520 -2.16710 2.91160 -33.26450 0.90570 0.73100 -2.15830 2.88380 -32.10760 0.90770 0.68540 -2.15010 2.83980 -30.23980 0.92470 0.56600 -2.14190 2.69790 -25.80150 0.94000 0.51480 -2.13310 2.64760 -22.91920 0.93420 0.52470 -2.12470 2.69280 -23.88350 0.93290 0.51750 -2.11650 2.67540 -23.94280 0.93370 0.52730 -2.10760 2.73430 -24.20210 0.90110 0.68970 -2.09930 2.66360 -22.09600 0.90910 0.64720 -2.09110 2.42880 -15.43280 0.92600 0.59350 -2.08220 2.52020 -17.74790 0.89700 0.72420 -2.07300 2.78860 -24.41520 0.83940 0.95480 -2.06490 2.74250 -23.62450 0.84190 0.96760 -2.05650 2.75220 -22.87050 0.82660 1.00660 "
    s = s + "-2.04870 2.67030 -19.48280 0.83290 0.97740 -2.04070 2.79810 -23.41100 0.81870 1.00310 -2.03250 2.81750 -24.38180 0.80970 1.02730 -2.02480 2.74820 -22.23840 0.81710 0.99710 -2.01760 2.81200 -23.59430 0.82970 0.86900 -2.00960 2.72070 -21.71700 0.82310 0.93320 -2.00260 2.72790 -21.75200 0.83770 0.83930 -1.99460 2.73180 -21.90530 0.82370 0.89250 -1.98640 2.67990 -20.68980 0.81370 0.96100 -1.97860 2.80360 -23.57130 0.79240 1.03610 -1.97110 2.63920 -19.44090 0.81300 0.96180 -1.96370 2.60840 -18.29310 0.81420 0.96610 -1.95650 2.60460 -17.56930 0.80520 1.02780 -1.94840 2.73260 -20.08400 0.77450 1.13930 -1.94050 2.74850 -20.12290 0.75420 1.23790 -1.93300 2.63030 -16.45550 0.75770 1.21980 -1.92510 2.63830 -16.61320 0.73900 1.31990 -1.91780 2.68890 -17.89220 0.73170 1.33170 -1.90930 2.59590 -15.54820 0.69760 1.55310 -1.90180 2.55340 -14.27810 0.69560 1.57230 -1.89490 2.34290 -8.84650 0.72590 1.44560 -1.88770 2.36360 -9.77060 0.72290 1.44100 -1.88030 2.32900 -8.02510 0.72400 1.40510 -1.87210 2.20890 -5.07190 0.69990 1.58080 -1.86460 2.28040 -6.47980 0.67280 1.71020 -1.85820 2.02820 0.02140 0.71050 1.57730 -1.85170 2.01960 0.49060 0.72580 1.48380 -1.84370 2.10700 -1.11440 0.69050 1.63210 -1.83610 2.06630 0.89280 0.67000 1.76700 -1.82840 2.08760 0.08920 0.65220 1.85280 "
    s = s + "-1.82080 2.10240 -1.06920 0.62660 2.02680 -1.81350 2.05860 0.51030 0.62490 2.03090 -1.80630 2.02380 1.28080 0.62190 2.05660 -1.79860 2.19770 -2.00900 0.58250 2.19490 -1.79110 2.11730 -0.42910 0.57560 2.27630 -1.78350 2.13890 -0.66640 0.55180 2.37820 -1.77650 2.19290 -1.39370 0.54070 2.42300 -1.76850 2.28770 -3.54670 0.49990 2.60750 -1.76180 2.34240 -5.23270 0.49650 2.63670 -1.75460 2.39020 -7.07310 0.49570 2.61420 -1.74720 2.23440 -2.86410 0.49460 2.66200 -1.74030 2.19470 -2.58430 0.49260 2.71540 -1.73330 2.15410 -1.54730 0.49210 2.74240 -1.72640 2.15860 -2.28780 0.48940 2.77210 -1.72000 2.16380 -2.46530 0.50370 2.67730 -1.71260 2.14390 -1.24890 0.48920 2.75580 -1.70570 2.02960 2.14370 0.49860 2.71830 -1.69870 1.87250 5.23940 0.51940 2.65460 -1.69140 1.84340 6.47290 0.50650 2.72040 -1.68460 1.80610 7.20920 0.50050 2.79980 -1.67700 1.83290 6.81870 0.47170 2.96540 -1.67030 1.92460 4.85610 0.45510 3.06220 -1.66240 1.80450 8.79820 0.44630 3.12380 -1.65520 1.75560 11.04680 0.43830 3.17470 -1.64790 1.83250 9.59920 0.41480 3.31790 -1.64100 1.72990 12.20280 0.41200 3.40440 -1.63330 1.75960 12.48500 0.38970 3.51320 -1.62630 1.56330 16.65390 0.40540 3.52330 -1.61860 1.54550 18.05100 0.38600 3.63600 -1.61080 1.50610 18.07590 0.37020 3.75020 "
    s = s + "-1.60310 1.43280 19.41560 0.35550 3.89440 -1.59630 1.47000 18.09870 0.35060 3.97170 -1.58960 1.34720 21.33630 0.36870 3.93720 -1.58270 1.14340 25.98640 0.39190 3.91660 -1.57430 1.11600 26.76300 0.36250 4.09720 -1.56750 0.94180 32.03140 0.37410 4.11720 -1.56000 0.86110 34.19250 0.36050 4.27750 -1.55350 0.95370 31.27350 0.37470 4.20540 -1.54580 1.00010 30.29620 0.35330 4.35080 -1.53830 1.03330 29.52970 0.33360 4.50860 -1.53010 1.09790 28.25160 0.30730 4.64760 -1.52240 1.15840 26.92340 0.29350 4.73600 -1.51510 1.18160 26.32030 0.28670 4.80890 -1.50800 1.04510 29.38400 0.29790 4.88340 -1.50080 0.96820 30.74030 0.30530 4.95400 -1.49310 1.01230 29.11810 0.29510 5.08760 -1.48570 1.05280 27.93410 0.28450 5.21470 -1.47860 0.96890 30.33630 0.29370 5.27550 -1.46980 0.77080 35.49760 0.27810 5.52140 -1.46200 0.63960 39.28610 0.28580 5.58780 -1.45390 0.56490 41.82480 0.28380 5.69160 -1.44500 0.66960 39.25920 0.24720 5.94570 -1.43700 0.76400 37.62750 0.22420 6.14210 -1.42930 0.80340 37.05150 0.22610 6.23810 -1.42090 0.80000 37.30210 0.21620 6.40380 -1.41200 0.93910 34.11390 0.18140 6.67390 -1.40390 1.13210 30.19040 0.16100 6.83740 -1.39460 1.41210 23.69360 0.10170 7.23460 -1.38580 1.42190 23.73200 0.08830 7.44020 -1.37710 1.25320 28.66710 0.11130 7.47590 "
    s = s + "-1.36810 1.14760 32.76440 0.10770 7.67830 -1.35950 1.32090 29.48130 0.09520 7.86440 -1.35070 0.95440 39.08500 0.13450 7.95480 -1.34160 0.73940 44.03780 0.16670 8.03610 -1.33290 0.78540 42.42980 0.17600 8.19100 -1.32480 0.51830 49.35270 0.23820 8.19100 -1.31460 0.33760 53.87670 0.23290 8.56590 -1.30420 0.66530 46.71990 0.17030 9.12960 -1.29490 0.76730 44.30300 0.19440 9.22610 -1.28460 0.73490 45.54030 0.19610 9.59380 -1.27360 0.89300 43.75790 0.17700 9.94570 -1.26280 0.74390 47.26530 0.20780 10.18820 -1.25340 0.70350 47.94640 0.26610 10.30230 -1.24210 0.98040 42.78550 0.24840 10.72480 -1.23010 0.93390 42.63610 0.27610 11.07550 -1.21800 0.75900 46.22370 0.32660 11.36350 -1.20610 1.13060 38.35760 0.33980 11.70130 -1.19470 1.06250 40.80160 0.43260 11.75040 -1.18280 1.20780 38.18340 0.51790 11.81080 -1.17090 1.16590 37.10380 0.63780 12.02810 -1.15740 1.63750 26.25500 0.70240 12.24910 -1.14270 1.67890 27.86430 0.80740 12.43780 -1.12920 1.54720 32.24480 1.02330 12.14170 -1.10910 1.84450 24.42070 1.07800 13.01300 -1.09350 1.46520 32.77000 1.43050 12.50600 -1.07540 1.77830 19.33720 1.74630 12.33490 -1.04660 3.35740 -19.39750 1.78290 13.55350 -1.02480 0.63580 51.20150 2.72250 11.49930 -0.99480 1.56540 31.95390 3.51970 10.01300 -0.98230 2.00400 24.24100 3.54230 10.49460 "
    s = s + "-0.97450 2.06340 22.60440 3.83320 9.84530 -0.96700 -0.04000 72.97870 4.46130 8.02930 -0.95540 -0.87670 97.57280 4.85640 7.12620 -0.94070 -0.45480 91.22840 5.13020 6.70380 -0.92140 1.78890 32.46530 5.28960 6.74180 -0.90760 3.36700 6.94740 6.06250 3.72830 -0.88460 6.28100 -89.37180 6.89310 1.14560 -0.82950 6.48000 -72.64390 7.88910 0.28370 -0.79430 11.17960 -149.30740 8.75320 -1.51890 -0.72760 16.48860 -196.61170 9.03360 0.32720 -0.68640 26.05120 -355.19170 8.70630 2.69940 "
    return(rowshape(strtoreal(tokens(s)), 221))
}

real matrix _xtfp_rs1()
{
    string scalar s
    s = "-5.92890 -62.19210 -123.83560 9.40040 -25.16750 -5.59210 21.54660 -2176.75900 1.36070 -11.01080 -5.31080 10.28280 -1447.54650 -0.43480 1.87620 -5.12500 12.31230 -1392.61870 -1.66080 9.80940 -4.97490 0.32540 -881.59330 0.66240 -0.29210 -4.87310 1.51640 -768.49070 0.79680 -2.27350 -4.80150 -4.78540 -547.73520 1.31690 -2.31960 -4.73340 -7.00160 -427.85390 1.28330 -1.37310 -4.67110 -9.05650 -327.25710 1.17700 0.13290 -4.63360 -8.75050 -310.88100 1.31860 -0.73300 -4.59900 -6.14250 -359.44420 1.12400 0.04310 -4.56750 -5.87180 -343.77760 1.19560 -0.36480 -4.53610 -6.31970 -316.75900 1.28930 -0.68810 -4.42760 -6.19420 -255.76120 1.49870 -1.03520 -4.34750 -4.90070 -246.38620 1.58480 -1.40490 -4.27410 -4.77040 -212.05380 1.62290 -1.61640 -4.20790 -4.05910 -200.55920 1.37830 -0.33890 -4.15700 -4.34250 -168.41420 1.40420 -0.19920 -4.10770 -4.59010 -140.68510 1.32430 0.54140 -4.06960 -4.11420 -135.96900 1.33740 0.53370 -4.03260 -4.37170 -107.46270 1.34510 0.52390 -3.99620 -4.60610 -86.04100 1.25190 1.26080 -3.96240 -3.82490 -89.49930 1.12580 1.66460 -3.92820 -3.74570 -81.06360 1.00320 2.32010 -3.90140 -3.50730 -77.52140 1.01770 2.24040 -3.87130 -3.10170 -77.62010 0.85830 3.08110 -3.84640 -3.44300 -57.81250 0.86160 3.19790 -3.82310 -3.13390 -53.70040 0.84800 3.21140 -3.79820 -2.92560 -52.69700 0.74200 3.81100 -3.77830 -2.79280 -49.22980 0.75420 3.78580 "
    s = s + "-3.75760 -2.79190 -40.63420 0.75430 3.83690 -3.73750 -2.85280 -31.96900 0.75140 3.86400 -3.71960 -2.91030 -25.11730 0.79410 3.69180 -3.70180 -2.55720 -31.38520 0.78160 3.73160 -3.68150 -2.41990 -28.79430 0.69810 4.16600 -3.66500 -2.26400 -29.06900 0.70290 4.11750 -3.64930 -2.67210 -13.69170 0.75030 4.08150 -3.63310 -2.74890 -6.16700 0.75760 4.04510 -3.61760 -2.31480 -11.09390 0.73140 4.03940 -3.60280 -2.12320 -13.34690 0.73200 4.03100 -3.58790 -1.87480 -14.69060 0.70000 4.17810 -3.57490 -1.64200 -18.32180 0.72500 4.00030 -3.56130 -1.27170 -26.10890 0.73140 3.86260 -3.54540 -1.08080 -27.50900 0.66360 4.16940 -3.53030 -1.03120 -26.40920 0.61630 4.44260 -3.51690 -1.13590 -20.37310 0.61840 4.48520 -3.50290 -1.01270 -20.73940 0.59080 4.58040 -3.49130 -1.03470 -16.65990 0.59660 4.65290 -3.47860 -0.80810 -20.61880 0.57630 4.72290 -3.46520 -1.05920 -10.12630 0.56540 4.89380 -3.45210 -1.10140 -5.32160 0.55010 4.95720 -3.43820 -0.82340 -9.05630 0.47050 5.31860 -3.42540 -1.01920 -1.37840 0.46090 5.42660 -3.41310 -1.20490 6.99130 0.46450 5.41860 -3.40240 -1.19290 7.98770 0.48690 5.27900 -3.39090 -0.97730 4.97870 0.45120 5.42610 -3.37920 -0.89030 4.50410 0.42290 5.59380 -3.36860 -0.67340 0.61540 0.39520 5.71850 -3.35800 -0.71960 3.67710 0.40110 5.72910 -3.34760 -0.34300 -2.98610 0.36450 5.82580 "
    s = s + "-3.33630 -0.08100 -8.47310 0.31650 6.03230 -3.32540 0.12890 -11.28310 0.27850 6.18120 -3.31510 0.27540 -13.61550 0.25940 6.26160 -3.30570 0.47240 -17.04180 0.26730 6.16900 -3.29560 0.43200 -13.07340 0.27410 6.13100 -3.28580 0.52900 -13.78300 0.26290 6.19310 -3.27610 0.24690 -4.93700 0.28510 6.20200 -3.26690 0.37980 -5.57600 0.27570 6.22730 -3.25680 0.27090 -1.01420 0.27400 6.27790 -3.24690 0.26600 1.07360 0.24840 6.46300 -3.23770 0.14090 5.33450 0.27140 6.36130 -3.22770 0.24600 4.29280 0.23830 6.51190 -3.21860 0.12360 9.63410 0.24870 6.48430 -3.20980 0.16490 8.71470 0.24760 6.51430 -3.19940 0.25710 8.84050 0.19200 6.78120 -3.18960 0.25790 10.89770 0.17030 6.88710 -3.18000 0.51540 6.97340 0.12890 7.00330 -3.17150 0.50530 7.80240 0.13830 6.96380 -3.16260 0.67830 5.83770 0.11840 6.99050 -3.15330 0.70210 6.26070 0.10140 7.09360 -3.14520 0.75660 6.37610 0.11160 7.01950 -3.13580 0.85420 5.07390 0.07660 7.17670 -3.12660 0.80890 8.25640 0.06110 7.27080 -3.11770 0.87340 6.65370 0.03610 7.44400 -3.10900 0.70170 12.92710 0.04070 7.47020 -3.10020 0.88820 10.44630 0.00230 7.63570 -3.09140 0.73780 16.10050 0.00250 7.66410 -3.08220 0.84850 14.73580 -0.03480 7.85010 -3.07380 0.79010 18.18120 -0.05510 8.00470 -3.06580 0.74730 20.58880 -0.03730 7.91010 "
    s = s + "-3.05750 0.70280 22.06920 -0.04000 7.97650 -3.04970 0.73150 22.47200 -0.04700 8.02090 -3.04110 0.48980 29.51050 -0.02830 7.96880 -3.03300 0.44570 31.39210 -0.02650 7.98270 -3.02530 0.50490 30.61180 -0.02180 7.94030 -3.01750 0.53570 32.16300 -0.04000 8.05210 -3.00960 0.51640 34.43220 -0.03800 8.04140 -3.00050 0.67280 31.36860 -0.07840 8.19970 -2.99290 0.70130 31.76820 -0.06860 8.11870 -2.98430 0.44640 39.69120 -0.07870 8.28900 -2.97620 0.41410 42.94640 -0.09370 8.38570 -2.96740 0.46440 43.74730 -0.13870 8.64210 -2.96020 0.50310 43.70830 -0.13850 8.64310 -2.95160 0.41170 47.38390 -0.15020 8.71410 -2.94350 0.42010 47.81920 -0.16150 8.80000 -2.93520 0.44640 48.73140 -0.18420 8.90800 -2.92830 0.29960 53.15990 -0.15440 8.80890 -2.92050 0.15400 57.04890 -0.14780 8.85960 -2.91240 -0.02430 63.25340 -0.14730 8.91170 -2.90390 0.13550 61.29280 -0.19500 9.12870 -2.89670 0.08590 63.87870 -0.18020 9.05790 -2.88830 0.08310 63.90980 -0.20270 9.20080 -2.88010 0.30870 58.53730 -0.24380 9.38450 -2.87260 0.27390 60.37100 -0.25380 9.48680 -2.86490 0.30800 60.57020 -0.25720 9.49300 -2.85670 0.36990 59.18520 -0.28750 9.68700 -2.84980 0.31000 61.77040 -0.27190 9.63000 -2.84300 0.18780 65.54460 -0.25280 9.60120 -2.83460 0.21360 65.87730 -0.27140 9.69450 -2.82730 0.02850 71.83960 -0.24860 9.61860 "
    s = s + "-2.81970 0.09200 70.76260 -0.27370 9.79790 -2.81130 0.03590 72.01480 -0.29450 9.94480 -2.80360 0.05760 73.76560 -0.31280 10.04610 -2.79600 0.02190 74.33040 -0.30830 10.08060 -2.78820 -0.03140 75.03260 -0.31260 10.16350 -2.78020 -0.19250 80.61560 -0.31050 10.19620 -2.77290 -0.35110 85.14040 -0.29960 10.22060 -2.76520 -0.31970 85.93110 -0.31630 10.33740 -2.75810 -0.09320 80.32490 -0.32810 10.35780 -2.75060 -0.20390 84.24660 -0.32140 10.37740 -2.74280 -0.16590 84.00810 -0.34130 10.49990 -2.73490 -0.16960 85.32660 -0.35180 10.56130 -2.72760 -0.10450 84.07700 -0.35540 10.57260 -2.71970 -0.14970 86.05340 -0.37180 10.71420 -2.71190 -0.17810 87.90710 -0.38450 10.83340 -2.70380 -0.33930 93.00450 -0.38740 10.91310 -2.69610 -0.38760 95.50160 -0.39660 11.00860 -2.68860 -0.32880 94.08060 -0.40830 11.09110 -2.68100 -0.30170 94.46150 -0.41200 11.11920 -2.67320 -0.41280 98.33380 -0.41020 11.16870 -2.66470 -0.36310 97.72390 -0.44430 11.37920 -2.65750 -0.23480 95.61400 -0.44560 11.36830 -2.65000 -0.18340 95.37850 -0.45310 11.42450 -2.64270 -0.00550 91.64620 -0.45750 11.39900 -2.63500 -0.01440 93.77810 -0.47230 11.51120 -2.62700 0.06420 92.05410 -0.49730 11.68750 -2.61890 0.09490 91.36600 -0.51050 11.78700 -2.61160 -0.07050 96.35640 -0.49470 11.80490 -2.60350 0.06190 94.19360 -0.51810 11.92690 -2.59530 0.10290 93.37860 -0.53650 12.06180 "
    s = s + "-2.58670 0.25700 90.31390 -0.57410 12.25580 -2.57890 0.26890 89.29890 -0.57420 12.32190 -2.57040 0.23530 90.85390 -0.58510 12.43910 -2.56150 0.20710 92.42280 -0.60990 12.61810 -2.55360 0.17960 94.50420 -0.60360 12.63790 -2.54530 0.18440 95.31280 -0.61250 12.71670 -2.53740 0.15370 96.51720 -0.61830 12.81730 -2.52800 0.12930 96.91870 -0.64890 13.06200 -2.51900 0.03510 100.07480 -0.66430 13.24150 -2.51130 -0.13160 105.04060 -0.64280 13.26100 -2.50240 -0.12100 106.02820 -0.66690 13.46310 -2.49300 0.12350 102.03390 -0.73740 13.84590 -2.48480 0.31660 97.24660 -0.74680 13.93710 -2.47690 0.29650 97.06450 -0.73750 14.02220 -2.46870 0.48180 93.89800 -0.74380 14.04600 -2.46020 0.76140 87.14430 -0.77580 14.25850 -2.45160 0.96790 81.78130 -0.80190 14.47000 -2.44150 0.98640 82.35220 -0.85190 14.86000 -2.43220 1.09840 80.45360 -0.88660 15.12640 -2.42270 1.21800 80.25200 -0.92270 15.36940 -2.41240 1.29020 80.11890 -0.96650 15.68440 -2.40340 1.00490 88.16490 -0.93980 15.73340 -2.39250 1.30340 82.42080 -1.01060 16.11560 -2.38340 1.16520 85.64300 -1.00720 16.32060 -2.37350 1.26390 83.67010 -1.03690 16.59560 -2.36310 1.25220 84.81910 -1.06210 16.84220 -2.35340 1.29490 82.70420 -1.07010 17.05580 -2.34290 1.56120 76.65630 -1.13100 17.51690 -2.33250 1.67620 74.49360 -1.16260 17.83070 -2.32320 1.31060 85.14070 -1.10210 17.74640 "
    s = s + "-2.31140 1.23790 89.52470 -1.16000 18.29140 -2.30060 0.99010 95.79920 -1.14290 18.47680 -2.29120 1.05760 94.39410 -1.11020 18.46610 -2.27930 1.07230 95.01330 -1.14790 18.89770 -2.26790 0.75310 105.40190 -1.14760 19.25790 -2.25620 1.01440 99.05890 -1.17240 19.55250 -2.24460 1.20800 95.23190 -1.18840 19.84250 -2.23270 1.07170 98.59770 -1.17930 20.13220 -2.21870 1.17380 96.30530 -1.21590 20.57890 -2.20780 1.59570 85.26480 -1.21100 20.74530 -2.19330 1.56140 89.17960 -1.23600 21.09610 -2.18030 1.98640 79.17400 -1.25500 21.37710 -2.16510 1.90400 82.51890 -1.29000 21.97810 -2.15010 2.08900 77.98200 -1.30380 22.30030 -2.13560 1.64340 93.34620 -1.26620 22.64190 -2.11960 1.42700 99.88320 -1.26120 23.18170 -2.10520 1.37720 102.77620 -1.22010 23.50170 -2.09030 1.42240 100.41340 -1.16050 23.84500 -2.07370 1.42180 102.28910 -1.11770 24.22030 -2.05530 0.86280 115.23300 -1.03850 24.67260 -2.03380 0.65200 121.09240 -1.03460 25.57110 -2.01230 0.98740 116.62150 -1.03930 26.23780 -1.99230 0.53740 128.56330 -0.85040 26.35070 -1.97100 0.40260 135.04640 -0.66590 26.43830 -1.94670 -0.13550 144.18440 -0.42760 26.65600 -1.91680 1.38830 110.43690 -0.41550 27.57880 -1.88900 1.88450 95.73820 -0.04550 27.28230 -1.85540 1.96250 92.76150 0.45910 27.01760 -1.79940 3.91540 41.73670 0.72170 28.48880 -1.79020 4.97320 10.34630 0.94140 27.69700 "
    s = s + "-1.77850 2.17390 77.80610 1.43090 26.94360 -1.76990 2.07720 79.47290 1.91100 25.29670 -1.74740 2.06610 92.57590 2.01180 26.09110 -1.72490 4.22370 28.80810 2.10090 26.51740 -1.70040 2.89140 61.79790 2.49830 26.76880 -1.68130 -1.45530 167.43250 3.79150 23.18720 -1.63780 1.49920 94.46150 4.08050 24.45040 -1.57520 -14.94810 586.91590 6.85490 18.95260 -1.52510 -7.84130 428.10710 7.34160 20.31480 -1.38440 3.37100 159.37320 5.58250 33.62430 -1.35520 -35.26250 1213.10730 11.49560 15.77780 "
    return(rowshape(strtoreal(tokens(s)), 221))
}

real matrix _xtfp_rs2()
{
    string scalar s
    s = "-6.29620 -78.79380 19.75070 2.01460 4.70340 -6.28350 -82.96850 207.22040 7.57910 -24.01650 -6.07060 -41.54670 -767.61110 4.28140 -14.58590 -5.89570 -19.17430 -1071.19930 1.45020 -2.94760 -5.73100 -11.52580 -1044.10650 1.01100 -0.84590 -5.64490 -17.87860 -713.09670 2.09150 -4.85950 -5.57100 -14.15060 -765.56720 1.84890 -3.37160 -5.51920 -14.19870 -702.98930 2.07870 -4.03050 -5.49020 -18.28630 -538.43550 3.02220 -7.80440 -5.43510 -13.39820 -618.67090 2.23650 -4.80040 -5.38340 -10.31280 -645.45520 1.54420 -1.78480 -5.34490 -9.15600 -643.82080 1.36310 -1.22380 -5.31430 -8.00850 -629.88450 1.15440 -0.48040 -5.18060 -9.75090 -477.83940 0.89650 1.98950 -5.10080 -11.34190 -360.60740 1.22040 1.22650 -5.02400 -11.87670 -288.34780 1.17740 1.49790 -4.96420 -10.18080 -291.35460 1.02470 2.03570 -4.90510 -11.04770 -222.34470 0.76610 4.02990 -4.86150 -9.46230 -234.41120 0.63470 4.66100 -4.81770 -10.03000 -186.66250 0.60790 5.02180 -4.78100 -9.24990 -181.63590 0.52200 5.39900 -4.74500 -8.86660 -176.25860 0.40880 6.22410 -4.71520 -8.93900 -150.84140 0.41200 6.42320 -4.68680 -8.02030 -160.07860 0.35820 6.58390 -4.65870 -7.77530 -146.25590 0.28970 6.92160 -4.63200 -8.43080 -110.79610 0.31780 6.98200 -4.60800 -7.60170 -120.09930 0.23790 7.30260 -4.58250 -7.26070 -114.09170 0.15110 7.70560 -4.56120 -7.36990 -100.37290 0.21600 7.45990 -4.53940 -7.67020 -77.26290 0.20150 7.74450 "
    s = s + "-4.51780 -7.19690 -78.23910 0.14140 7.99310 -4.49650 -6.64650 -81.30550 0.00210 8.76700 -4.47570 -6.04100 -87.84060 -0.10880 9.24180 -4.45640 -7.01370 -49.63240 -0.07420 9.40460 -4.43770 -7.23870 -34.47990 -0.08990 9.59970 -4.41990 -6.85480 -36.26540 -0.16460 9.99310 -4.40380 -7.12630 -23.30110 -0.10910 9.79720 -4.38660 -6.96760 -18.20350 -0.15850 10.13360 -4.36890 -6.73800 -17.59280 -0.22880 10.52010 -4.35350 -6.87470 -5.98520 -0.23580 10.63560 -4.33890 -6.54030 -7.53110 -0.25550 10.69270 -4.32430 -5.98540 -12.88210 -0.31400 10.89410 -4.31010 -5.48930 -19.60500 -0.36290 11.07630 -4.29520 -5.53840 -11.57050 -0.40740 11.40390 -4.28170 -5.65070 -1.25370 -0.39700 11.42040 -4.26820 -5.57630 1.51320 -0.40430 11.47540 -4.25530 -5.22850 -1.10360 -0.43020 11.51410 -4.24310 -5.25870 6.58900 -0.42000 11.49560 -4.22940 -5.49140 18.55960 -0.43650 11.69170 -4.21660 -5.48440 23.42660 -0.44980 11.80340 -4.20540 -5.50540 28.99680 -0.42370 11.71770 -4.19330 -5.31860 29.48610 -0.46340 11.93310 -4.18140 -5.08830 27.69570 -0.50090 12.13880 -4.17020 -5.02990 29.78290 -0.50810 12.23160 -4.15880 -5.09340 37.61620 -0.50430 12.24440 -4.14710 -5.05790 41.37220 -0.52330 12.35220 -4.13590 -4.98140 44.51960 -0.55300 12.51890 -4.12420 -4.99530 48.51820 -0.57860 12.68880 -4.11380 -5.20930 58.79920 -0.55630 12.66800 -4.10280 -5.44800 68.68840 -0.53770 12.61670 "
    s = s + "-4.09170 -5.34750 71.10950 -0.58050 12.87140 -4.08260 -5.28150 72.33720 -0.54780 12.70870 -4.07130 -5.29380 77.23130 -0.57110 12.86630 -4.05980 -4.98540 74.31580 -0.64360 13.17080 -4.04940 -4.66940 69.48220 -0.69260 13.38750 -4.03990 -4.55870 70.31840 -0.69510 13.40840 -4.03050 -4.62250 74.85820 -0.67710 13.36340 -4.02130 -4.60260 78.18440 -0.67730 13.40750 -4.01120 -4.30980 73.51830 -0.72110 13.61540 -4.00180 -4.25290 74.03910 -0.72090 13.61620 -3.99260 -4.16640 75.90570 -0.74620 13.76920 -3.98340 -4.45270 86.90210 -0.72830 13.80150 -3.97440 -4.57120 94.33610 -0.71630 13.80140 -3.96420 -4.60620 98.24740 -0.74950 14.03110 -3.95550 -4.61130 101.34610 -0.74710 14.07580 -3.94560 -4.85500 111.98340 -0.74810 14.16930 -3.93670 -4.67240 109.74880 -0.76090 14.19670 -3.92800 -4.57340 109.19150 -0.77500 14.29730 -3.91930 -4.45240 109.81350 -0.79640 14.39800 -3.90980 -4.44500 112.13800 -0.81740 14.51960 -3.90050 -4.52930 118.89790 -0.83620 14.67770 -3.89210 -4.63080 124.93460 -0.81850 14.64810 -3.88280 -4.54770 125.80500 -0.85170 14.85080 -3.87490 -4.54020 129.55840 -0.84960 14.89370 -3.86620 -4.25660 124.19080 -0.87950 14.99720 -3.85830 -4.01260 120.16030 -0.89360 15.04510 -3.85060 -4.09690 124.60420 -0.86990 14.98850 -3.84240 -4.31360 133.22500 -0.86490 15.10360 -3.83340 -4.25490 133.70940 -0.90520 15.37190 -3.82560 -4.25880 136.81820 -0.91420 15.48530 "
    s = s + "-3.81740 -4.18470 136.93520 -0.93440 15.60510 -3.80820 -3.90780 131.89120 -1.00590 15.99470 -3.80050 -3.97900 136.61210 -0.99350 15.97280 -3.79170 -4.01010 139.66650 -1.00650 16.05600 -3.78440 -4.17340 145.05070 -0.98370 16.05060 -3.77600 -4.18170 147.73880 -0.99680 16.16140 -3.76820 -4.09920 147.45940 -1.01380 16.28040 -3.76050 -4.32220 156.61860 -1.00200 16.33680 -3.75140 -4.21070 157.18090 -1.05560 16.63310 -3.74320 -4.15650 158.24840 -1.06780 16.73110 -3.73550 -4.19370 160.52460 -1.07460 16.83580 -3.72860 -3.97280 157.24470 -1.07590 16.79680 -3.72100 -3.72590 152.65110 -1.10830 16.97730 -3.71340 -3.68050 153.23640 -1.11870 17.08970 -3.70570 -3.64410 153.59690 -1.13170 17.19270 -3.69760 -3.67670 156.28340 -1.14250 17.30160 -3.68950 -3.54600 154.80860 -1.17410 17.50620 -3.68150 -3.70860 162.26400 -1.16780 17.52330 -3.67290 -3.68410 163.48210 -1.20210 17.77440 -3.66520 -3.70420 165.61270 -1.20070 17.81030 -3.65580 -3.57120 164.24500 -1.25050 18.04840 -3.64880 -3.54710 165.27670 -1.23830 17.99360 -3.64020 -3.56050 167.83020 -1.26250 18.16080 -3.63270 -3.38960 166.20240 -1.28850 18.30580 -3.62510 -3.50700 171.19850 -1.29200 18.42390 -3.61690 -3.47440 173.43030 -1.31740 18.60950 -3.60960 -3.39730 173.33530 -1.31550 18.62960 -3.60140 -3.33630 173.45870 -1.33850 18.77410 -3.59290 -3.25460 172.89080 -1.36420 18.89420 -3.58510 -3.20460 173.25950 -1.37740 19.01430 "
    s = s + "-3.57710 -3.22180 175.22020 -1.38660 19.09740 -3.56880 -3.12360 174.33530 -1.41660 19.27210 -3.56060 -3.08320 175.82820 -1.42910 19.32770 -3.55380 -3.05110 176.67420 -1.40980 19.23960 -3.54660 -3.21600 183.38620 -1.38600 19.18760 -3.53840 -3.09260 183.03180 -1.43180 19.47950 -3.53040 -3.25900 188.94560 -1.43180 19.59900 -3.52290 -3.19280 187.42640 -1.44500 19.73490 -3.51490 -2.89390 181.46830 -1.47880 19.87050 -3.50760 -2.77110 179.50840 -1.49350 19.99520 -3.50000 -2.64900 178.57740 -1.50420 20.04830 -3.49160 -2.44900 175.74530 -1.55540 20.33540 -3.48280 -2.59080 181.35480 -1.57370 20.53480 -3.47530 -2.66360 185.57160 -1.56560 20.55840 -3.46760 -2.60180 185.82820 -1.57250 20.63260 -3.45950 -2.47750 184.66440 -1.60480 20.82410 -3.45190 -2.33770 182.76960 -1.61950 20.90230 -3.44400 -2.29410 184.38960 -1.63640 21.02130 -3.43610 -2.13980 181.67540 -1.66170 21.18070 -3.42780 -2.03780 181.01260 -1.68100 21.27550 -3.42030 -2.18720 187.07030 -1.65720 21.21500 -3.41150 -2.16350 188.92950 -1.68750 21.41170 -3.40360 -2.27360 193.51050 -1.68120 21.47540 -3.39490 -2.00680 188.05270 -1.72950 21.70540 -3.38690 -1.77000 184.41220 -1.76310 21.89460 -3.37910 -1.88620 188.88620 -1.75700 21.98060 -3.37090 -1.80590 188.40180 -1.78090 22.17270 -3.36240 -1.80950 190.81370 -1.79450 22.32160 -3.35370 -1.83320 194.69150 -1.81800 22.51200 -3.34520 -1.83300 195.52890 -1.83520 22.69480 "
    s = s + "-3.33830 -1.93460 199.37070 -1.79340 22.56850 -3.33080 -1.92900 200.66990 -1.78840 22.62880 -3.32170 -1.87350 200.92620 -1.82740 22.91090 -3.31270 -1.48430 192.18320 -1.88340 23.18200 -3.30380 -1.13740 185.04880 -1.94090 23.46510 -3.29420 -1.32320 191.48320 -1.96470 23.73220 -3.28550 -1.06170 185.69700 -2.01000 24.02590 -3.27570 -0.95310 183.65610 -2.07010 24.43700 -3.26670 -1.02160 186.61530 -2.07490 24.58330 -3.25940 -1.05130 188.68970 -2.03400 24.45640 -3.24980 -0.92560 188.30330 -2.07440 24.69220 -3.24130 -0.93990 188.15760 -2.06690 24.77290 -3.23250 -0.93450 190.56570 -2.07520 24.91830 -3.22260 -0.65170 185.76660 -2.14840 25.35570 -3.21330 -0.63760 188.25280 -2.16860 25.53450 -3.20290 -0.97960 198.41220 -2.18500 25.87430 -3.19230 -0.93540 198.41040 -2.22910 26.21080 -3.18230 -1.11170 203.69300 -2.22630 26.34970 -3.17270 -1.28600 209.12110 -2.23390 26.60700 -3.16340 -1.19520 207.72040 -2.24580 26.78260 -3.15430 -1.47930 215.81430 -2.22010 26.91110 -3.14390 -1.20900 210.81070 -2.26370 27.20260 -3.13390 -1.28220 214.17670 -2.26680 27.39970 -3.12440 -1.55520 223.42230 -2.23330 27.40710 -3.11370 -1.39150 222.22700 -2.27960 27.75430 -3.10370 -1.32110 221.08150 -2.26950 27.83100 -3.09220 -1.27990 222.05250 -2.29700 28.09100 -3.08110 -1.24950 222.91360 -2.33140 28.45000 -3.07060 -0.88310 215.69840 -2.36640 28.72840 -3.06050 -0.80490 215.43350 -2.36990 28.94080 "
    s = s + "-3.04910 -0.77840 216.27940 -2.39340 29.28790 -3.03900 -0.83490 218.70930 -2.36230 29.37090 -3.02630 -0.52610 213.61650 -2.43030 29.85410 -3.01330 -0.59430 218.25530 -2.46600 30.24490 -3.00130 -0.49920 218.49260 -2.47310 30.43040 -2.98870 -0.48190 219.01390 -2.49350 30.82620 -2.97610 -0.33260 219.15480 -2.51800 31.14170 -2.96320 -0.16330 216.26150 -2.55170 31.58520 -2.94930 0.19100 209.76740 -2.62150 32.20400 -2.93500 -0.39950 225.58480 -2.57910 32.49720 -2.92110 -0.13910 220.47620 -2.58930 32.82650 -2.90790 -0.08970 223.46710 -2.55960 32.96170 -2.89370 -0.36980 233.37630 -2.50640 33.06540 -2.87770 -0.38800 238.08030 -2.51840 33.54730 -2.86220 -0.33390 236.52960 -2.48790 33.84550 -2.84520 -0.14860 232.85980 -2.51690 34.51140 -2.82850 0.04930 227.42580 -2.53270 35.19230 -2.81080 0.34150 222.43350 -2.54670 35.82680 -2.79350 1.18870 204.19650 -2.54090 36.12720 -2.77720 1.23360 205.16460 -2.39880 35.98950 -2.75960 1.67890 194.71290 -2.24390 35.71770 -2.73580 2.74400 170.79140 -2.29940 36.57470 -2.71000 2.75850 177.04350 -2.23380 37.07350 -2.68420 3.30050 165.08970 -2.21720 38.07910 -2.65620 3.17590 172.09200 -2.06570 38.66450 -2.62520 4.07420 159.46430 -1.87220 38.63280 -2.58590 3.67480 169.88590 -1.59270 38.98760 -2.54990 4.62940 145.18890 -0.96510 37.74180 -2.49740 5.70000 113.95400 -0.34900 37.64810 -2.48910 4.80900 131.47750 0.06850 36.69760 "
    s = s + "-2.47480 4.42400 139.71380 0.32270 36.36350 -2.45440 7.51040 60.35070 0.24300 36.90390 -2.43910 7.48500 67.28880 0.62890 36.06160 -2.40910 7.93160 51.60130 0.52500 38.57960 -2.38360 9.92440 -10.60100 0.69380 39.30730 -2.35820 9.50700 -8.32220 1.42500 37.51720 -2.33800 13.53700 -120.51330 2.37570 35.09100 -2.26040 10.74870 -123.90000 4.01290 33.58690 -2.18290 11.49650 -149.46190 5.06830 31.07880 -2.18360 12.24820 -146.22780 8.68580 19.43670 -2.15370 15.60450 -261.20410 10.86220 7.24730 "
    return(rowshape(strtoreal(tokens(s)), 221))
}

real matrix _xtfp_rs3()
{
    string scalar s
    s = "-6.91000 67.66360 -5216.85680 -10.63640 26.57220 -6.90380 52.93670 -4503.27670 -5.27510 6.42190 -6.60520 21.30750 -3300.94600 -5.82910 25.93260 -6.50490 9.11750 -2726.22670 -2.87670 15.10360 -6.35070 13.36610 -2431.34740 -2.82050 13.34150 -6.23470 8.91940 -2072.21010 -2.68740 13.22400 -6.15270 8.20010 -1944.02230 -2.75520 14.36500 -6.08390 2.69170 -1679.64720 -2.65240 16.04960 -6.04540 5.59070 -1653.90910 -2.46820 14.11530 -5.99680 7.97070 -1646.92610 -2.75410 15.11380 -5.96990 2.03070 -1437.97560 -2.01760 13.18840 -5.94370 0.66740 -1357.44440 -1.75830 12.46570 -5.91660 -1.63470 -1246.94090 -1.43770 11.28520 -5.79610 -2.47010 -1080.53090 -1.39610 11.93450 -5.70830 -4.55210 -915.33220 -1.26990 12.23430 -5.63810 -5.01110 -812.75140 -1.24390 12.54860 -5.57940 -4.65490 -758.00130 -1.32070 13.19650 -5.53190 -5.14530 -692.75010 -1.23350 13.19430 -5.47990 -7.31530 -579.15620 -1.21240 13.77800 -5.43720 -7.88430 -523.66410 -1.28170 14.77890 -5.40240 -8.00090 -490.65930 -1.23850 14.80700 -5.36880 -7.99360 -456.24040 -1.23680 14.88600 -5.33460 -8.16450 -416.38020 -1.34870 15.82150 -5.30440 -9.02630 -361.42380 -1.39060 16.64060 -5.27770 -9.96770 -307.63960 -1.31030 16.63600 -5.25180 -11.41090 -245.89000 -1.17430 16.43590 -5.22900 -11.37860 -223.04630 -1.13990 16.32760 -5.20470 -11.63860 -196.25510 -1.17590 16.87350 -5.18220 -11.42770 -181.09860 -1.23000 17.26950 -5.16120 -10.88570 -178.45510 -1.30940 17.74670 "
    s = s + "-5.14190 -10.86010 -161.95470 -1.30280 17.88470 -5.12360 -11.13900 -139.58580 -1.25430 17.82540 -5.10560 -11.67170 -110.46620 -1.19380 17.75500 -5.08670 -11.59220 -98.10680 -1.20940 17.86720 -5.06970 -12.12810 -69.16290 -1.15450 17.80700 -5.05350 -12.45370 -46.37160 -1.11210 17.73840 -5.03670 -11.94510 -50.60250 -1.20300 18.28540 -5.01950 -11.39810 -52.21510 -1.30330 18.78500 -5.00260 -11.28670 -44.09240 -1.38120 19.29590 -4.98730 -10.83560 -45.99730 -1.44080 19.61970 -4.97150 -10.87820 -34.57500 -1.43730 19.61970 -4.95810 -11.12820 -17.97510 -1.38770 19.50440 -4.94450 -11.69490 8.37730 -1.32870 19.40860 -4.92930 -11.69470 16.92010 -1.36620 19.72800 -4.91650 -11.87370 31.05400 -1.32410 19.60450 -4.90370 -11.40120 28.31430 -1.33730 19.54700 -4.89030 -11.33610 34.61280 -1.35290 19.68400 -4.87640 -11.15460 40.71940 -1.41690 20.03210 -4.86440 -11.52880 59.74900 -1.38810 20.07850 -4.85100 -11.34000 62.30320 -1.44180 20.38540 -4.83740 -10.95780 61.33060 -1.52470 20.80510 -4.82630 -10.98820 71.25580 -1.51560 20.82190 -4.81460 -10.83120 72.97990 -1.55790 21.17820 -4.80280 -10.85060 82.95530 -1.56060 21.19930 -4.79160 -10.61850 83.41950 -1.58280 21.35890 -4.77940 -10.57140 88.92970 -1.63260 21.72810 -4.76900 -10.63270 96.19700 -1.61460 21.74370 -4.75760 -10.27140 93.62260 -1.68400 22.11110 -4.74630 -10.39640 103.58980 -1.69690 22.30570 -4.73490 -10.41740 109.50760 -1.72630 22.58580 "
    s = s + "-4.72360 -10.37260 114.55230 -1.75560 22.82220 -4.71400 -10.37850 120.68530 -1.75680 22.94250 -4.70350 -10.09280 118.50900 -1.80800 23.23560 -4.69430 -9.99480 121.17350 -1.78520 23.15810 -4.68450 -10.42890 137.97220 -1.74340 23.13770 -4.67460 -10.49840 146.93100 -1.75020 23.24430 -4.66480 -10.68410 158.32490 -1.72830 23.18930 -4.65570 -10.64810 164.72150 -1.71690 23.13930 -4.64550 -10.45830 165.18540 -1.75370 23.38130 -4.63640 -10.31320 165.79200 -1.75700 23.42300 -4.62760 -10.45570 176.47130 -1.73940 23.42200 -4.61750 -10.15890 173.01080 -1.79720 23.71760 -4.60700 -10.10870 175.93020 -1.85070 24.11820 -4.59720 -9.75810 169.87640 -1.91240 24.47340 -4.58830 -9.75680 174.06290 -1.92390 24.60720 -4.57950 -9.44250 170.47580 -1.94790 24.67830 -4.57010 -9.25760 170.26550 -1.98410 24.87260 -4.56120 -9.31350 178.37260 -1.99160 24.98910 -4.55230 -9.23800 181.30820 -2.00760 25.11240 -4.54380 -9.31500 187.00770 -1.99370 25.12550 -4.53510 -9.18960 188.33000 -2.01510 25.26600 -4.52610 -8.93690 185.67730 -2.05090 25.44760 -4.51740 -8.65500 181.89970 -2.09060 25.67040 -4.50860 -8.52320 183.74890 -2.11760 25.80590 -4.50060 -8.49890 189.08670 -2.12290 25.88250 -4.49210 -8.21820 185.83620 -2.15850 26.07430 -4.48300 -8.19370 190.97850 -2.18700 26.24280 -4.47500 -8.08860 193.50720 -2.19050 26.27490 -4.46750 -8.07220 197.30130 -2.17820 26.23400 -4.45950 -8.17920 204.86940 -2.16450 26.24030 "
    s = s + "-4.45030 -7.83420 200.06220 -2.24230 26.65490 -4.44230 -7.59650 197.54140 -2.26320 26.75520 -4.43400 -7.50620 199.18000 -2.27540 26.84520 -4.42530 -7.24990 196.01300 -2.31670 27.05070 -4.41800 -7.36480 202.04140 -2.30100 27.09910 -4.40950 -7.43460 207.76150 -2.31440 27.24970 -4.40150 -7.37800 209.71110 -2.31960 27.31410 -4.39340 -7.14880 207.36510 -2.34530 27.43560 -4.38520 -7.16200 211.64430 -2.35790 27.55700 -4.37650 -6.93570 208.76830 -2.40310 27.78360 -4.36800 -6.89940 211.61400 -2.43320 27.99470 -4.35970 -6.90440 215.61080 -2.43980 28.06170 -4.35180 -6.93170 219.84490 -2.45080 28.17800 -4.34410 -6.83910 220.98740 -2.46980 28.31720 -4.33590 -6.60120 219.21440 -2.50950 28.51210 -4.32830 -6.42500 218.74900 -2.52640 28.61520 -4.32000 -6.37710 220.58310 -2.53850 28.69200 -4.31070 -5.96150 212.94260 -2.62240 29.07510 -4.30300 -5.93720 216.47490 -2.62390 29.11910 -4.29450 -5.49520 208.42310 -2.68930 29.41930 -4.28620 -5.53610 212.18010 -2.69660 29.53760 -4.27800 -5.35990 210.62760 -2.73010 29.73840 -4.27020 -5.48710 216.00260 -2.72540 29.82030 -4.26210 -5.64300 223.41310 -2.72440 29.88990 -4.25400 -5.39080 219.84230 -2.76840 30.15570 -4.24600 -5.23210 219.41220 -2.79720 30.29910 -4.23830 -5.45470 228.53640 -2.77840 30.29990 -4.23030 -5.33820 228.20950 -2.80360 30.46410 -4.22250 -5.14140 225.75030 -2.82990 30.61390 -4.21490 -4.99130 224.21630 -2.85480 30.79990 "
    s = s + "-4.20760 -5.18390 233.10450 -2.82730 30.75120 -4.20050 -5.18130 235.22150 -2.81790 30.75180 -4.19260 -5.08610 236.05760 -2.83320 30.86770 -4.18450 -5.20520 242.90520 -2.83160 30.93880 -4.17620 -5.16120 245.28940 -2.85250 31.07800 -4.16820 -4.85650 239.40920 -2.89160 31.25950 -4.16060 -4.99860 246.89620 -2.88770 31.34280 -4.15320 -4.86500 247.37310 -2.90840 31.50410 -4.14540 -4.69080 246.02740 -2.93640 31.66980 -4.13820 -4.67880 249.51580 -2.93160 31.72020 -4.13060 -4.52080 248.38930 -2.94070 31.77880 -4.12310 -4.52780 251.79880 -2.93580 31.78950 -4.11480 -4.56530 257.52560 -2.94940 31.92670 -4.10710 -4.36070 253.38740 -2.96990 32.06530 -4.09880 -4.25910 253.98940 -2.99140 32.20140 -4.09160 -4.33230 259.42210 -2.98040 32.26190 -4.08380 -4.44950 265.09520 -2.97550 32.36140 -4.07530 -4.35270 265.86040 -3.01920 32.67190 -4.06710 -4.17250 265.01510 -3.04510 32.84420 -4.05890 -4.04170 265.03010 -3.07710 33.05840 -4.05170 -3.91830 263.36110 -3.07310 33.12590 -4.04430 -3.63990 259.02240 -3.08950 33.20630 -4.03560 -3.51270 257.73420 -3.12770 33.48740 -4.02700 -3.50450 261.42000 -3.14060 33.58180 -4.01900 -3.53770 264.68130 -3.14840 33.74100 -4.01120 -3.50310 266.96700 -3.14400 33.77300 -4.00310 -3.45270 269.21730 -3.15230 33.87210 -3.99480 -3.16750 264.87650 -3.19290 34.09740 -3.98610 -3.33330 272.20860 -3.20980 34.38290 -3.97790 -3.42120 277.58350 -3.21120 34.50140 "
    s = s + "-3.96970 -3.26500 277.79090 -3.23960 34.71780 -3.96160 -3.26500 283.07280 -3.23990 34.77110 -3.95310 -2.86210 275.04350 -3.27900 34.94420 -3.94470 -2.95090 279.92570 -3.27990 35.09650 -3.93720 -2.89480 280.51860 -3.25920 35.07060 -3.92710 -2.91790 285.18770 -3.30720 35.43560 -3.91870 -2.67780 282.55230 -3.33180 35.59140 -3.90940 -2.89530 292.70920 -3.34100 35.80830 -3.90050 -2.63010 288.94020 -3.36270 35.92800 -3.89160 -2.67670 293.67800 -3.36130 36.03600 -3.88290 -2.74070 298.04870 -3.35670 36.12970 -3.87490 -2.68470 299.38410 -3.32970 36.06600 -3.86620 -2.46040 296.80180 -3.35600 36.30390 -3.85640 -2.70210 306.17040 -3.35630 36.49910 -3.84760 -2.74220 308.58290 -3.33810 36.52810 -3.83800 -2.27420 299.75510 -3.39080 36.79080 -3.82960 -2.09300 296.98300 -3.37750 36.80370 -3.82010 -2.32520 306.19830 -3.36780 37.00600 -3.81090 -2.16470 304.63750 -3.36840 37.09620 -3.80080 -2.54570 318.28220 -3.33920 37.17380 -3.79110 -2.48140 321.45630 -3.34560 37.31890 -3.77970 -2.04260 313.01580 -3.43130 37.80740 -3.76990 -2.09320 317.67610 -3.43200 37.96350 -3.75890 -1.94870 316.11520 -3.46160 38.24320 -3.74920 -2.15480 323.61720 -3.43320 38.36160 -3.73890 -2.00550 322.57870 -3.45380 38.61110 -3.72740 -1.89720 323.57690 -3.47850 38.83460 -3.71570 -1.77000 324.31780 -3.51780 39.20600 -3.70370 -1.80220 329.88510 -3.54510 39.56330 -3.69200 -1.61720 328.49500 -3.58200 39.92310 "
    s = s + "-3.68030 -1.00480 316.36210 -3.64950 40.28720 -3.67050 -1.00690 319.74310 -3.60640 40.29230 -3.65740 -0.69320 315.36640 -3.67550 40.81920 -3.64390 -0.57980 317.40560 -3.73160 41.28280 -3.62990 -0.55290 320.25930 -3.78240 41.82510 -3.61840 -0.30360 315.58880 -3.78250 42.05120 -3.60450 0.16700 306.49940 -3.84250 42.51120 -3.59040 0.88270 292.60390 -3.94230 43.15670 -3.57720 0.93030 296.54880 -3.94770 43.48370 -3.56320 0.74930 304.91910 -3.93230 43.82350 -3.54810 1.21370 294.53450 -3.99300 44.41520 -3.53410 1.49620 288.44910 -4.00270 44.84990 -3.52000 1.53280 292.48990 -3.98290 45.14220 -3.50480 1.57990 292.93500 -3.99060 45.76710 -3.48710 2.16430 283.04950 -4.06030 46.41590 -3.47210 2.50990 279.32480 -4.01280 46.54670 -3.45380 1.51650 308.77890 -3.92100 46.93750 -3.43630 1.35720 318.69470 -3.83620 47.06920 -3.41850 1.52520 320.60060 -3.75230 47.21830 -3.39840 2.36200 304.23100 -3.80400 48.05200 -3.37670 3.22360 287.37220 -3.83590 48.86030 -3.35060 2.95720 299.97860 -3.82470 49.81570 -3.32580 2.66870 317.55390 -3.71370 50.20950 -3.30080 2.76230 316.61110 -3.58240 50.79790 -3.26950 3.21850 312.77700 -3.53220 51.83090 -3.23290 3.71360 302.33930 -3.49720 53.34960 -3.19600 4.46310 284.31080 -3.23600 53.83550 -3.15490 3.11900 310.97160 -2.45450 52.48210 -3.09990 3.04380 321.39410 -1.66410 51.85210 -3.08220 3.70140 307.93030 -1.68760 52.67400 "
    s = s + "-3.06990 5.54500 257.94410 -1.56070 52.65220 -3.04740 4.97720 283.42930 -1.44780 53.35640 -3.03280 4.86810 296.20800 -1.14680 53.12550 -3.01100 7.53810 242.40590 -1.09000 53.42330 -2.97950 8.27050 230.26810 -0.85740 54.36120 -2.94850 9.25850 201.00030 -0.22920 52.99540 -2.91460 6.42550 281.32690 1.25130 49.31620 -2.86060 9.59020 190.73670 3.00880 44.35400 -2.78540 -7.75940 702.26470 6.17950 34.91070 -2.79120 -39.49820 1640.39520 14.26230 2.43120 -2.74040 -65.61350 2347.03500 18.18710 -4.50650 "
    return(rowshape(strtoreal(tokens(s)), 221))
}

end