-- Copyright (c) 1991-2002, The Numerical ALgorithms Group Ltd.
-- All rights reserved.
-- Copyright (C) 2007, Gabriel Dos Reis.
-- All rights reserved.
--
-- Redistribution and use in source and binary forms, with or without
-- modification, are permitted provided that the following conditions are
-- met:
--
--     - Redistributions of source code must retain the above copyright
--       notice, this list of conditions and the following disclaimer.
--
--     - Redistributions in binary form must reproduce the above copyright
--       notice, this list of conditions and the following disclaimer in
--       the documentation and/or other materials provided with the
--       distribution.
--
--     - Neither the name of The Numerical ALgorithms Group Ltd. nor the
--       names of its contributors may be used to endorse or promote products
--       derived from this software without specific prior written permission.
--
-- THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS
-- IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
-- TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A
-- PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER
-- OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
-- EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
-- PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
-- PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
-- LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
-- NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
-- SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.


import '"macros"
)package "BOOT"

e02adf() ==
  htInitPage('"E02ADF - Least-squares curve fit, by polynomials, arbitrary data points", nil)
  htMakePage '(
    (domainConditions 
       (isDomain PI (PositiveInteger)))
    (text . "\windowlink{Manual Page}{manpageXXe02adf} for this routine ")
    (text . "\newline ")
    (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|e02adf| '|NagFittingPackage|)} for this routine")
    (text . "\newline \horizontalline ")
    (text . "Determines weighted least-squares polynomial approximations of ")
    (text . "degrees 0,1,...,k to the set of points {\it (} ")
    (text . "\htbitmap{xr}, \htbitmap{yr}{\it )} ")
    (text . "with weights \htbitmap{wr}, for r = 1,2,...,m. ")
    (text . "The polynomials are in the Chebyshev series form, the ")
    (text . "approximation of degree {\it i} being represented as ")
    (text . "\newline \center{\htbitmap{e02adf}} , where ")
    (text . "\htbitmap{xbar} is the normalised argument, which is ")
    (text . "related to the original variable {\it x} by the transformation ")
    (text . "\blankline  \center{\htbitmap{e02adf1}} ")
    (text . ",\htbitmap{xmin} and \htbitmap{xmax} being ")
    (text . "the values of \htbitmap{xr} respectively ")
    (text . "\blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} \newline ")
    (text . "Number of data points {\it m}:")
    (text . "\newline \tab{2} ")
    (bcStrings (6 11 m PI))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2}")
    (text . "Maximum degree required {\it k}:")
    (text . "\newline \tab{2} ")
    (bcStrings (6 3 k PI))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2} \newline ")
    (text . "First dimension of A, {\it nrows} \htbitmap{great=} {\it k+1}: ")
    (text . "\newline \tab{2} ")
    (bcStrings (6 50 nrows I))
    (text . "\blankline")
    (text . "\newline ")
    (text . "\menuitemstyle{} \tab{2} ")
    (text . "\newline \tab{2} ")
    (text . "Ifail value: ")
    (radioButtons ifail
        ("" " -1, Print error messages" minusOne)
        ("" "  1, Suppress error messages" one)))
  htMakeDoneButton('"Continue", 'e02adfSolve)
  htShowPage()

e02adfSolve htPage ==
  m :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'm)
    objValUnwrap htpLabelSpadValue(htPage, 'm)
  k :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'k)
    objValUnwrap htpLabelSpadValue(htPage, 'k)
  nrows :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'nrows)
    objValUnwrap htpLabelSpadValue(htPage, 'nrows)
  error := htpButtonValue(htPage,'ifail)
  ifail :=
    error = 'one => '1
    '-1
  (m = '11 and k ='3) => e02adfDefaultSolve(htPage,k,nrows,ifail)
  labelList :=
    "append"/[f(i) for i in 1..m] where f(i) ==
      prefix := ('"\newline \tab{2} ")
      middle := ('"\tab{22} ")
      post   := ('" \tab{42} ")
      xnam := INTERN STRCONC ('"x",STRINGIMAGE i)
      ynam := INTERN STRCONC ('"y",STRINGIMAGE i)
      znam := INTERN STRCONC ('"z",STRINGIMAGE i)
      num := INTERN STRCONC (STRINGIMAGE (QUOTIENT(i,10)),".",STRINGIMAGE (REM(i,10)))
      [['text,:prefix],['bcStrings,[10, num, xnam, 'F]], 
       ['text,:middle],['bcStrings,[10, 0.0, ynam, 'F]],
        ['text,:post],['bcStrings,[10, 0.0, znam, 'F]]]
  equationPart := [
     '(domainConditions 
        (isDomain P (Polynomial $EmptyMode))
         (isDomain S (String))
          (isDomain F (Float))
           (isDomain PI (PositiveInteger))),
            :labelList]
  page := htInitPage("E02ADF - Least-squares curve fit, by polynomials, arbitrary data points", htpPropertyList htPage)
  htSay '"\menuitemstyle{}\tab{2} Values of \space{1} "
  htSay '"\htbitmap{xr}: \tab{20} \menuitemstyle{}\tab{22} "
  htSay '"Values of \space{1} \htbitmap{fr}: \tab{40}"
  htSay '"\menuitemstyle{}\tab{42} Values of \space{1} "
  htSay '"\htbitmap{dr}: "
  htMakePage equationPart
  htSay '"\blankline "
  htMakeDoneButton('"Continue",'e02adfGen)
  htpSetProperty(page,'m,m)
  htpSetProperty(page,'k,k)
  htpSetProperty(page,'nrows,nrows)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()


e02adfDefaultSolve  (htPage,k,nrows,ifail) ==
  m := '11
  page := htInitPage('"E02ADF - Least-squares curve fit, by polynomials, arbitrary data points", htpPropertyList htPage)
  htMakePage '(
    (domainConditions 
       (isDomain F (Float)))
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2} Values of \space{1} ")
    (text . "\htbitmap{xr}: \tab{20} \menuitemstyle{}\tab{22} ")
    (text . "Values of \space{1} \htbitmap{yr}: \tab{40} ")
    (text . "\menuitemstyle{}\tab{42} Values of \space{1} ")
    (text . "\htbitmap{wr}: ")
    (text . "\newline \tab{2}")
    (bcStrings (10 "1.00" x1 F))
    (text . "\tab{22} ")
    (bcStrings (10 "10.40" y1 F))
    (text . "\tab{42} ")
    (bcStrings (10 "1.00" z1 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "2.10" x2 F))
    (text . "\tab{22} ")
    (bcStrings (10 "7.90" y2 F))
    (text . "\tab{42} ")
    (bcStrings (10 "1.00" z2 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "3.10" x3 F))
    (text . "\tab{22} ")
    (bcStrings (10 "4.70" y3 F))
    (text . "\tab{42} ")
    (bcStrings (10 "1.00" z3 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "3.90" x4 F))
    (text . "\tab{22} ")
    (bcStrings (10 "2.50" y4 F))
    (text . "\tab{42} ")
    (bcStrings (10 "1.00" z4 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "4.90" x5 F))
    (text . "\tab{22} ")
    (bcStrings (10 "1.20" y5 F))
    (text . "\tab{42} ")
    (bcStrings (10 "1.00" z5 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "5.80" x6 F))
    (text . "\tab{22} ")
    (bcStrings (10 "2.20" y6 F))
    (text . "\tab{42} ")
    (bcStrings (10 "0.80" z6 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "6.50" x7 F))
    (text . "\tab{22} ")
    (bcStrings (10 "5.10" y7 F))
    (text . "\tab{42} ")
    (bcStrings (10 "0.80" z7 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "7.10" x8 F))
    (text . "\tab{22} ")
    (bcStrings (10 "9.20" y8 F))
    (text . "\tab{42} ")
    (bcStrings (10 "0.70" z8 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "7.80" x9 F))
    (text . "\tab{22} ")
    (bcStrings (10 "16.10" y9 F))
    (text . "\tab{42} ")
    (bcStrings (10 "0.50" z9 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "8.40" x10 F))
    (text . "\tab{22} ")
    (bcStrings (10 "24.50" y10 F))
    (text . "\tab{42} ")
    (bcStrings (10 "0.30" z10 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "9.00" x11 F))
    (text . "\tab{22} ")
    (bcStrings (10 "35.30" y11 F))
    (text . "\tab{42} ")
    (bcStrings (10 "0.20" z11 F))
    (text . "\blankline"))
  htMakeDoneButton('"Continue",'e02adfGen)      
  htpSetProperty(page,'m,m)
  htpSetProperty(page,'k,k)
  htpSetProperty(page,'nrows,nrows)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()

e02adfGen htPage ==
  m := htpProperty(htPage,'m)
  k := htpProperty(htPage,'k)
  nrows := htpProperty(htPage,'nrows)
  ifail := htpProperty(htPage,'ifail)
  alist := htpInputAreaAlist htPage
  kplus1 := k + 1
  y := alist
  while y repeat
    right := STRCONC ((first y).1," ")
    y := rest y    
    mid :=  STRCONC ((first y).1," ")
    y := rest y
    left :=  STRCONC ((first y).1," ")
    y := rest y
    xlist := [left,:xlist]
    ylist := [mid,:ylist]
    wlist := [right,:wlist]
  xstring := bcwords2liststring xlist
  ystring := bcwords2liststring ylist
  wstring := bcwords2liststring wlist
  prefix := STRCONC('"e02adf(",STRINGIMAGE m,", ",STRINGIMAGE kplus1,", ")
  prefix := STRCONC(prefix,STRINGIMAGE nrows,", [",xstring,"],[",ystring,"],[")
  prefix := STRCONC(prefix,wstring,"],",STRINGIMAGE ifail,")")
  linkGen prefix

e02aef() ==
  htInitPage('"E02AEF - Evaluation of fitted polynomial in one variable from Chebyshev series form", nil)
  htMakePage '(
    (domainConditions 
      (isDomain F (Float))
       (isDomain PI (PositiveInteger)))
    (text . "\windowlink{Manual Page}{manpageXXe02aef} for this routine ")
    (text . "\newline ")
    (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|e02aef| '|NagFittingPackage|)} for this routine")
    (text . "\newline \horizontalline ")
    (text . "Evaluates a polynomial in Chabyshev series representation ")
    (text . "\newline \center{\htbitmap{e02aef}} ")
    (text . "\blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} \newline ")
    (text . "Number of terms in the series {\it n}:")
    (text . "\newline \tab{2} ")
    (bcStrings (6 4 n PI))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2}")
    (text . "\htbitmap{xbar}: ")
    (text . " \newline \tab{2} ")
    (bcStrings (6 "-1.0" xcap F))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{} \tab{2} ")
    (text . "\newline \tab{2} ")
    (text . "Ifail value: ")
    (radioButtons ifail
        ("" " -1, Print error messages" minusOne)
        ("" "  1, Suppress error messages" one)))
  htMakeDoneButton('"Continue", 'e02aefSolve)
  htShowPage()

e02aefSolve htPage ==
  n :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'n)
    objValUnwrap htpLabelSpadValue(htPage, 'n)
  xcap := htpLabelInputString(htPage,'xcap)
  error := htpButtonValue(htPage,'ifail)
  ifail :=
    error = 'one => '1
    '-1
  n = '4 => e02aefDefaultSolve(htPage,xcap,ifail)
  labelList :=
    "append"/[f(i) for i in 1..(n+1)] where f(i) ==
      prefix := ('"\newline \tab{15} ")
      anam := INTERN STRCONC ('"a",STRINGIMAGE i)
      [['text,:prefix],['bcStrings,[10, 0.0, anam, 'F]]]
  equationPart := [
     '(domainConditions 
        (isDomain P (Polynomial $EmptyMode))
         (isDomain S (String))
          (isDomain F (Float))
           (isDomain PI (PositiveInteger))),
            :labelList]
  page := htInitPage("E02AEF - Evaluation of fitted polynomial in one variable from Chebyshev series from", nil)
  htSay '"\menuitemstyle{}\tab{2} Enter the coefficients of {\it a(n+1)}:"
  htSay '"\blankline "
  htMakePage equationPart
  htSay '"\blankline "
  htMakeDoneButton('"Continue",'e02aefGen)
  htpSetProperty(page,'n,n)
  htpSetProperty(page,'xcap,xcap)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()


e02aefDefaultSolve  (htPage,xcap,ifail) ==
  n := '4
  page := htInitPage('"E02AEF - Evaluation of fitted polynomial in one variable from Chebyshev series form", nil)
  htMakePage '(
    (domainConditions 
       (isDomain F (Float)))
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2} ")
    (text . "Enter the coeffients of {\it a(n+1)}: ")
    (text . "\blankline ")
    (text . "\newline \tab{15} ")
    (bcStrings (10 "2.0000" a1 F))
    (text . "\newline \tab{15} ")
    (bcStrings (10 "0.5000" a2 F))
    (text . "\newline \tab{15} ")
    (bcStrings (10 "0.2500" a3 F))
    (text . "\newline \tab{15} ")
    (bcStrings (10 "0.1250" a4 F))
    (text . "\newline \tab{15} ")
    (bcStrings (10 "0.0625" a5 F))
    (text . "\blankline "))
  htMakeDoneButton('"Continue",'e02aefGen)      
  htpSetProperty(page,'n,n)
  htpSetProperty(page,'xcap,xcap)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()

e02aefGen htPage ==
  n := htpProperty(htPage,'n)
  xcap := htpProperty(htPage,'xcap)
  ifail := htpProperty(htPage,'ifail)
  alist := htpInputAreaAlist htPage
  nplus1 := n + 1
  y := alist
  while y repeat
    right := STRCONC ((first y).1," ")
    y := rest y    
    arrayList := [right,:arrayList]
  astring := bcwords2liststring arrayList
  prefix := STRCONC('"e02aef(",STRINGIMAGE nplus1,", [",astring ,"], ")
  prefix := STRCONC(prefix,STRINGIMAGE xcap,", ",STRINGIMAGE ifail,")")
  linkGen prefix

e02agf() ==
  htInitPage('"E02AGF - Least-squares polynomial fit, values and derivatives may be constrained, arbitrary data values",nil)
  htMakePage '(
    (domainConditions 
       (isDomain PI (PositiveInteger)))
    (text . "\windowlink{Manual Page}{manpageXXe02agf} for this routine ")
    (text . "\newline ")
    (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|e02agf| '|NagFittingPackage|)} for this routine")
    (text . "\newline \horizontalline ")
    (text . "Determines constrained least-squares polynomial approximations ")
    (text . "to the set of points {\it (\htbitmap{xr},\htbitmap{yr})} with ")
    (text . "weights \htbitmap{wr}, for r = 1,2,...,m. The values of the ")
    (text . "approximations and any number of their derivatives must be ")
    (text . "specified at a further set of points \htbitmap{xii}, ")
    (text . "for i = 1,2,...,{\it mf}. The total number of interpolating ")
    (text . "conditions is given by \center{\htbitmap{e02agf}} where ")
    (text . "\htbitmap{pi} is the highest order derivative ")
    (text . "specified at point \htbitmap{xii}.  The values ")
    (text . "\htbitmap{xr} and \htbitmap{xii} all lie ")
    (text . "in the interval [\htbitmap{xmin},")
    (text . "\htbitmap{xmax}]. The polynomials are given in ")
    (text . "Chebyshev series form, the approximation of degree {\it i} being")
    (text . " represented as\blankline  \center{\htbitmap{e02agf1}}")
    (text . "\newline, where \htbitmap{xbar} is the normalised ")
    (text . "argument, related to the original variable {\it x} by the ")
    (text . "transformation \newline \center{\htbitmap{e02adf1}} ")
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2}")
    (text . "\newline Number of data points {\it m}:")
    (text . "\tab{32} \menuitemstyle{}\tab{34}")
    (text . "Maximum degree required {\it k}:")
    (text . "\newline\tab{2} ")
    (bcStrings (6 5 m PI))
    (text . "\tab{34} ")
    (bcStrings (6 4 k PI))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2} \newline ")
    (text . "First dimension of A, {\it nrows \htbitmap{great=} k+1}: ")
    (text . "\newline \tab{2} ")
    (bcStrings (6 6 nrows I))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2}")
    (text . "\newline \htbitmap{xmin}: ")
    (text . "\tab{32} \menuitemstyle{}\tab{34}")
    (text . "\htbitmap{xmax}: ")
    (text . "\newline\tab{2} ")
    (bcStrings (6 "0.0" xmin F))
    (text . "\tab{34} ")
    (bcStrings (6 "4.0" xmax F))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2}")
    (text . "\newline Dimension of {\it xf} & {\it ip}, {\it mf}: ")
    (text . "\tab{32} \menuitemstyle{}\tab{34}")
    (text . "Dimension of {\it yf}, {\it lyf}:")
    (text . "\newline\tab{2} ")
    (bcStrings (6 2 mf PI))
    (text . "\tab{34} ")
    (bcStrings (6 15 lyf PI))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{} \tab{2} ")
    (text . "\newline \tab{2} ")
    (text . "Ifail value: ")
    (radioButtons ifail
        ("" " -1, Print error messages" minusOne)
        ("" "  1, Suppress error messages" one)))
  htMakeDoneButton('"Continue", 'e02agfSolve)
  htShowPage()

e02agfSolve htPage ==
  m :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'm)
    objValUnwrap htpLabelSpadValue(htPage, 'm)
  k :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'k)
    objValUnwrap htpLabelSpadValue(htPage, 'kplus1)
  nrows :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'nrows)
    objValUnwrap htpLabelSpadValue(htPage, 'nrows)
  xmin := htpLabelInputString(htPage,'xmin)
  xmax := htpLabelInputString(htPage,'xmax)
  mf :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'mf)
    objValUnwrap htpLabelSpadValue(htPage, 'mf)
  lyf :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'lyf)
    objValUnwrap htpLabelSpadValue(htPage, 'lyf)
  error := htpButtonValue(htPage,'ifail)
  ifail :=
    error = 'one => '1
    '-1
  (m = '5 and k ='4 and mf = '2 and lyf = '15) => e02agfDefaultSolve(htPage,nrows,xmin,xmax,ifail)
  labelList :=
    "append"/[f(i) for i in 1..m] where f(i) ==
      prefix := ('"\newline \tab{2} ")
      middle := ('"\tab{22} ")
      post   := ('" \tab{42} ")
      xnam := INTERN STRCONC ('"x",STRINGIMAGE i)
      ynam := INTERN STRCONC ('"y",STRINGIMAGE i)
      znam := INTERN STRCONC ('"z",STRINGIMAGE i)
      num := INTERN STRCONC (STRINGIMAGE (QUOTIENT(i,10)),".",STRINGIMAGE (REM(i,10)))
      [['text,:prefix],['bcStrings,[10, num, xnam, 'F]], 
       ['text,:middle],['bcStrings,[10, 0.0, ynam, 'F]],
        ['text,:post],['bcStrings,[10, 0.0, znam, 'F]]]
  xfList := 
    "append"/[g(j) for j in 1..mf] where g(j) ==
      xfnam := INTERN STRCONC ('"xf",STRINGIMAGE j)
      [['bcStrings,[6, 0.0, xfnam, 'F]]]
  prefix := ('"\blankline \newline \menuitemstyle{}\tab{2} Values of ")
  prefix := STRCONC(prefix,"{\it xf}: \newline \tab{2} ")
  xfList := [['text,:prefix],:xfList] 
  ipList := 
    "append"/[h(k) for k in 1..mf] where h(k) ==
      ipnam := INTERN STRCONC ('"ip",STRINGIMAGE k)
      [['bcStrings,[6, 0, ipnam, 'PI]]]
  prefix := ('"\blankline \newline \menuitemstyle{}\tab{2} Values of ")
  prefix := STRCONC(prefix,"{\it ip}: \newline \tab{2} ")
  ipList := [['text,:prefix],:ipList] 
  yfList :=
    "append"/[i(l) for l in 1..lyf] where i(l) ==
      prefix := ('"\newline \tab{2} ")
      yfnam := INTERN STRCONC ('"lyf",STRINGIMAGE l)
      [['text,:prefix],['bcStrings,[10, 0.0, yfnam, 'F]]]
  prefix := ('"\blankline \newline \menuitemstyle{}\tab{2} Values of ")
  prefix := STRCONC(prefix,"{\it yf}: \newline \tab{2} ")
  yfList := [['text,:prefix],:yfList] 
  equationPart := [
     '(domainConditions 
        (isDomain P (Polynomial $EmptyMode))
         (isDomain S (String))
          (isDomain F (Float))
           (isDomain PI (PositiveInteger))),
            :labelList,:xfList,:ipList,:yfList]
  page := htInitPage("E02AGF - Least-squares polynomial fit, values and derivatives may be constrained, arbitrary data values",nil)
  htSay '"\menuitemstyle{}\tab{2} Values of \space{1} "
  htSay '"\htbitmap{xr}: \tab{20} \menuitemstyle{}\tab{22} "
  htSay '"Values of \space{1} \htbitmap{yr}: \tab{40}"
  htSay '"\menuitemstyle{}\tab{42} Values of \space{1} "
  htSay '"\htbitmap{wr}: "
  htMakePage equationPart
  htSay '"\blankline "
  htMakeDoneButton('"Continue",'e02agfGen)
  htpSetProperty(page,'m,m)
  htpSetProperty(page,'k,k)
  htpSetProperty(page,'nrows,nrows)
  htpSetProperty(page,'nrows,nrows)
  htpSetProperty(page,'xmin,xmin)
  htpSetProperty(page,'xmax,xmax)
  htpSetProperty(page,'mf,mf)
  htpSetProperty(page,'lyf,lyf)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()


e02agfDefaultSolve  (htPage,nrows,xmin,xmax,ifail) == 
  m := '5
  k := '4
  mf := '2
  lyf := '15
  page := htInitPage('"E02AGF - Least-squares polynomial fit, values and derivativesby polynomials, arbitrary data points", htpPropertyList htPage)
  htMakePage '(
    (domainConditions 
       (isDomain F (Float)))
    (text . "\menuitemstyle{}\tab{2} Values of \space{1} ")
    (text . "\htbitmap{xr}: \tab{20} \menuitemstyle{}\tab{22} ")
    (text . "Values of \space{1} \htbitmap{yr}: \tab{40} ")
    (text . "\menuitemstyle{}\tab{42} Values of \space{1} ")
    (text . "\htbitmap{wr}: ")
    (text . "\newline \tab{2}")
    (bcStrings (10 "0.5" x1 F))
    (text . "\tab{22} ")
    (bcStrings (10 "0.03" y1 F))
    (text . "\tab{42} ")
    (bcStrings (10 "1.0" z1 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "1.0" x2 F))
    (text . "\tab{22} ")
    (bcStrings (10 "-0.75" y2 F))
    (text . "\tab{42} ")
    (bcStrings (10 "1.0" z2 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "2.0" x3 F))
    (text . "\tab{22} ")
    (bcStrings (10 "-1.0" y3 F))
    (text . "\tab{42} ")
    (bcStrings (10 "1.0" z3 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "2.5" x4 F))
    (text . "\tab{22} ")
    (bcStrings (10 "-0.1" y4 F))
    (text . "\tab{42} ")
    (bcStrings (10 "1.0" z4 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "3.0" x5 F))
    (text . "\tab{22} ")
    (bcStrings (10 "1.75" y5 F))
    (text . "\tab{42} ")
    (bcStrings (10 "1.0" z5 F))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2} Values of {\it xf}: ")
    (text . "\newline \tab{2} ")
    (bcStrings (6 "0.0" xf1 F))
    (bcStrings (6 "4.0" xf2 F))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2} Values of {\it ip}: ")
    (text . "\newline \tab{2} ")
    (bcStrings (6 1 ip1 PI))
    (bcStrings (6 0 ip2 PI))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2} Values of {\it yf}: ")
    (text . "\newline \tab{2} ")
    (bcStrings (10 "1.0" lyf1 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "-2.0" lyf2 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "9.0" lyf3 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "0.0" lyf4 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "0.0" lyf5 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "0.0" lyf6 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "0.0" lyf7 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "0.0" lyf8 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "0.0" lyf9 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "0.0" lyf10 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "0.0" lyf11 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "0.0" lyf12 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "0.0" lyf13 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "0.0" lyf14 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "0.0" lyf15 F)))
  htMakeDoneButton('"Continue",'e02agfGen)      
  htpSetProperty(page,'m,m)
  htpSetProperty(page,'k,k)
  htpSetProperty(page,'nrows,nrows)
  htpSetProperty(page,'xmin,xmin)
  htpSetProperty(page,'xmax,xmax)
  htpSetProperty(page,'mf,mf)
  htpSetProperty(page,'lyf,lyf)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()

e02agfGen htPage ==
  m := htpProperty(htPage,'m)
  k := htpProperty(htPage,'k)
  nrows := htpProperty(htPage,'nrows)
  xmin := htpProperty(htPage,'xmin)
  xmax := htpProperty(htPage,'xmax)
  mf := htpProperty(htPage,'mf)
  lyf := htpProperty(htPage,'lyf)
  ifail := htpProperty(htPage,'ifail)
  alist := htpInputAreaAlist htPage
  kplus1 := k + 1
  ipsum := 0
  y := alist
  for i in 1..lyf repeat
    yf := STRCONC((first y).1," ")
    yfList := [yf,:yfList]
    y := rest y
  yfstring := bcwords2liststring yfList
  for i in 1..mf repeat
    iptest := (first y).1
    iptestval := READ_-FROM_-STRING(iptest)
    ipsum := ipsum + iptestval
    ip := STRCONC(iptest," ")
    iptestList := [iptestval,:iptestList]
    ipList := [ip,:ipList]
    y := rest y
  ipstring := bcwords2liststring ipList
  ipmax := APPLY ('MAX, iptestList)
  n := mf + ipsum
  for i in 1..mf repeat
    xf := STRCONC((first y).1," ")
    xfList := [xf,:xfList]
    y := rest y
  xfstring := bcwords2liststring xfList
  while y repeat
    right := STRCONC ((first y).1," ")
    y := rest y    
    mid :=  STRCONC ((first y).1," ")
    y := rest y
    left :=  STRCONC ((first y).1," ")
    y := rest y
    xlist := [left,:xlist]
    ylist := [mid,:ylist]
    wlist := [right,:wlist]
  xstring := bcwords2liststring xlist
  ystring := bcwords2liststring ylist
  wstring := bcwords2liststring wlist
  wrktest1 := 4*m + 3*kplus1
  wrktest2 := 8*n + 5*ipmax + mf +10
  wrktestlist := [wrktest1,wrktest2]
  wrkmax := APPLY ('MAX, wrktestlist)
  lwrk := wrkmax + 2*n + 2
  liwrk := 2*mf + 2
  prefix := STRCONC('"e02agf(",STRINGIMAGE m,", ",STRINGIMAGE kplus1,", ")
  prefix := STRCONC(prefix,STRINGIMAGE nrows,", ",xmin,", ",xmax,", [",xstring)
  prefix := STRCONC(prefix,"],[",ystring,"],[",wstring,"],",STRINGIMAGE mf)
  prefix := STRCONC(prefix,", [",xfstring,"],[",yfstring,"],")
  prefix := STRCONC(prefix,STRINGIMAGE lyf,", [",ipstring,"]::Matrix Integer,")
  prefix := STRCONC(prefix,STRINGIMAGE lwrk,", ",STRINGIMAGE liwrk,", ")
  prefix := STRCONC(prefix,STRINGIMAGE ifail,")")
  linkGen prefix

e02ahf() ==
  htInitPage('"E02AHF - Derivative of fitted polynomial in Chebyshev series",nil)
  htMakePage '(
    (domainConditions
      (isDomain F (Float)) 
       (isDomain PI (PositiveInteger)))
    (text . "\windowlink{Manual Page}{manpageXXe02ahf} for this routine ")
    (text . "\newline ")
    (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|e02ahf| '|NagFittingPackage|)} for this routine")
    (text . "\newline \horizontalline ")
    (text . "Determines the indefinite integral of the Chebyshev series ")
    (text . "representation \newline \center{\htbitmap{e02ahf1}} ")
    (text . "of a polynomial, where \htbitmap{xbar} is the ")
    (text . "normalised argument, related to the original variable x by the ")
    (text . "transformation \blankline \center{\htbitmap{e02adf1}}")
    (text . ",\htbitmap{xmin} and \htbitmap{xmax} being ")
    (text . "minimum and maximum values of {\it x} respectively. The integral")
    (text . " polynomial has the form ")
    (text . "\blankline \center{\htbitmap{e02ahf}}")
    (text . "\blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} \newline ")
    (text . "Degree of the polynomial {\it n}:")
    (text . "\newline \tab{2} ")
    (bcStrings (6 6 n PI))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2}")
    (text . "\newline \htbitmap{xmin}: ")
    (text . "\tab{32} \menuitemstyle{}\tab{34}")
    (text . "\htbitmap{xmax}: ")
    (text . "\newline\tab{2} ")
    (bcStrings (6 "-0.5" xmin F))
    (text . "\tab{34} ")
    (bcStrings (6 "2.5" xmax F))
--    (text . "\blankline ")
--    (text . "\newline ")
--    (text . "\menuitemstyle{}\tab{2}")
--    (text . "\newline Dimension of array {\it a}, {\it la}: ")
--    (text . "\tab{32} \menuitemstyle{}\tab{34}")
--    (text . "Dimension of {\it adif}, {\it ladif}: ") 
--    (text . "\newline\tab{2} ")
--    (bcStrings (6 7 la PI))
--    (text . "\tab{34} ")
--    (bcStrings (6 7 ladif PI))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2}")
    (text . "\newline Increment of array {\it a}, {\it ia1}: ")
    (text . "\tab{32} \menuitemstyle{}\tab{34}")
    (text . "\newline Increment of array {\it adif}, {\it ladif1}: ")
    (text . "\newline\tab{2} ")
    (bcStrings (6 1 iaone PI))
    (text . "\tab{34} ")
    (bcStrings (6 1 ladifone PI))
    (text . "\blankline")
    (text . "\newline \menuitemstyle{} \tab{2} ")
    (text . "\newline \tab{2} ")
    (text . "Ifail value: ")
    (radioButtons ifail
        ("" " -1, Print error messages" minusOne)
        ("" "  1, Suppress error messages" one)))
  htMakeDoneButton('"Continue", 'e02ahfSolve)
  htShowPage()

e02ahfSolve htPage ==
  n :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'n)
    objValUnwrap htpLabelSpadValue(htPage, 'n)
  xmin := htpLabelInputString(htPage,'xmin)
  xmax := htpLabelInputString(htPage,'xmax)
  iaone :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'iaone)
    objValUnwrap htpLabelSpadValue(htPage, 'iaone)
  ladifone :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'ladifone)
    objValUnwrap htpLabelSpadValue(htPage, 'ladifone)
  la := 1+n*iaone
--    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'la)
--    objValUnwrap htpLabelSpadValue(htPage, 'la)
  ladif :=1+n*ladifone
--    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'ladif)
--    objValUnwrap htpLabelSpadValue(htPage, 'ladif)
  error := htpButtonValue(htPage,'ifail)
  ifail :=
    error = 'one => '1
    '-1
  (n = '6 and (la ='7 and ladif = '7)) => 
     e02ahfDefaultSolve(htPage,xmin,xmax,iaone,ladifone,ifail)
  labelList :=
    "append"/[f(i) for i in 1..la] where f(i) ==
      prefix := ('"\newline \tab{15} ")
      anam := INTERN STRCONC ('"a",STRINGIMAGE i)
      [['text,:prefix],['bcStrings,[10, "0.0000", anam, 'F]]]
  equationPart := [
     '(domainConditions 
        (isDomain P (Polynomial $EmptyMode))
         (isDomain S (String))
          (isDomain F (Float))
           (isDomain PI (PositiveInteger))),
            :labelList]
  page := htInitPage("E02AHF - Derivative of fitted polynomial in Chebyshev series",nil)
  htSay '"\menuitemstyle{}\tab{2} Coefficients of {\it a(la)}:"
  htMakePage equationPart
  htSay '"\blankline "
  htMakeDoneButton('"Continue",'e02ahfGen)
  htpSetProperty(page,'n,n)
  htpSetProperty(page,'la,la)
  htpSetProperty(page,'ladif,ladif)
  htpSetProperty(page,'xmin,xmin)
  htpSetProperty(page,'xmax,xmax)
  htpSetProperty(page,'iaone,iaone)
  htpSetProperty(page,'ladifone,ladifone)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()


e02ahfDefaultSolve  (htPage,xmin,xmax,iaone,ladifone,ifail) ==
  n := '6
  la := '7
  ladif := '7
  page := htInitPage('"E02AHF - Derivative of fitted polynomial in Chebyshev series",nil)
  htMakePage '(
    (domainConditions 
      (isDomain PI (Positive Integer))
       (isDomain F (Float)))
    (text . "\menuitemstyle{}\tab{2} ")
    (text . "Coefficients of {\it a(la)}: ")
    (text . "\newline \tab{15}")
    (bcStrings (10 "2.53213" a1 F))
    (text . "\newline \tab{15}")
    (bcStrings (10 "1.13032" a2 F))
    (text . "\newline \tab{15}")
    (bcStrings (10 "0.27150" a3 F))
    (text . "\newline \tab{15}")
    (bcStrings (10 "0.04434" a4 F))
    (text . "\newline \tab{15}")
    (bcStrings (10 "0.00547" a5 F))
    (text . "\newline \tab{15}")
    (bcStrings (10 "0.00054" a6 F))
    (text . "\newline \tab{15}")
    (bcStrings (10 "0.00004" a7 F))
    (text . "\blankline"))
  htMakeDoneButton('"Continue",'e02ahfGen)      
  htpSetProperty(page,'n,n)
  htpSetProperty(page,'la,la)
  htpSetProperty(page,'ladif,ladif)
  htpSetProperty(page,'xmin,xmin)
  htpSetProperty(page,'xmax,xmax)
  htpSetProperty(page,'iaone,iaone)
  htpSetProperty(page,'ladifone,ladifone)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()

e02ahfGen htPage ==
  n := htpProperty(htPage,'n)
  la := htpProperty(htPage,'la)
  ladif := htpProperty(htPage,'ladif)
  xmin := htpProperty(htPage,'xmin)
  xmax := htpProperty(htPage,'xmax)
  iaone := htpProperty(htPage,'iaone)
  ladifone := htpProperty(htPage,'ladifone)
  ifail := htpProperty(htPage,'ifail)
  alist := htpInputAreaAlist htPage
  np1 := n + 1
  y := alist
  while y repeat
    right := STRCONC ((first y).1," ")
    y := rest y
    arrayList := [right,:arrayList]
  astring := bcwords2liststring arrayList
  prefix := STRCONC('"e02ahf(",STRINGIMAGE np1,", ",xmin,", ",xmax,", [")
  prefix := STRCONC(prefix,astring,"], ",STRINGIMAGE iaone,", ")
  prefix := STRCONC(prefix,STRINGIMAGE la,", ",STRINGIMAGE ladifone,", ")
  prefix := STRCONC(prefix,STRINGIMAGE ladif,", ",STRINGIMAGE ifail,")")
  linkGen prefix

e02ajf() ==
  htInitPage('"E02AJF - Integral of fitted polynomial in Chebyshev series form",nil)
  htMakePage '(
    (domainConditions
      (isDomain F (Float)) 
       (isDomain PI (PositiveInteger)))
    (text . "\windowlink{Manual Page}{manpageXXe02ajf} for this routine ")
    (text . "\newline ")
    (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|e02ajf| '|NagFittingPackage|)} for this routine")
    (text . "\newline \horizontalline ")
    (text . "Determines the indefinite integral of the Chebyshev series ")
    (text . "representation \newline \center{\htbitmap{e02ahf1}} ")
    (text . "of a polynomial, where \htbitmap{xbar} is the normalis")
    (text . "ed argument, related to the original variable {\it x} by the ")
    (text . "transformation \blankline \center{\htbitmap{e02adf1}}")
    (text . ",\htbitmap{xmin} and \htbitmap{xmax} being ")
    (text . "minimum and maximum values of {\it x} respectively. The integral")
    (text . " polynomial has the form ")
    (text . "\blankline \center{\htbitmap{e02ajf}}")
    (text . "and the integration is with respect to the original variable ")
    (text . "{\it x} \blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} \newline ")
    (text . "Degree of the polynomial {\it n}:")
    (text . "\newline \tab{2} ")
    (bcStrings (6 6 n PI))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2}")
    (text . "\newline \htbitmap{xmin}: ")
    (text . "\tab{32} \menuitemstyle{}\tab{34}")
    (text . "\htbitmap{xmax}: ")
    (text . "\newline\tab{2} ")
    (bcStrings (6 "-0.5" xmin F))
    (text . "\tab{34} ")
    (bcStrings (6 "2.5" xmax F))
--    (text . "\blankline ")
--    (text . "\newline ")
--    (text . "\menuitemstyle{}\tab{2}")
--    (text . "\newline Dimension of array {\it a}, {\it la}: ")
--    (text . "\tab{32} \menuitemstyle{}\tab{34}")
--    (text . "Dimension of {\it aint}, {\it laint}: ")
--    (text . "\newline\tab{2} ")
--    (bcStrings (6 7 la PI))
--    (text . "\tab{34} ")
--    (bcStrings (6 8 laint PI))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2}")
    (text . "\newline Index increment of {\it a}, {\it ia1}: ")
    (text . "\tab{32} \menuitemstyle{}\tab{34}")
    (text . "Increment of {\it aint}, {\it iaint1}: ")
    (text . "\newline \tab{2} ")
    (bcStrings (6 1 iaone PI))
    (text . "\tab{34} ")
    (bcStrings (6 1 iaintone PI))
    (text . "\blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} \newline ")
    (text . "Constant of integration {\it qatm1}: ")
    (text . "\newline \tab{2} ")
    (bcStrings (6 "0.0" qatmone F))
    (text . "\blankline")
    (text . "\newline \menuitemstyle{} \tab{2} ")
    (text . "\newline \tab{2} ")
    (text . "Ifail value: ")
    (radioButtons ifail
        ("" " -1, Print error messages" minusOne)
        ("" "  1, Suppress error messages" one)))
  htMakeDoneButton('"Continue", 'e02ajfSolve)
  htShowPage()

e02ajfSolve htPage ==
  n :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'n)
    objValUnwrap htpLabelSpadValue(htPage, 'n)
  xmin := htpLabelInputString(htPage,'xmin)
  xmax := htpLabelInputString(htPage,'xmax)
  iaone :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'iaone)
    objValUnwrap htpLabelSpadValue(htPage, 'iaone)
  iaintone :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'iaintone)
    objValUnwrap htpLabelSpadValue(htPage, 'iaintone)
  la := 1+n*iaone
--    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'la)
--    objValUnwrap htpLabelSpadValue(htPage, 'la)
  laint := n*iaintone + 1
  qatmone := htpLabelInputString(htPage,'qatmone)
  error := htpButtonValue(htPage,'ifail)
  ifail :=
    error = 'one => '1
    '-1
  (n = '6 and (la ='7 and laint = '7)) => 
     e02ajfDefaultSolve(htPage,xmin,xmax,iaone,iaintone,qatmone,ifail)
  labelList :=
    "append"/[f(i) for i in 1..la] where f(i) ==
      prefix := ('"\newline \tab{15} ")
      anam := INTERN STRCONC ('"a",STRINGIMAGE i)
      [['text,:prefix],['bcStrings,[10, "0.0000", anam, 'F]]]
  equationPart := [
     '(domainConditions 
        (isDomain P (Polynomial $EmptyMode))
         (isDomain S (String))
          (isDomain F (Float))
           (isDomain PI (PositiveInteger))),
            :labelList]
  page := htInitPage("E02AJF - Integral of fitted polynomial in Chebyshev series form",nil)
  htSay '"\menuitemstyle{}\tab{2} Coefficients of {\it a(la)}: "
  htMakePage equationPart
  htSay '"\blankline "
  htMakeDoneButton('"Continue",'e02ajfGen)
  htpSetProperty(page,'n,n)
  htpSetProperty(page,'la,la)
  htpSetProperty(page,'laint,laint)
  htpSetProperty(page,'xmin,xmin)
  htpSetProperty(page,'xmax,xmax)
  htpSetProperty(page,'iaone,iaone)
  htpSetProperty(page,'iaintone,iaintone)
  htpSetProperty(page,'qatmone,qatmone)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()


e02ajfDefaultSolve  (htPage,xmin,xmax,iaone,iaintone,qatmone,ifail) ==
  n := '6
  la := '7
  laint := '8
  page := htInitPage('"E02AJF - Integral of fitted polynomial in Chebyshev series form",nil)
  htMakePage '(
    (domainConditions 
      (isDomain PI (Positive Integer))
       (isDomain F (Float)))
    (text . "\menuitemstyle{}\tab{2} ")
    (text . "Coefficients of {\it a(la)}: ")
    (text . "\newline \tab{15}")
    (bcStrings (10 "2.53213" a1 F))
    (text . "\newline \tab{15}")
    (bcStrings (10 "1.13032" a2 F))
    (text . "\newline \tab{15}")
    (bcStrings (10 "0.27150" a3 F))
    (text . "\newline \tab{15}")
    (bcStrings (10 "0.04434" a4 F))
    (text . "\newline \tab{15}")
    (bcStrings (10 "0.00547" a5 F))
    (text . "\newline \tab{15}")
    (bcStrings (10 "0.00054" a6 F))
    (text . "\newline \tab{15}")
    (bcStrings (10 "0.00004" a7 F))
    (text . "\blankline"))
  htMakeDoneButton('"Continue",'e02ajfGen)      
  htpSetProperty(page,'n,n)
  htpSetProperty(page,'la,la)
  htpSetProperty(page,'laint,laint)
  htpSetProperty(page,'xmin,xmin)
  htpSetProperty(page,'xmax,xmax)
  htpSetProperty(page,'iaone,iaone)
  htpSetProperty(page,'iaintone,iaintone)
  htpSetProperty(page,'qatmone,qatmone)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()

e02ajfGen htPage ==
  n := htpProperty(htPage,'n)
  la := htpProperty(htPage,'la)
  laint := htpProperty(htPage,'laint)
  xmin := htpProperty(htPage,'xmin)
  xmax := htpProperty(htPage,'xmax)
  iaone := htpProperty(htPage,'iaone)
  iaintone := htpProperty(htPage,'iaintone)
  qatmone := htpProperty(htPage,'qatmone)
  ifail := htpProperty(htPage,'ifail)
  alist := htpInputAreaAlist htPage
  np1 := n + 1
  y := alist
  while y repeat
    right := STRCONC ((first y).1," ")
    y := rest y
    arrayList := [right,:arrayList]
  astring := bcwords2liststring arrayList
  prefix := STRCONC('"e02ajf(",STRINGIMAGE np1,", ",xmin,", ",xmax,", [")
  prefix := STRCONC(prefix,astring,"], ",STRINGIMAGE iaone,", ")
  prefix := STRCONC(prefix,STRINGIMAGE la,", ",qatmone,", ")
  prefix := STRCONC(prefix,STRINGIMAGE iaintone)
  prefix := STRCONC(prefix,", ",STRINGIMAGE laint,", ",STRINGIMAGE ifail,")")
  linkGen prefix

e02akf() ==
  htInitPage('"E02AKF - Evaluation of fitted polynomial in one variable, from Chebyshev series form",nil)
  htMakePage '(
    (domainConditions
      (isDomain F (Float)) 
       (isDomain PI (PositiveInteger)))
    (text . "\windowlink{Manual Page}{manpageXXe02akf} for this routine ")
    (text . "\newline ")
    (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|e02akf| '|NagFittingPackage|)} for this routine")
    (text . "\newline \horizontalline ")
    (text . "Evaluates at the point x the Chebyshev series representation ")
    (text . "representation \newline \center{\htbitmap{e02ahf1}} ")
    (text . "of a polynomial, where \htbitmap{xbar} is the normalis")
    (text . "ed argument, related to the original variable {\it x} by the ")
    (text . "transformation \blankline \center{\htbitmap{e02adf1}}")
    (text . ",\htbitmap{xmin} and \htbitmap{xmax} being ")
    (text . "minimum and maximum values of {\it x} respectively. ")
    (text . "\blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} \newline ")
    (text . "Degree of the polynomial {\it n}:")
    (text . "\tab{32} \menuitemstyle{}\tab{34}")
    (text . "Evaluation point {\it x}: ")
    (text . "\newline \tab{2} ")
    (bcStrings (6 6 n PI))
    (text . "\tab{34} ")
    (bcStrings (6 "-0.5" x F))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2}")
    (text . "\newline \htbitmap{xmin}: ")
    (text . "\tab{32} \menuitemstyle{}\tab{34}")
    (text . "\htbitmap{xmax}: ")
    (text . "\newline\tab{2} ")
    (bcStrings (6 "-0.5" xmin F))
    (text . "\tab{34} ")
    (bcStrings (6 "2.5" xmax F))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2}")
--    (text . "\newline Dimension of array {\it a}, {\it la} : ")
--    (text . "\tab{32} \menuitemstyle{}\tab{34}")
    (text . "Index increment of {\it a}, {\it ia1}: ")
    (text . "\newline\tab{2} ")
--    (bcStrings (6 7 la PI))
--    (text . "\tab{34} ")
    (bcStrings (6 1 iaone PI))
    (text . "\blankline")
    (text . "\newline \menuitemstyle{} \tab{2} ")
    (text . "\newline \tab{2} ")
    (text . "Ifail value: ")
    (radioButtons ifail
        ("" " -1, Print error messages" minusOne)
        ("" "  1, Suppress error messages" one)))
  htMakeDoneButton('"Continue", 'e02akfSolve)
  htShowPage()

e02akfSolve htPage ==
  n :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'n)
    objValUnwrap htpLabelSpadValue(htPage, 'n)
  x := htpLabelInputString(htPage,'x)
  xmin := htpLabelInputString(htPage,'xmin)
  xmax := htpLabelInputString(htPage,'xmax)
  iaone :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'iaone)
    objValUnwrap htpLabelSpadValue(htPage, 'iaone)
  la := 1+n*iaone
--    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'la)
--    objValUnwrap htpLabelSpadValue(htPage, 'la)
  error := htpButtonValue(htPage,'ifail)
  ifail :=
    error = 'one => '1
    '-1
  (n = '6 and la ='7) =>  e02akfDefaultSolve(htPage,xmin,xmax,x,iaone,ifail)
  labelList :=
    "append"/[f(i) for i in 1..la] where f(i) ==
      prefix := ('"\newline \tab{15} ")
      anam := INTERN STRCONC ('"a",STRINGIMAGE i)
      [['text,:prefix],['bcStrings,[10, "0.0000", anam, 'F]]]
  equationPart := [
     '(domainConditions 
        (isDomain P (Polynomial $EmptyMode))
         (isDomain S (String))
          (isDomain F (Float))
           (isDomain PI (PositiveInteger))),
            :labelList]
  page := htInitPage("E02AKF - Evaluation of fitted polynomial in one variable, from Chebyshev series form",nil)
  htSay '"\menuitemstyle{}\tab{2} Coefficients of {\it a(la)}:"
  htMakePage equationPart
  htSay '"\blankline "
  htMakeDoneButton('"Continue",'e02akfGen)
  htpSetProperty(page,'n,n)
  htpSetProperty(page,'x,x)
  htpSetProperty(page,'la,la)
  htpSetProperty(page,'xmin,xmin)
  htpSetProperty(page,'xmax,xmax)
  htpSetProperty(page,'iaone,iaone)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()


e02akfDefaultSolve  (htPage,xmin,xmax,x,iaone,ifail) ==
  n := '6
  la := '7
  page := htInitPage('"E02AKF - Evaluation of fitted polynomial in one variable, from Chebyshev series form",nil)
  htMakePage '(
    (domainConditions 
      (isDomain PI (Positive Integer))
       (isDomain F (Float)))
    (text . "\menuitemstyle{}\tab{2} ")
    (text . "Coefficients of {\it a(la)}: ")
    (text . "\newline \tab{15}")
    (bcStrings (10 "2.53213" a1 F))
    (text . "\newline \tab{15}")
    (bcStrings (10 "1.13032" a2 F))
    (text . "\newline \tab{15}")
    (bcStrings (10 "0.27150" a3 F))
    (text . "\newline \tab{15}")
    (bcStrings (10 "0.04434" a4 F))
    (text . "\newline \tab{15}")
    (bcStrings (10 "0.00547" a5 F))
    (text . "\newline \tab{15}")
    (bcStrings (10 "0.00054" a6 F))
    (text . "\newline \tab{15}")
    (bcStrings (10 "0.00004" a7 F))
    (text . "\blankline"))
  htMakeDoneButton('"Continue",'e02akfGen)      
  htpSetProperty(page,'n,n)
  htpSetProperty(page,'la,la)
  htpSetProperty(page,'x,x)
  htpSetProperty(page,'xmin,xmin)
  htpSetProperty(page,'xmax,xmax)
  htpSetProperty(page,'iaone,iaone)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()

e02akfGen htPage ==
  n := htpProperty(htPage,'n)
  x := htpProperty(htPage,'x)
  la := htpProperty(htPage,'la)
  xmin := htpProperty(htPage,'xmin)
  xmax := htpProperty(htPage,'xmax)
  iaone := htpProperty(htPage,'iaone)
  ifail := htpProperty(htPage,'ifail)
  alist := htpInputAreaAlist htPage
  np1 := n + 1
  y := alist
  while y repeat
    right := STRCONC ((first y).1," ")
    y := rest y
    arrayList := [right,:arrayList]
  astring := bcwords2liststring arrayList
  prefix := STRCONC('"e02akf(",STRINGIMAGE np1,", ",xmin,", ",xmax,", [")
  prefix := STRCONC(prefix,astring,"], ",STRINGIMAGE iaone,", ")
  prefix := STRCONC(prefix,STRINGIMAGE la,", ",x,", ",STRINGIMAGE ifail,")")
  linkGen prefix


e02baf() ==
  htInitPage('"E02BAF - Least-squares curve cubic spine fit",nil)
  htMakePage '(
    (domainConditions 
       (isDomain PI (PositiveInteger)))
    (text . "\windowlink{Manual Page}{manpageXXe02baf} for this routine ")
    (text . "\newline ")
    (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|e02baf| '|NagFittingPackage|)} for this routine")
    (text . "\newline \horizontalline ")
    (text . "Determines a least-squares cubic spline approximation to the ")
    (text . "set of points {\it (}\htbitmap{xr},")
    (text . "\htbitmap{yr}{\it )} with weights ")
    (text . "\htbitmap{wr}, for r = 1,2,...,m. ")
    (text . "The knots \htbitmap{lamdai}, for i = 1,2,...,")
    (text . "\htbitmap{ncap}+7, are prescribed by the user. The ")
    (text . "spline is given by the B-spline representation \blankline ")
    (text . "\center{\htbitmap{e02baf}} where ")
    (text . "\htbitmap{ncap} is the number of intervals of the ")
    (text . "spline. \blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} \newline ")
    (text . "Number of data points {\it m}:")
    (text . "\newline \tab{2} ")
    (bcStrings (6 14 m PI))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2}")
    (text . "Number of intervals in the spline \htbitmap{ncap}: ")
    (text . "\newline\tab{2} ")
    (bcStrings (6 5 ncap PI))
    (text . "\blankline")
    (text . "\newline")
    (text . "\menuitemstyle{} \tab{2} ")
    (text . "\newline \tab{2} ")
    (text . "Ifail value: ")
    (radioButtons ifail
        ("" " -1, Print error messages" minusOne)
        ("" "  1, Suppress error messages" one)))
  htMakeDoneButton('"Continue", 'e02bafSolve)
  htShowPage()

e02bafSolve htPage ==
  m :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'm)
    objValUnwrap htpLabelSpadValue(htPage, 'm)
  ncap :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'ncap)
    objValUnwrap htpLabelSpadValue(htPage, 'ncap)
  error := htpButtonValue(htPage,'ifail)
  ifail :=
    error = 'one => '1
    '-1
  (m = '14 and ncap ='5) => e02bafDefaultSolve(htPage,ifail)
  labelList :=
    "append"/[f(i) for i in 1..m] where f(i) ==
      prefix := ('"\newline \tab{2} ")
      middle := ('"\tab{22} ")
      post   := ('" \tab{42} ")
      xnam := INTERN STRCONC ('"x",STRINGIMAGE i)
      ynam := INTERN STRCONC ('"y",STRINGIMAGE i)
      znam := INTERN STRCONC ('"z",STRINGIMAGE i)
      num := INTERN STRCONC (STRINGIMAGE (QUOTIENT(i,10)),".",STRINGIMAGE (REM(i,10)))
      [['text,:prefix],['bcStrings,[10, num, xnam, 'F]], 
       ['text,:middle],['bcStrings,[10, 0.0, ynam, 'F]],
        ['text,:post],['bcStrings,[10, 0.0, znam, 'F]]]
  lamdaList := 
    "append"/[g(j) for j in 5..(ncap+3)] where g(j) ==
       anam := INTERN STRCONC ('"a",STRINGIMAGE j)
       [['bcStrings,[6, 0.0, anam, 'F]]]
  prefix := ('"\blankline \menuitemstyle{}\tab{2} Interior knots ")
  prefix := STRCONC(prefix,"\htbitmap{lamdai}, for i = 5,6,...,")
  prefix := STRCONC(prefix,"\htbitmap{ncap} + 3: \newline \tab{2}" )
  lamdaList := [['text,:prefix],:lamdaList]
  equationPart := [
     '(domainConditions 
        (isDomain P (Polynomial $EmptyMode))
         (isDomain S (String))
          (isDomain F (Float))
           (isDomain PI (PositiveInteger))),
            :labelList,:lamdaList]
  page := htInitPage("E02BAF - Least-squares curve cubic spline fit",nil)
  htSay '"\menuitemstyle{}\tab{2} Values of \space{1} "
  htSay '"\htbitmap{xr}: \tab{20} \menuitemstyle{}\tab{22} "
  htSay '"Values of \space{1} \htbitmap{fr}: \tab{40}"
  htSay '"\menuitemstyle{}\tab{42} Values of \space{1} "
  htSay '"\htbitmap{dr}: "
  htMakePage equationPart
  htSay '"\blankline "
  htMakeDoneButton('"Continue",'e02bafGen)
  htpSetProperty(page,'m,m)
  htpSetProperty(page,'ncap,ncap)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()


e02bafDefaultSolve (htPage,ifail) ==
  m := '14
  ncap := '5
  page := htInitPage('"E02BAF - Least-squares curve cubic spline fit",nil)
  htMakePage '(
    (domainConditions 
       (isDomain F (Float)))
    (text . "\menuitemstyle{}\tab{2} Values of \space{1} ")
    (text . "\htbitmap{xr}: \tab{20} \menuitemstyle{}\tab{22} ")
    (text . "Values of \space{1} \htbitmap{yr}: \tab{40} ")
    (text . "\menuitemstyle{}\tab{42} Values of \space{1} ")
    (text . "\htbitmap{wr}: ")
    (text . "\newline \tab{2}")
    (bcStrings (10 "0.20" x1 F))
    (text . "\tab{22} ")
    (bcStrings (10 "0.00" y1 F))
    (text . "\tab{42} ")
    (bcStrings (10 "0.20" z1 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "0.47" x2 F))
    (text . "\tab{22} ")
    (bcStrings (10 "2.00" y2 F))
    (text . "\tab{42} ")
    (bcStrings (10 "0.20" z2 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "0.74" x3 F))
    (text . "\tab{22} ")
    (bcStrings (10 "4.00" y3 F))
    (text . "\tab{42} ")
    (bcStrings (10 "0.30" z3 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "1.09" x4 F))
    (text . "\tab{22} ")
    (bcStrings (10 "6.00" y4 F))
    (text . "\tab{42} ")
    (bcStrings (10 "0.70" z4 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "1.60" x5 F))
    (text . "\tab{22} ")
    (bcStrings (10 "8.00" y5 F))
    (text . "\tab{42} ")
    (bcStrings (10 "0.90" z5 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "1.90" x6 F))
    (text . "\tab{22} ")
    (bcStrings (10 "8.62" y6 F))
    (text . "\tab{42} ")
    (bcStrings (10 "1.00" z6 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "2.60" x7 F))
    (text . "\tab{22} ")
    (bcStrings (10 "9.10" y7 F))
    (text . "\tab{42} ")
    (bcStrings (10 "1.00" z7 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "3.10" x8 F))
    (text . "\tab{22} ")
    (bcStrings (10 "8.90" y8 F))
    (text . "\tab{42} ")
    (bcStrings (10 "1.00" z8 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "4.00" x9 F))
    (text . "\tab{22} ")
    (bcStrings (10 "8.15" y9 F))
    (text . "\tab{42} ")
    (bcStrings (10 "0.80" z9 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "5.15" x10 F))
    (text . "\tab{22} ")
    (bcStrings (10 "7.00" y10 F))
    (text . "\tab{42} ")
    (bcStrings (10 "0.50" z10 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "6.17" x11 F))
    (text . "\tab{22} ")
    (bcStrings (10 "6.00" y11 F))
    (text . "\tab{42} ")
    (bcStrings (10 "0.70" z11 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "8.00" x12 F))
    (text . "\tab{22} ")
    (bcStrings (10 "4.54" y12 F))
    (text . "\tab{42} ")
    (bcStrings (10 "1.00" z12 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "10.00" x13 F))
    (text . "\tab{22} ")
    (bcStrings (10 "3.39" y13 F))
    (text . "\tab{42} ")
    (bcStrings (10 "1.00" z13 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "12.00" x14 F))
    (text . "\tab{22} ")
    (bcStrings (10 "2.56" y14 F))
    (text . "\tab{42} ")
    (bcStrings (10 "1.00" z14 F))
    (text . "\blankline ")
    (text . "\menuitemstyle{} \tab{2} ")
    (text . "Interior knots \htbitmap{lamdai}, for i = 5,6,...")
    (text . "\htbitmap{ncap} + 3: \newline \tab{2}")
    (bcStrings (6 "1.50" l1 F))
    (bcStrings (6 "2.60" l2 F))
    (bcStrings (6 "4.00" l3 F))
    (bcStrings (6 "8.00" l4 F))
    (text . "\blankline"))
  htMakeDoneButton('"Continue",'e02bafGen)      
  htpSetProperty(page,'m,m)
  htpSetProperty(page,'ncap,ncap)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()

e02bafGen htPage ==
  m := htpProperty(htPage,'m)
  ncap := htpProperty(htPage,'ncap)
  ifail := htpProperty(htPage,'ifail)
  alist := htpInputAreaAlist htPage
  ncap7 := ncap + 7
  y := alist
  for i in (ncap+4)..(ncap+7) repeat
    lambda := STRCONC( "0.0"," ")
    lambdaList := [lambda,:lambdaList]
  for i in 5..(ncap+3) repeat
    lambda := STRCONC ((first y).1," ")
    y := rest y
    lambdaList := [lambda,:lambdaList]
  for i in 1..4 repeat
    lambda := STRCONC( "0.0"," ")
    lambdaList := [lambda,:lambdaList]
  lambdaString := bcwords2liststring lambdaList
  while y repeat
    right := STRCONC ((first y).1," ")
    y := rest y    
    mid :=  STRCONC ((first y).1," ")
    y := rest y
    left :=  STRCONC ((first y).1," ")
    y := rest y
    xlist := [left,:xlist]
    ylist := [mid,:ylist]
    wlist := [right,:wlist]
  xstring := bcwords2liststring xlist
  ystring := bcwords2liststring ylist
  wstring := bcwords2liststring wlist
  prefix := STRCONC('"e02baf(",STRINGIMAGE m,", ",STRINGIMAGE ncap7,", [")
  prefix := STRCONC(prefix,xstring,"],[",ystring,"],[",wstring,"], [")
  prefix := STRCONC(prefix,lambdaString,"], ",STRINGIMAGE ifail,")")
  linkGen prefix


e02bbf() ==
  htInitPage('"E02BBF - Evaluation of fitted cubic spline, function only",nil)
  htMakePage '(
    (domainConditions 
       (isDomain PI (PositiveInteger)))
    (text . "\windowlink{Manual Page}{manpageXXe02bbf} for this routine ")
    (text . "\newline ")
    (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|e02bbf| '|NagFittingPackage|)} for this routine")
    (text . "\newline \horizontalline ")
    (text . "Evaluates at the point {\it x} a cubic spline from its B-spline ")
    (text . "B-spline representation ")
    (text . "\center{\htbitmap{e02baf}} where ")
    (text . "\htbitmap{ncap} is the number of intervals of the ")
    (text . "spline. The spline has knots \htbitmap{lamdai}, for ")
    (text . "i = 1,2,...,\htbitmap{ncap} + 7. \blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} \newline ")
    (text . "Number of intervals in the spline \htbitmap{ncap}: ")
    (text . "\newline\tab{2} ")
    (bcStrings (6 4 ncap PI))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2}")
    (text . "Evaluation point {\it x}:")
    (text . "\newline \tab{2} ")
    (bcStrings (6 "2.0" x F))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{} \tab{2} ")
    (text . "\newline \tab{2} ")
    (text . "Ifail value: ")
    (radioButtons ifail
        ("" " -1, Print error messages" minusOne)
        ("" "  1, Suppress error messages" one)))
  htMakeDoneButton('"Continue", 'e02bbfSolve)
  htShowPage()

e02bbfSolve htPage ==
  ncap :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'ncap)
    objValUnwrap htpLabelSpadValue(htPage, 'ncap)
  x := htpLabelInputString(htPage,'x)
  error := htpButtonValue(htPage,'ifail)
  ifail :=
    error = 'one => '1
    '-1
  ncap = '4  => e02bbfDefaultSolve(htPage,x,ifail)
  labelList :=
    "append"/[f(i) for i in 1..(ncap+7)] where f(i) ==
      prefix := ('"\newline \tab{2} ")
      middle := ('"\tab{22} ")
      lnam := INTERN STRCONC ('"l",STRINGIMAGE i)
      cnam := INTERN STRCONC ('"c",STRINGIMAGE i)
      num := INTERN STRCONC (STRINGIMAGE (QUOTIENT(i,10)),".",STRINGIMAGE (REM(i,10)))
      [['text,:prefix],['bcStrings,[10, num, lnam, 'F]], 
       ['text,:middle],['bcStrings,[10, 0.0, cnam, 'F]]]
  equationPart := [
     '(domainConditions 
        (isDomain P (Polynomial $EmptyMode))
         (isDomain S (String))
          (isDomain F (Float))
           (isDomain PI (PositiveInteger))),
            :labelList]
  page := htInitPage("E02BBF - Evaluation of fitted cubic spline, function only",nil)
  htSay '"\menuitemstyle{}\tab{2} Knots \htbitmap{lamdai}: "
  htSay '"\tab{20} \menuitemstyle{}\tab{22} "
  htSay '"Coefficients \space{1} \htbitmap{ci}: "
  htMakePage equationPart
  htSay '"\blankline "
  htMakeDoneButton('"Continue",'e02bbfGen)
  htpSetProperty(page,'ncap,ncap)
  htpSetProperty(page,'x,x)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()


e02bbfDefaultSolve (htPage,x,ifail) ==
  ncap := '4
  page := htInitPage('"E02BBF - Evaluation of fitted cubic spline, function only",nil)
  htMakePage '(
    (domainConditions 
       (isDomain F (Float)))
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2} Knots \space{1}")
    (text . "\htbitmap{lamdai}: \tab{20} \menuitemstyle{}\tab{22} ")
    (text . "Coefficients \space{1} \htbitmap{ci}: ")
    (text . "\newline \tab{2}")
    (bcStrings (10 "1.00" l1 F))
    (text . "\tab{22} ")
    (bcStrings (10 "1.00" c1 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "1.00" l2 F))
    (text . "\tab{22} ")
    (bcStrings (10 "2.00" c2 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "1.00" l3 F))
    (text . "\tab{22} ")
    (bcStrings (10 "4.00" c3 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "1.00" l4 F))
    (text . "\tab{22} ")
    (bcStrings (10 "7.00" c4 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "3.00" l5 F))
    (text . "\tab{22} ")
    (bcStrings (10 "6.00" c5 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "6.00" l6 F))
    (text . "\tab{22} ")
    (bcStrings (10 "4.00" c6 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "8.00" l7 F))
    (text . "\tab{22} ")
    (bcStrings (10 "3.00" c7 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "9.00" l8 F))
    (text . "\tab{22} ")
    (bcStrings (10 "0.00" c8 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "9.00" l9 F))
    (text . "\tab{22} ")
    (bcStrings (10 "0.00" c9 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "9.00" l10 F))
    (text . "\tab{22} ")
    (bcStrings (10 "0.00" c10 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "9.00" l11 F))
    (text . "\tab{22} ")
    (bcStrings (10 "0.00" c11 F))
    (text . "\blankline"))
  htMakeDoneButton('"Continue",'e02bbfGen)      
  htpSetProperty(page,'ncap,ncap)
  htpSetProperty(page,'x,x)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()

e02bbfGen htPage ==
  ncap := htpProperty(htPage,'ncap)
  x := htpProperty(htPage,'x)
  ifail := htpProperty(htPage,'ifail)
  alist := htpInputAreaAlist htPage
  ncap7 := ncap + 7
  y := alist
  while y repeat
    right := STRCONC ((first y).1," ")
    y := rest y    
    left :=  STRCONC ((first y).1," ")
    y := rest y
    lamlist := [left,:lamlist]
    clist := [right,:clist]
  lamstring := bcwords2liststring lamlist
  cstring := bcwords2liststring clist
  prefix := STRCONC('"e02bbf(",STRINGIMAGE ncap7,", [",lamstring,"],[")
  prefix := STRCONC(prefix,cstring,"], ",x,", ",STRINGIMAGE ifail,")")
  linkGen prefix


e02bcf() ==
  htInitPage('"E02BCF - Evaluation of fitted cubic spline, function and derivatives",nil)
  htMakePage '(
    (domainConditions 
       (isDomain PI (PositiveInteger)))
    (text . "\windowlink{Manual Page}{manpageXXe02bcf} for this routine ")
    (text . "\newline ")
    (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|e02bcf| '|NagFittingPackage|)} for this routine")
    (text . "\newline \horizontalline ")
    (text . "Evaluates at the point {\it x} a cubic spline and its first ")
    (text . "three derivatives from its B-spline representation ")
    (text . "\center{\htbitmap{e02baf}} where ")
    (text . "\htbitmap{ncap} is the number of intervals of the ")
    (text . "spline. The spline has knots \htbitmap{lamdai}, for ")
    (text . "i = 1,2,...,\htbitmap{ncap} + 7. \blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} \newline ")
    (text . "Number of intervals in the spline \htbitmap{ncap}: ")
    (text . "\newline\tab{2} ")
    (bcStrings (6 7 ncap PI))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2}")
    (text . "Evaluation point {\it x}:")
    (text . "\newline \tab{2} ")
    (bcStrings (6 "2.0" x F))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{} \tab{2} ")
    (text . "\newline \tab{2} ")
    (text . "{\it LEFT} specifies whether LH or RH derivatives are required: ")
    (radioButtons deriv
        ("" "  Left-hand derivative" left)
        ("" "  Right-hand derivative" right))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{} \tab{2} ")
    (text . "\newline \tab{2} ")
    (text . "Ifail value: ")
    (radioButtons ifail
        ("" " -1, Print error messages" minusOne)
        ("" "  1, Suppress error messages" one)))
  htMakeDoneButton('"Continue", 'e02bcfSolve)
  htShowPage()

e02bcfSolve htPage ==
  ncap :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'ncap)
    objValUnwrap htpLabelSpadValue(htPage, 'ncap)
  x := htpLabelInputString(htPage,'x)
  temp := htpButtonValue(htPage,'deriv)
  deriv :=
    temp = 'left => '1
    '2
  error := htpButtonValue(htPage,'ifail)
  ifail :=
    error = 'one => '1
    '-1
  ncap = '7  => e02bcfDefaultSolve(htPage,x,deriv,ifail)
  labelList :=
    "append"/[f(i) for i in 1..(ncap+7)] where f(i) ==
      prefix := ('"\newline \tab{2} ")
      middle := ('"\tab{22} ")
      lnam := INTERN STRCONC ('"l",STRINGIMAGE i)
      cnam := INTERN STRCONC ('"c",STRINGIMAGE i)
      num := INTERN STRCONC (STRINGIMAGE (QUOTIENT(i,10)),".",STRINGIMAGE (REM(i,10)))
      [['text,:prefix],['bcStrings,[10, num, lnam, 'F]], 
       ['text,:middle],['bcStrings,[10, 0.0, cnam, 'F]]]
  equationPart := [
     '(domainConditions 
        (isDomain P (Polynomial $EmptyMode))
         (isDomain S (String))
          (isDomain F (Float))
           (isDomain PI (PositiveInteger))),
            :labelList]
  page := htInitPage("E02BCF - Evaluation of fitted cubic spline, function and derivatives",nil)
  htSay '"\menuitemstyle{}\tab{2} Knots \htbitmap{lamdai}: "
  htSay '"\tab{20} \menuitemstyle{}\tab{22} "
  htSay '"Coefficients \space{1} \htbitmap{ci}: "
  htMakePage equationPart
  htSay '"\blankline "
  htMakeDoneButton('"Continue",'e02bcfGen)
  htpSetProperty(page,'ncap,ncap)
  htpSetProperty(page,'x,x)
  htpSetProperty(page,'deriv,deriv)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()


e02bcfDefaultSolve (htPage,x,deriv,ifail) ==
  ncap := '7
  page := htInitPage('"E02BCF - Evaluation of fitted cubic spline, function and derivatives",nil)
  htMakePage '(
    (domainConditions 
       (isDomain F (Float)))
    (text . "\menuitemstyle{}\tab{2} Knots \space{1}")
    (text . "\htbitmap{lamdai}: \tab{20} \menuitemstyle{}\tab{22} ")
    (text . "Coefficients \space{1} \htbitmap{ci}: ")
    (text . "\newline \tab{2}")
    (bcStrings (10 "0.0" l1 F))
    (text . "\tab{22} ")
    (bcStrings (10 "10.00" c1 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "0.00" l2 F))
    (text . "\tab{22} ")
    (bcStrings (10 "12.00" c2 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "0.00" l3 F))
    (text . "\tab{22} ")
    (bcStrings (10 "13.00" c3 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "0.00" l4 F))
    (text . "\tab{22} ")
    (bcStrings (10 "15.00" c4 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "1.00" l5 F))
    (text . "\tab{22} ")
    (bcStrings (10 "22.00" c5 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "3.00" l6 F))
    (text . "\tab{22} ")
    (bcStrings (10 "26.00" c6 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "3.00" l7 F))
    (text . "\tab{22} ")
    (bcStrings (10 "24.00" c7 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "3.00" l8 F))
    (text . "\tab{22} ")
    (bcStrings (10 "18.00" c8 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "4.00" l9 F))
    (text . "\tab{22} ")
    (bcStrings (10 "14.00" c9 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "4.00" l10 F))
    (text . "\tab{22} ")
    (bcStrings (10 "12.00" c10 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "6.00" l11 F))
    (text . "\tab{22} ")
    (bcStrings (10 "0.00" c11 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "6.00" l12 F))
    (text . "\tab{22} ")
    (bcStrings (10 "0.00" c12 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "6.00" l13 F))
    (text . "\tab{22} ")
    (bcStrings (10 "0.00" c13 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "6.00" l14 F))
    (text . "\tab{22} ")
    (bcStrings (10 "0.00" c14 F))
    (text . "\blankline"))
  htMakeDoneButton('"Continue",'e02bcfGen)      
  htpSetProperty(page,'ncap,ncap)
  htpSetProperty(page,'x,x)
  htpSetProperty(page,'deriv,deriv)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()

e02bcfGen htPage ==
  ncap := htpProperty(htPage,'ncap)
  x := htpProperty(htPage,'x)
  deriv := htpProperty(htPage,'deriv)
  ifail := htpProperty(htPage,'ifail)
  alist := htpInputAreaAlist htPage
  ncap7 := ncap + 7
  y := alist
  while y repeat
    right := STRCONC ((first y).1," ")
    y := rest y    
    left :=  STRCONC ((first y).1," ")
    y := rest y
    lamlist := [left,:lamlist]
    clist := [right,:clist]
  lamstring := bcwords2liststring lamlist
  cstring := bcwords2liststring clist
  prefix := STRCONC('"e02bcf(",STRINGIMAGE ncap7,", [",lamstring,"],[")
  prefix := STRCONC(prefix,cstring,"], ",x,", ",STRINGIMAGE deriv)
  prefix := STRCONC(prefix,", ",STRINGIMAGE ifail,")")
  linkGen prefix



e02bdf() ==
  htInitPage('"E02BDF - Evaluation of fitted cubic spline, definite integral",nil)
  htMakePage '(
    (domainConditions 
       (isDomain PI (PositiveInteger)))
    (text . "\windowlink{Manual Page}{manpageXXe02bdf} for this routine ")
    (text . "\newline ")
    (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|e02bdf| '|NagFittingPackage|)} for this routine")
    (text . "\newline \horizontalline ")
    (text . "Evaluates the definite integral of a cubic spline from its ")
    (text . "B-spline representation \center{\htbitmap{e02baf}} ")
    (text . "where \htbitmap{ncap} is the number of intervals of ")
    (text . "the spline. The spline has knots \htbitmap{lamdai}, ")
    (text . "for i = 1,2,...,\htbitmap{ncap} + 7, and the integral ")
    (text . "is evaluated over the range \htbitmap{e02bdf} ")
    (text . "\blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} \newline ")
    (text . "Number of intervals in the spline \htbitmap{ncap}: ")
    (text . "\newline\tab{2} ")
    (bcStrings (6 7 ncap PI))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{} \tab{2} ")
    (text . "\newline \tab{2} ")
    (text . "Ifail value: ")
    (radioButtons ifail
        ("" " -1, Print error messages" minusOne)
        ("" "  1, Suppress error messages" one)))
  htMakeDoneButton('"Continue", 'e02bdfSolve)
  htShowPage()

e02bdfSolve htPage ==
  ncap :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'ncap)
    objValUnwrap htpLabelSpadValue(htPage, 'ncap)
  error := htpButtonValue(htPage,'ifail)
  ifail :=
    error = 'one => '1
    '-1
  ncap = '7 => e02bdfDefaultSolve(htPage,ifail)
  labelList :=
    "append"/[f(i) for i in 1..(ncap+7)] where f(i) ==
      prefix := ('"\newline \tab{2} ")
      middle := ('"\tab{22} ")
      lnam := INTERN STRCONC ('"l",STRINGIMAGE i)
      cnam := INTERN STRCONC ('"c",STRINGIMAGE i)
      num := INTERN STRCONC (STRINGIMAGE (QUOTIENT(i,10)),".",STRINGIMAGE (REM(i,10)))
      [['text,:prefix],['bcStrings,[10, num, lnam, 'F]], 
       ['text,:middle],['bcStrings,[10, 0.0, cnam, 'F]]]
  equationPart := [
     '(domainConditions 
        (isDomain P (Polynomial $EmptyMode))
         (isDomain S (String))
          (isDomain F (Float))
           (isDomain PI (PositiveInteger))),
            :labelList]
  page := htInitPage("E02BDF - Evaluation of fitted cubic spline, definite integral",nil)
  htSay '"\menuitemstyle{}\tab{2} Knots \htbitmap{lamdai}: "
  htSay '"\tab{20} \menuitemstyle{}\tab{22} "
  htSay '"Coefficients \space{1} \htbitmap{ci}: "
  htMakePage equationPart
  htSay '"\blankline "
  htMakeDoneButton('"Continue",'e02bdfGen)
  htpSetProperty(page,'ncap,ncap)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()


e02bdfDefaultSolve(htPage,ifail) ==
  ncap := '7
  page := htInitPage('"E02BDF - Evaluation of fitted cubic spline, definite integral",nil)
  htMakePage '(
    (domainConditions 
       (isDomain EM $EmptyMode)
       (isDomain F (Float))
       (isDomain I (Integer)))
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2} Knots \space{1}")
    (text . "\htbitmap{lamdai}: \tab{20} \menuitemstyle{}\tab{22} ")
    (text . "Coefficients \space{1} \htbitmap{ci}: ")
    (text . "\newline \tab{2}")
    (bcStrings (10 "0.0" l1 F))
    (text . "\tab{22} ")
    (bcStrings (10 "10.00" c1 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "0.00" l2 F))
    (text . "\tab{22} ")
    (bcStrings (10 "12.00" c2 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "0.00" l3 F))
    (text . "\tab{22} ")
    (bcStrings (10 "13.00" c3 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "0.00" l4 F))
    (text . "\tab{22} ")
    (bcStrings (10 "15.00" c4 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "1.00" l5 F))
    (text . "\tab{22} ")
    (bcStrings (10 "22.00" c5 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "3.00" l6 F))
    (text . "\tab{22} ")
    (bcStrings (10 "26.00" c6 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "3.00" l7 F))
    (text . "\tab{22} ")
    (bcStrings (10 "24.00" c7 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "3.00" l8 F))
    (text . "\tab{22} ")
    (bcStrings (10 "18.00" c8 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "4.00" l9 F))
    (text . "\tab{22} ")
    (bcStrings (10 "14.00" c9 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "4.00" l10 F))
    (text . "\tab{22} ")
    (bcStrings (10 "12.00" c10 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "6.00" l11 F))
    (text . "\tab{22} ")
    (bcStrings (10 "0.00" c11 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "6.00" l12 F))
    (text . "\tab{22} ")
    (bcStrings (10 "0.00" c12 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "6.00" l13 F))
    (text . "\tab{22} ")
    (bcStrings (10 "0.00" c13 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "6.00" l14 F))
    (text . "\tab{22} ")
    (bcStrings (10 "0.00" c14 F)))
  htpSetProperty(page,'ncap,ncap)
  htpSetProperty(page,'ifail,ifail)
  htMakeDoneButton('"Continue",'e02bdfGen)      
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()

e02bdfGen htPage ==
  ncap := htpProperty(htPage,'ncap)
  ifail := htpProperty(htPage,'ifail)
  alist := htpInputAreaAlist htPage
  ncap7 := ncap + 7
  y := alist
  while y repeat
    right := STRCONC ((first y).1," ")
    y := rest y    
    left :=  STRCONC ((first y).1," ")
    y := rest y
    lamlist := [left,:lamlist]
    clist := [right,:clist]
  lamstring := bcwords2liststring lamlist
  cstring := bcwords2liststring clist
  prefix := STRCONC('"e02bdf(",STRINGIMAGE ncap7,", [",lamstring,"],[")
  prefix := STRCONC(prefix,cstring,"], ",STRINGIMAGE ifail,")")
  linkGen prefix



e02bef() ==
  htInitPage('"E02BEF - Least-squares curve cubic spline fit, automatic knot placement",nil)
  htMakePage '(
    (domainConditions 
       (isDomain PI (PositiveInteger)))
    (text . "\windowlink{Manual Page}{manpageXXe02bef} for this routine ")
    (text . "\newline ")
    (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|e02bef| '|NagFittingPackage|)} for this routine")
    (text . "\newline \horizontalline ")
    (text . "\newline ")
    (text . "Determines a cubic spline approximation to the set of points ")
    (text . "{\it ( \htbitmap{xr},\htbitmap{yr}) } ")
    (text . "with weights \htbitmap{wr}, for r = 1,2,...,m. ")
    (text . "The knots \htbitmap{lamdai}, for i = 1,2,...,n, ")
    (text . "are chosen by the routine, but a single parameter S must be ")
    (text . "specified to control the trade-off between closeness of fit and ")
    (text . "smoothness of fit. This affects the number of knots required ")
    (text . "by the spline, which is given in the B-spline representation ")
    (text . "\center{\htbitmap{e02bef}}, where n-1 is the number of")
    (text . " intervals of the spline. ")
    (text . "\blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} \newline ")
    (text . "Number of data points {\it m}: ")
    (text . "\newline\tab{2} ")
    (bcStrings (6 15 m PI))
    (text . "\blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} \newline ")
    (text . "Smoothing factor {\it s}: ")
    (text . "\newline\tab{2} ")
    (bcStrings (6 "1.0" s F))
    (text . "\blankline ")
    (text . "\newline \menuitemstyle{}\tab{2}  ")
    (text . "Over-estimate of number n of knots {\it nest}:\newline\tab{2} ")
    (bcStrings (6 54 nest PI))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{} \tab{2} ")
    (text . "\newline \tab{2} ")
    (text . "Start value: ")
    (radioButtons start
        ("" "  Cold Start - no values needed for {\it n,lamda,wrk} or {\it iwrk}" cold)
        ("" "  Warm Start - uses knots found in a previous call" warm))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{} \tab{2} ")
    (text . "\newline \tab{2} ")
    (text . "Ifail value: ")
    (radioButtons ifail
        ("" " -1, Print error messages" minusOne)
        ("" "  1, Suppress error messages" one)))
  htMakeDoneButton('"Continue", 'e02befSolve)
  htShowPage()

e02befSolve htPage ==
  m :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'm)
    objValUnwrap htpLabelSpadValue(htPage, 'm)
  nest :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'nest)
    objValUnwrap htpLabelSpadValue(htPage, 'nest)
  lwrk := 4*m +16*nest + 41
  s := htpLabelInputString(htPage,'s)
  initial := htpButtonValue(htPage,'start)
  start :=
    initial = 'cold => '1
    '2
  error := htpButtonValue(htPage,'ifail)
  ifail :=
    error = 'one => '1
    '-1
  (m = 15 and start = 1) => e02befDefaultSolve (htPage,nest,lwrk,s,ifail)
  start = 1 => e02befColdSolve (htPage,m,nest,lwrk,s,ifail)
  -- warm start not really possible from hyperdoc
  -- as inputing a workspace array of dimension 1105 is asking too much
  -- user should use the command line, using the previous calculated 
  -- parameters
  htInitPage('"E02BEF - Least-squares curve cubic spline fit, automatic knot placement",nil)
  htMakePage '(
    (domainConditions 
       (isDomain PI (PositiveInteger)))
    (text . "\blankline ")
    (text . "{\center{\it Hyperdoc interface not available for warm starts.}}")
    (text . "\newline ")
    (text . "{\center{\it Please use the command line.}}"))
  htMakeDoneButton('"Continue",'e02bef)
  htShowPage()



e02befColdSolve(htPage,m,nest,lwrk,s,ifail) ==
  labelList :=
    "append"/[f(i) for i in 1..m] where f(i) ==
      prefix := ('"\newline \tab{2} ")
      middle := ('"\tab{22} ")
      post   := ('" \tab{42} ")
      xnam := INTERN STRCONC ('"x",STRINGIMAGE i)
      ynam := INTERN STRCONC ('"y",STRINGIMAGE i)
      znam := INTERN STRCONC ('"z",STRINGIMAGE i)
      num := INTERN STRCONC (STRINGIMAGE (QUOTIENT(i,10)),".",STRINGIMAGE (REM(i,10)))
      [['text,:prefix],['bcStrings,[10, num, xnam, 'F]], 
       ['text,:middle],['bcStrings,[10, 0.0, ynam, 'F]],
        ['text,:post],['bcStrings,[10, 0.0, znam, 'F]]]
  equationPart := [
     '(domainConditions 
        (isDomain P (Polynomial $EmptyMode))
         (isDomain S (String))
          (isDomain F (Float))
           (isDomain PI (PositiveInteger))),
            :labelList]
  page := htInitPage('"E02BEF - Least-squares curve cubic spline fit, automatic knot placement",nil)
  htSay '"\menuitemstyle{}\tab{2} Values of \space{1} "
  htSay '"\htbitmap{xr}: \tab{20} \menuitemstyle{}\tab{22} "
  htSay '"Values of \space{1} \htbitmap{yr}: \tab{40}"
  htSay '"\menuitemstyle{}\tab{42} Values of \space{1} "
  htSay '"\htbitmap{wr}: "
  htMakePage equationPart
  htSay '"\blankline "
  htMakeDoneButton('"Continue",'e02befColdGen)
  htpSetProperty(page,'m,m)
  htpSetProperty(page,'nest,nest)
  htpSetProperty(page,'lwrk,lwrk)
  htpSetProperty(page,'s,s)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()


e02befDefaultSolve (htPage,nest,lwrk,s,ifail) ==
  m := 15
  page := htInitPage('"E02BEF - Least-squares curve cubic spline fit, automatic knot placement",nil)
  htMakePage '(
    (domainConditions 
       (isDomain F (Float)))
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2} Values of \space{1} ")
    (text . "\htbitmap{xr}: \tab{20} \menuitemstyle{}\tab{22} ")
    (text . "Values of \space{1} \htbitmap{yr}: \tab{40} ")
    (text . "\menuitemstyle{}\tab{42} Values of \space{1} ")
    (text . "\htbitmap{wr}: ")
    (text . "\newline \tab{2}")
    (bcStrings (10 "0.00" x1 F))
    (text . "\tab{22} ")
    (bcStrings (10 "-1.1" y1 F))
    (text . "\tab{42} ")
    (bcStrings (10 "1.00" z1 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "0.50" x2 F))
    (text . "\tab{22} ")
    (bcStrings (10 "-0.372" y2 F))
    (text . "\tab{42} ")
    (bcStrings (10 "2.00" z2 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "1.00" x3 F))
    (text . "\tab{22} ")
    (bcStrings (10 "0.431" y3 F))
    (text . "\tab{42} ")
    (bcStrings (10 "1.50" z3 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "1.50" x4 F))
    (text . "\tab{22} ")
    (bcStrings (10 "1.69" y4 F))
    (text . "\tab{42} ")
    (bcStrings (10 "1.00" z4 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "2.00" x5 F))
    (text . "\tab{22} ")
    (bcStrings (10 "2.11" y5 F))
    (text . "\tab{42} ")
    (bcStrings (10 "3.00" z5 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "2.50" x6 F))
    (text . "\tab{22} ")
    (bcStrings (10 "3.10" y6 F))
    (text . "\tab{42} ")
    (bcStrings (10 "1.00" z6 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "3.00" x7 F))
    (text . "\tab{22} ")
    (bcStrings (10 "4.23" y7 F))
    (text . "\tab{42} ")
    (bcStrings (10 "0.50" z7 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "4.00" x8 F))
    (text . "\tab{22} ")
    (bcStrings (10 "4.35" y8 F))
    (text . "\tab{42} ")
    (bcStrings (10 "1.00" z8 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "4.50" x9 F))
    (text . "\tab{22} ")
    (bcStrings (10 "4.81" y9 F))
    (text . "\tab{42} ")
    (bcStrings (10 "2.00" z9 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "5.00" x10 F))
    (text . "\tab{22} ")
    (bcStrings (10 "4.61" y10 F))
    (text . "\tab{42} ")
    (bcStrings (10 "2.50" z10 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "5.50" x11 F))
    (text . "\tab{22} ")
    (bcStrings (10 "4.79" y11 F))
    (text . "\tab{42} ")
    (bcStrings (10 "1.00" z11 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "6.00" x12 F))
    (text . "\tab{22} ")
    (bcStrings (10 "5.23" y12 F))
    (text . "\tab{42} ")
    (bcStrings (10 "3.00" z12 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "7.00" x13 F))
    (text . "\tab{22} ")
    (bcStrings (10 "6.35" y13 F))
    (text . "\tab{42} ")
    (bcStrings (10 "1.00" z13 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "7.50" x14 F))
    (text . "\tab{22} ")
    (bcStrings (10 "7.19" y14 F))
    (text . "\tab{42} ")
    (bcStrings (10 "2.00" z14 F))
    (text . "\newline \tab{2} ")
    (bcStrings (10 "8.00" x15 F))
    (text . "\tab{22} ")
    (bcStrings (10 "7.97" y15 F))
    (text . "\tab{42} ")
    (bcStrings (10 "1.00" z15 F))
    (text . "\blankline"))
  htMakeDoneButton('"Continue",'e02befColdGen)          
  htpSetProperty(page,'m,m)
  htpSetProperty(page,'nest,nest)
  htpSetProperty(page,'lwrk,lwrk)
  htpSetProperty(page,'s,s)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()

e02befColdGen htPage ==
  m := htpProperty(htPage,'m)
  nest := htpProperty(htPage,'nest)
  lwrk := htpProperty(htPage,'lwrk)
  s := htpProperty(htPage,'s)
  cold := '"c"
  ifail := htpProperty(htPage,'ifail)
  alist := htpInputAreaAlist htPage
  y := alist
  while y repeat
    right := STRCONC ((first y).1," ")
    y := rest y    
    mid :=  STRCONC ((first y).1," ")
    y := rest y
    left :=  STRCONC ((first y).1," ")
    y := rest y
    xlist := [left,:xlist]
    ylist := [mid,:ylist]
    wlist := [right,:wlist]
  xstring := bcwords2liststring xlist
  ystring := bcwords2liststring ylist
  wstring := bcwords2liststring wlist
  -- additional entries needed to get it running
  -- but as Start = c they are not used
  -- mmax := 50
  -- nest := mmax + 4 (54)
  -- lwrk := 4*mmax + 16*nest+41 (1105)
  prefix := STRCONC('"e02bef(_"",cold,"_",",STRINGIMAGE m,", [",xstring,"],[")
  prefix := STRCONC(prefix,ystring,"],[",wstring,"], ",STRINGIMAGE s,", ")
  prefix := STRCONC(prefix,STRINGIMAGE nest,", ",STRINGIMAGE lwrk)
--  prefix := STRCONC(prefix,",0, [[0.0 for i in 1..",STRINGIMAGE nest,"]],")
--  prefix := STRCONC(prefix,STRINGIMAGE ifail,", [[0.0 for i in 1..")
--  prefix := STRCONC(prefix,STRINGIMAGE lwrk,"]], [[0 for i in 1..")
--  prefix := STRCONC(prefix,STRINGIMAGE nest,"]] :: Matrix Integer)")
  prefix := STRCONC(prefix,",0, new(1,",STRINGIMAGE nest,",0.0)$Matrix DoubleFloat,")
  prefix := STRCONC(prefix,STRINGIMAGE ifail,", new(1,",STRINGIMAGE lwrk,",0.0)$Matrix DoubleFloat, ")
  prefix := STRCONC(prefix," new(1,",STRINGIMAGE nest,",0)$Matrix Integer)")
  linkGen prefix

e02def() == 
  htInitPage('"E02DEF - Evaluation of a fitted bicubic spline at a vector of points",nil)
  htMakePage '(
    (domainConditions 
       (isDomain PI (PositiveInteger)))
    (text . "\windowlink{Manual Page}{manpageXXe02def} for this routine ")
    (text . "\newline ")
    (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|e02def| '|NagFittingPackage|)} for this routine")
    (text . "\newline \horizontalline ")
    (text . "Evaluates a bicubic spline at the (\htbitmap{xr},")
    (text . "\htbitmap{yr}), for r = 1,2,...,m, from its B-spline ")
    (text . "representation \htbitmap{e02daf} ")
    (text . "\blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} \newline ")
    (text . "Number of evaluation points, {\it m}: ")
    (text . "\newline\tab{2} ")
    (bcStrings (6 7 m PI))
    (text . "\blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} \newline ")
    (text . "Number of (interior & exterior) knots ")
    (text . "\lambda, \htbitmap{px} \htbitmap{great=} 8: \newline\tab{2} ")
    (bcStrings (6 11 px PI))
    (text . "\blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} \newline ")
    (text . "Number of (interior & exterior) knots ")
    (text . "\mu, \htbitmap{py} \htbitmap{great=} 8: \newline\tab{2} ")
    (bcStrings (6 10 py PI))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{} \tab{2} ")
    (text . "\newline \tab{2} ")
    (text . "Ifail value: ")
    (radioButtons ifail
        ("" " -1, Print error messages" minusOne)
        ("" "  1, Suppress error messages" one)))
  htMakeDoneButton('"Continue", 'e02defSolve)
  htShowPage()

e02defSolve htPage ==
  m :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'm)
    objValUnwrap htpLabelSpadValue(htPage, 'm)
  px :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'px)
    objValUnwrap htpLabelSpadValue(htPage, 'px)
  py :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'py)
    objValUnwrap htpLabelSpadValue(htPage, 'py)
  error := htpButtonValue(htPage,'ifail)
  ifail :=
    error = 'one => '1
    '-1
  ((m = '7 and px = '11) and py = '10) => e02defDefaultSolve(htPage,ifail)
  labelList :=
    "append"/[fxy(i) for i in 1..m] where fxy(i) ==
      prefix := ('"\newline \tab{2} ")
      middle := ('"\tab{22} ")
      xnam := INTERN STRCONC ('"x",STRINGIMAGE i)
      ynam := INTERN STRCONC ('"y",STRINGIMAGE i)
      [['text,:prefix],['bcStrings,[8, 0.0, xnam, 'F]], 
       ['text,:middle],['bcStrings,[8, 0.0, ynam, 'F]]]
  lamList :=
    "append"/[flam(i) for i in 1..px] where flam(i) ==
      lnam := INTERN STRCONC ('"l",STRINGIMAGE i)
      [['bcStrings,[8, 0.0, lnam, 'F]]]
  prefix := ('"\blankline \menuitemstyle{} \tab{2} {\it \lambda(nxest)}: \newline")
  lamList := [['text,:prefix],:lamList]
  muList :=
    "append"/[fmu(i) for i in 1..(py)] where fmu(i) ==
      mnam := INTERN STRCONC ('"m",STRINGIMAGE i)
      [['bcStrings,[8, 0.0, mnam, 'F]]]
  prefix := ('"\blankline \menuitemstyle{} \tab{2} {\it \mu(nyest)}:")
  prefix := STRCONC(prefix,"\newline ")
  muList := [['text,:prefix],:muList]
  cList :=
    "append"/[fp(i) for i in 1..((px-4)*(py-4))] where fp(i) ==
      pnam := INTERN STRCONC ('"p",STRINGIMAGE i)
      [['bcStrings,[8, 0.0, pnam, 'F]]]
  prefix := ('"\blankline \menuitemstyle{}\tab{2}Enter values of ")
  prefix := STRCONC(prefix,"{\it c((nxest*4)-(nyest*4))}: \newline ")
  cList := [['text,:prefix],:cList]
  equationPart := [
     '(domainConditions 
        (isDomain P (Polynomial $EmptyMode))
         (isDomain S (String))
          (isDomain F (Float))
           (isDomain PI (PositiveInteger))),
            :labelList,:lamList,:muList,:cList]
  page := htInitPage('"E02DEF - Evaluation of a fitted bicubic spline at a vector of points",nil)
  htSay '"\menuitemstyle{}\tab{2} Values of \htbitmap{xr}: "
  htSay '"\tab{20} \menuitemstyle{}\tab{22} Values of \htbitmap{yr}:"
  htMakePage equationPart
  htSay '"\blankline "
  htMakeDoneButton('"Continue",'e02defGen)
  htpSetProperty(page,'m,m)
  htpSetProperty(page,'px,px)
  htpSetProperty(page,'py,py)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()


e02defDefaultSolve (htPage,ifail) ==
  m := '7
  px := '11
  py := '10
  page := htInitPage('"E02DEF - Evaluation of a fitted bicubic spline at a vector of points",nil)
  htMakePage '(
    (domainConditions 
       (isDomain F (Float)))
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2} Values of \htbitmap{xr}:")
    (text . "\tab{20} \menuitemstyle{} \tab{22} Values of ")
    (text . "\htbitmap{yr}: ")
    (text . "\newline \tab{2} ")
    (bcStrings (8 "1" x1 F))  
    (text . "\tab{22}")
    (bcStrings (8 "0" y1 F))  
    (text . "\newline \tab{2} ")
    (bcStrings (8 "1.1" x2 F))  
    (text . "\tab{22}")
    (bcStrings (8 "0.1" y2 F))  
    (text . "\newline \tab{2} ")
    (bcStrings (8 "1.5" x3 F))  
    (text . "\tab{22}")
    (bcStrings (8 "0.7" y3 F))  
    (text . "\newline \tab{2} ")
    (bcStrings (8 "1.6" x4 F))  
    (text . "\tab{22}")
    (bcStrings (8 "0.4" y4 F))  
    (text . "\newline \tab{2} ")
    (bcStrings (8 "1.9" x5 F))  
    (text . "\tab{22}")
    (bcStrings (8 "0.3" y5 F))  
    (text . "\newline \tab{2} ")
    (bcStrings (8 "1.9" x6 F))  
    (text . "\tab{22}")
    (bcStrings (8 "0.8" y6 F))  
    (text . "\newline \tab{2} ")
    (bcStrings (8 "2" x7 F))  
    (text . "\tab{22}")
    (bcStrings (8 "1" y7 F))  
    (text . "\blankline ")
    (text . "\menuitemstyle{}\tab{2} {\it \lambda(nxest)}:")
    (text . "\newline ")
    (bcStrings (8 "1.0" l1 F))
    (bcStrings (8 "1.0" l2 F))
    (bcStrings (8 "1.0" l3 F))
    (bcStrings (8 "1.0" l4 F))
    (bcStrings (8 "1.3" l5 F))
    (bcStrings (8 "1.5" l6 F))
    (bcStrings (8 "1.6" l7 F))
    (bcStrings (8 "2" l8 F))
    (bcStrings (8 "2" l9 F))
    (bcStrings (8 "2" l10 F))
    (bcStrings (8 "2" l11 F))
    (text . "\blankline ")
    (text . "\menuitemstyle{}\tab{2} {\it \mu(nyest)}:")
    (text . "\newline ")
    (bcStrings (8 "0" mu1 F))
    (bcStrings (8 "0" mu2 F))
    (bcStrings (8 "0" mu3 F))
    (bcStrings (8 "0" mu4 F))
    (bcStrings (8 "0.4" mu5 F))
    (bcStrings (8 "0.7" mu6 F))
    (bcStrings (8 "1" mu7 F))
    (bcStrings (8 "1" mu8 F))
    (bcStrings (8 "1" mu9 F))
    (bcStrings (8 "1" mu10 F))
    (text . "\blankline \menuitemstyle{}\tab{2} ")
    (text . "Enter values for {\it c((nxest-4)*(nyest-4))}:")
    (text . "\newline ")
    (bcStrings (8 "1" c1 F))    
    (bcStrings (8 "1.1333" c2 F))       
    (bcStrings (8 "1.3667" c3 F))       
    (bcStrings (8 "1.7" c4 F))  
    (bcStrings (8 "1.9" c5 F))  
    (bcStrings (8 "2" c6 F))    
    (bcStrings (8 "1.2" c7 F))  
    (bcStrings (8 "1.3333" c8 F))       
    (bcStrings (8 "1.5667" c9 F))       
    (bcStrings (8 "1.9" c10 F)) 
    (bcStrings (8 "2.1" c11 F)) 
    (bcStrings (8 "2.2" c12 F)) 
    (bcStrings (8 "1.5833" c13 F))      
    (bcStrings (8 "1.7167" c14 F))      
    (bcStrings (8 "1.95" c15 F))        
    (bcStrings (8 "2.2833" c16 F))      
    (bcStrings (8 "2.4833" c17 F))      
    (bcStrings (8 "2.5833" c18 F))      
    (bcStrings (8 "2.1433" c19 F))      
    (bcStrings (8 "2.2767" c20 F))      
    (bcStrings (8 "2.51" c21 F))        
    (bcStrings (8 "2.8433" c22 F))      
    (bcStrings (8 "3.0433" c23 F))      
    (bcStrings (8 "3.1433" c24 F))      
    (bcStrings (8 "2.8667" c25 F))      
    (bcStrings (8 "3" c26 F))   
    (bcStrings (8 "3.2333" c27 F))      
    (bcStrings (8 "3.5667" c28 F))      
    (bcStrings (8 "3.7667" c29 F))      
    (bcStrings (8 "3.8667" c30 F))      
    (bcStrings (8 "3.4667" c31 F))      
    (bcStrings (8 "3.6" c32 F)) 
    (bcStrings (8 "3.8333" c33 F))      
    (bcStrings (8 "4.1667" c34 F))      
    (bcStrings (8 "4.3667" c35 F))      
    (bcStrings (8 "4.4667" c36 F))      
    (bcStrings (8 "4" c37 F))   
    (bcStrings (8 "4.1333" c38 F))      
    (bcStrings (8 "4.3667" c39 F))      
    (bcStrings (8 "4.7" c40 F)) 
    (bcStrings (8 "4.9" c41 F)) 
    (bcStrings (8 "5" c42 F))   
    (text . "\blankline "))
  htMakeDoneButton('"Continue",'e02defGen)      
  htpSetProperty(page,'m,m)
  htpSetProperty(page,'px,px)
  htpSetProperty(page,'py,py)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()

e02defGen htPage ==
  m := htpProperty(htPage,'m)
  px := htpProperty(htPage,'px)
  py := htpProperty(htPage,'py)
  ifail := htpProperty(htPage,'ifail)
  alist := htpInputAreaAlist htPage
  y := alist
  -- c
  for i in 1..((px-4)*(py-4)) repeat
    right := STRCONC ((first y).1," ")
    y := rest y
    cList := [right,:cList]
  cstring := bcwords2liststring cList
  -- mu 
  for i in 1..py repeat
    right := STRCONC ((first y).1," ")
    y := rest y    
    muList := [right,:muList]
  mustring := bcwords2liststring muList
  -- lamda
  for i in 1..px repeat
    right := STRCONC ((first y).1," ")
    y := rest y    
    lamList := [right,:lamList]
  lamstring := bcwords2liststring lamList
  -- x & y
  while y repeat
    one := STRCONC((first y).1," ")
    y := rest y
    two := STRCONC((first y).1," ")
    y := rest y
    xlist := [two,:xlist]
    ylist := [one,:ylist]
  xstring := bcwords2liststring xlist
  ystring := bcwords2liststring ylist
  prefix := STRCONC('"e02def(",STRINGIMAGE m,", ",STRINGIMAGE px,", ")
  prefix := STRCONC(prefix,STRINGIMAGE py,",[",xstring,"],[",ystring,"],[")
  prefix := STRCONC(prefix,lamstring,"],[",mustring,"],[",cstring,"],")
  prefix := STRCONC(prefix,STRINGIMAGE ifail,")")
  linkGen prefix


e02dff() == 
  htInitPage('"E02DFF - Evaluation of a fitted bicubic spline at a mesh of points",nil)
  htMakePage '(
    (domainConditions 
       (isDomain PI (PositiveInteger)))
    (text . "\windowlink{Manual Page}{manpageXXe02dff} for this routine ")
    (text . "\newline ")
    (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|e02dff| '|NagFittingPackage|)} for this routine")
    (text . "\newline \horizontalline ")
    (text . "Evaluates a bicubic spline at all the points on a rectangular ")
    (text . "grid defined by \htbitmap{mx} points ")
    (text . "\htbitmap{xq}on the x-axis and \htbitmap{my}")
    (text . "points \htbitmap{yr} on the y-axis, from its B-spline ")
    (text . "representation \center{\htbitmap{e02daf}} \newline with knot sets ")
    (text . "\{\lambda\} and \{\mu\}. ")
    (text . "\blankline ")
    (text . "\menuitemstyle{}\tab{2}")
    (text . "\newline Grid points on x-axis \htbitmap{mx}: ")
    (text . "\tab{32} \menuitemstyle{}\tab{34}")
    (text . "Grid points on y-axis \htbitmap{my}: ")
    (text . "\newline\tab{2} ")
    (bcStrings (6 7 mx PI))
    (text . "\tab{34} ")
    (bcStrings (6 6 my PI))
    (text . "\blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} ")
    (text . "Number of (interior & exterior) knots \lambda, ")
    (text . "\htbitmap{px} \htbitmap{great=} 8: \newline\tab{2} ")
    (bcStrings (6 11 px PI))
    (text . "\blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} ")
    (text . "Number of (interior & exterior) knots \mu, ")
    (text . "\htbitmap{py} \htbitmap{great=} 8: \newline\tab{2} ")
    (bcStrings (6 10 py PI))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{} \tab{2} ")
    (text . "\newline \tab{2} ")
    (text . "Ifail value: ")
    (radioButtons ifail
        ("" " -1, Print error messages" minusOne)
        ("" "  1, Suppress error messages" one)))
  htMakeDoneButton('"Continue", 'e02dffSolve)
  htShowPage()

e02dffSolve htPage ==
  mx :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'mx)
    objValUnwrap htpLabelSpadValue(htPage, 'mx)
  my :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'my)
    objValUnwrap htpLabelSpadValue(htPage, 'my)
  px :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'px)
    objValUnwrap htpLabelSpadValue(htPage, 'px)
  py :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'py)
    objValUnwrap htpLabelSpadValue(htPage, 'py)
  nwrk1 := 4*mx + px
  nwrk2 := 4*my + py
  nwrklist := [nwrk1,nwrk2]
  nwrkmin := APPLY ('MIN, nwrklist)
  lwrk := nwrkmin
  liwrk :=
    nwrkmin = nwrk2 => my + py -4
    mx + px -4
  error := htpButtonValue(htPage,'ifail)
  ifail :=
    error = 'one => '1
    '-1
  ((mx = '7 and my = '6) and (px = '11 and py = '10)) => 
    e02dffDefaultSolve(htPage,lwrk,liwrk,ifail)
  xList :=
    "append"/[fx(i) for i in 1..mx] where fx(i) ==
      xnam := INTERN STRCONC ('"x",STRINGIMAGE i)
      [['bcStrings,[8, 0.0, xnam, 'F]]]
  yList :=
    "append"/[fy(i) for i in 1..my] where fy(i) ==
      ynam := INTERN STRCONC ('"y",STRINGIMAGE i)
      [['bcStrings,[8, 0.0, ynam, 'F]]]
  prefix := ('"\blankline \menuitemstyle{} \tab{2} Enter values of ")
  prefix := STRCONC(prefix,"\htbitmap{yr} : \newline")
  yList := [['text,:prefix],:yList]
  lamList :=
    "append"/[flam(i) for i in 1..px] where flam(i) ==
      lnam := INTERN STRCONC ('"l",STRINGIMAGE i)
      [['bcStrings,[8, 0.0, lnam, 'F]]]
  prefix := ('"\blankline \menuitemstyle{}\tab{2} {\it \lambda(nxest)}:\newline")
  lamList := [['text,:prefix],:lamList]
  muList :=
    "append"/[fmu(i) for i in 1..(py)] where fmu(i) ==
      mnam := INTERN STRCONC ('"m",STRINGIMAGE i)
      [['bcStrings,[8, 0.0, mnam, 'F]]]
  prefix := ('"\blankline \menuitemstyle{} \tab{2} {\it mu(nyest)}:")
  prefix := STRCONC(prefix,"\newline ")
  muList := [['text,:prefix],:muList]
  cList :=
    "append"/[fp(i) for i in 1..((px-4)*(py-4))] where fp(i) ==
      pnam := INTERN STRCONC ('"p",STRINGIMAGE i)
      [['bcStrings,[8, 0.0, pnam, 'F]]]
  prefix := ('"\blankline \menuitemstyle{}\tab{2} Enter values of ")
  prefix := STRCONC(prefix,"{\it c((px-4)*(py-4))}: \newline")
  cList := [['text,:prefix],:cList]
  equationPart := [
     '(domainConditions 
        (isDomain P (Polynomial $EmptyMode))
         (isDomain S (String))
          (isDomain F (Float))
           (isDomain PI (PositiveInteger))),
            :xList,:yList,:lamList,:muList,:cList]
  page := htInitPage('"E02DFF - Evaluation of a fitted bicubic spline at a mesh of points",nil)
  htSay '"\menuitemstyle{}\tab{2} Values of \htbitmap{xr}:\newline "
  htMakePage equationPart
  htSay '"\blankline "
  htMakeDoneButton('"Continue",'e02dffGen)
  htpSetProperty(page,'mx,mx)
  htpSetProperty(page,'my,my)
  htpSetProperty(page,'px,px)
  htpSetProperty(page,'py,py)
  htpSetProperty(page,'lwrk,lwrk)
  htpSetProperty(page,'liwrk,liwrk)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()


e02dffDefaultSolve (htPage,lwrk,liwrk,ifail) ==
  mx := '7
  my := '6
  px := '11
  py := '10
  page := htInitPage('"E02DFF - Evaluation of a fitted bicubic spline at a mesh of points",nil)
  htMakePage '(
    (domainConditions 
       (isDomain F (Float)))
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2} Enter values of \htbitmap{xr}:")
    (text . "\newline ")
    (bcStrings (8 "1" x1 F))  
    (bcStrings (8 "1.1" x2 F))  
    (bcStrings (8 "1.3" x3 F))  
    (bcStrings (8 "1.4" x4 F))  
    (bcStrings (8 "1.5" x5 F))  
    (bcStrings (8 "1.7" x6 F))  
    (bcStrings (8 "2" x7 F))  
    (text . "\blankline ")
    (text . "\menuitemstyle{}\tab{2} Enter values of \htbitmap{yr}:")
    (text . "\newline ")
    (bcStrings (8 "0" y1 F))  
    (bcStrings (8 "0.2" y2 F))  
    (bcStrings (8 "0.4" y3 F))  
    (bcStrings (8 "0.6" y4 F))  
    (bcStrings (8 "0.8" y5 F))  
    (bcStrings (8 "1" y6 F))  
    (text . "\blankline ")
    (text . "\menuitemstyle{}\tab{2} Enter values of {\it \lambda(nxest)}:")
    (text . "\newline ")
    (bcStrings (8 "1" l1 F))
    (bcStrings (8 "1" l2 F))
    (bcStrings (8 "1" l3 F))
    (bcStrings (8 "1" l4 F))
    (bcStrings (8 "1.3" l5 F))
    (bcStrings (8 "1.5" l6 F))
    (bcStrings (8 "1.6" l7 F))
    (bcStrings (8 "2" l8 F))
    (bcStrings (8 "2" l9 F))
    (bcStrings (8 "2" l10 F))
    (bcStrings (8 "2" l11 F))
    (text . "\blankline ")
    (text . "\menuitemstyle{}\tab{2} Enter values of {\it \mu(nyest)}:")
    (text . "\newline ")
    (bcStrings (8 "0" mu1 F))
    (bcStrings (8 "0" mu2 F))
    (bcStrings (8 "0" mu3 F))
    (bcStrings (8 "0" mu4 F))
    (bcStrings (8 "0.4" mu5 F))
    (bcStrings (8 "0.7" mu6 F))
    (bcStrings (8 "1" mu7 F))
    (bcStrings (8 "1" mu8 F))
    (bcStrings (8 "1" mu9 F))
    (bcStrings (8 "1" mu10 F))
    (text . "\blankline ")
    (text . "\menuitemstyle{}\tab{2} Enter values of {\it c((px-4)*(py-4))}:")
    (text . "\newline ")
    (bcStrings (8 "1" c1 F))    
    (bcStrings (8 "1.1333" c2 F))       
    (bcStrings (8 "1.3667" c3 F))       
    (bcStrings (8 "1.7" c4 F))  
    (bcStrings (8 "1.9" c5 F))  
    (bcStrings (8 "2" c6 F))    
    (bcStrings (8 "1.2" c7 F))  
    (bcStrings (8 "1.3333" c8 F))       
    (bcStrings (8 "1.5667" c9 F))       
    (bcStrings (8 "1.9" c10 F)) 
    (bcStrings (8 "2.1" c11 F)) 
    (bcStrings (8 "2.2" c12 F)) 
    (bcStrings (8 "1.5833" c13 F))      
    (bcStrings (8 "1.7167" c14 F))      
    (bcStrings (8 "1.95" c15 F))        
    (bcStrings (8 "2.2833" c16 F))      
    (bcStrings (8 "2.4833" c17 F))      
    (bcStrings (8 "2.5833" c18 F))      
    (bcStrings (8 "2.1433" c19 F))      
    (bcStrings (8 "2.2767" c20 F))      
    (bcStrings (8 "2.51" c21 F))        
    (bcStrings (8 "2.8433" c22 F))      
    (bcStrings (8 "3.0433" c23 F))      
    (bcStrings (8 "3.1433" c24 F))      
    (bcStrings (8 "2.8667" c25 F))      
    (bcStrings (8 "3" c26 F))   
    (bcStrings (8 "3.2333" c27 F))      
    (bcStrings (8 "3.5667" c28 F))      
    (bcStrings (8 "3.7667" c29 F))      
    (bcStrings (8 "3.8667" c30 F))      
    (bcStrings (8 "3.4667" c31 F))      
    (bcStrings (8 "3.6" c32 F)) 
    (bcStrings (8 "3.8333" c33 F))      
    (bcStrings (8 "4.1667" c34 F))      
    (bcStrings (8 "4.3667" c35 F))      
    (bcStrings (8 "4.4667" c36 F))      
    (bcStrings (8 "4" c37 F))   
    (bcStrings (8 "4.1333" c38 F))      
    (bcStrings (8 "4.3667" c39 F))      
    (bcStrings (8 "4.7" c40 F)) 
    (bcStrings (8 "4.9" c41 F)) 
    (bcStrings (8 "5" c42 F))   
    (text . "\blankline"))
  htMakeDoneButton('"Continue",'e02dffGen)      
  htpSetProperty(page,'mx,mx)
  htpSetProperty(page,'my,my)
  htpSetProperty(page,'px,px)
  htpSetProperty(page,'py,py)
  htpSetProperty(page,'lwrk,lwrk)
  htpSetProperty(page,'liwrk,liwrk)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()

e02dffGen htPage ==
  mx := htpProperty(htPage,'mx)
  my := htpProperty(htPage,'my)
  px := htpProperty(htPage,'px)
  py := htpProperty(htPage,'py)
  lwrk := htpProperty(htPage,'lwrk)
  liwrk := htpProperty(htPage,'liwrk)
  ifail := htpProperty(htPage,'ifail)
  alist := htpInputAreaAlist htPage
  y := alist
  -- c
  for i in 1..((px-4)*(py-4)) repeat
    right := STRCONC ((first y).1," ")
    y := rest y
    cList := [right,:cList]
  cstring := bcwords2liststring cList
  -- mu 
  for i in 1..py repeat
    right := STRCONC ((first y).1," ")
    y := rest y    
    muList := [right,:muList]
  mustring := bcwords2liststring muList
  -- lamda
  for i in 1..px repeat
    right := STRCONC ((first y).1," ")
    y := rest y    
    lamList := [right,:lamList]
  lamstring := bcwords2liststring lamList
  -- y
  for i in 1..my repeat
    right := STRCONC ((first y).1," ")
    y := rest y    
    yList := [right,:yList]
  ystring := bcwords2liststring yList
  -- x 
  for i in 1..mx repeat
    right := STRCONC ((first y).1," ")
    y := rest y    
    xList := [right,:xList]
  xstring := bcwords2liststring xList
  prefix := STRCONC('"e02dff(",STRINGIMAGE mx,", ",STRINGIMAGE my,", ")
  prefix := STRCONC(prefix,STRINGIMAGE px,", ",STRINGIMAGE py,",[")
  prefix := STRCONC(prefix,xstring,"],[",ystring,"],[",lamstring,"],[")
  prefix := STRCONC(prefix,mustring,"],[",cstring,"],",STRINGIMAGE lwrk,", ")
  prefix := STRCONC(prefix,STRINGIMAGE liwrk,", ",STRINGIMAGE ifail,")")
  linkGen prefix

e02gaf() ==
  htInitPage('"E02GAF - \htbitmap{l1}-approximation by general linear function",nil)
  htMakePage '(
    (domainConditions 
      (isDomain F (Float))
       (isDomain PI (PositiveInteger)))
    (text . "\windowlink{Manual Page}{manpageXXe02gaf} for this routine ")
    (text . "\newline ")
    (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|e02gaf| '|NagFittingPackage|)} for this routine")
    (text . "\newline \horizontalline ")
    (text . "Calculates an \htbitmap{l1} solution to the over determined system")
    (text . " of linear equations {\it Ax = b}, where A is an {\it m} by {\it n")
    (text . "} matrix, {\it x} is an {\it n} element vector, and {\it b} is an ")
    (text . "{\it m} element vector.  The matrix {\it A} need not be of full ")
    (text . "rank. \blankline ")
    (text . "\menuitemstyle{}\tab{2} \newline ")
    (text . "Number of rows of {\it A}, {\it m}: ")
    (text . "\newline\tab{2} ")
    (bcStrings (6 5 m PI))
    (text . "\blankline ")
    (text . "\menuitemstyle{}\tab{2} \newline ")
    (text . "Number of columns of {\it A}, {\it n}: ")
    (text . "\newline\tab{2} ")
    (bcStrings (6 3 n PI))
--    (text . "\blankline ")
--    (text . "\menuitemstyle{}\tab{2} \newline ")
--    (text . "First dimension of {\it A(la,n+2)}, {\it la}\htbitmap{great=}")
--    (text . " {\it m + 2}: \newline\tab{2} ")
--    (bcStrings (6 7 la PI))
    (text . "\blankline ")
    (text . "\menuitemstyle{}\tab{2} \newline ")
    (text . "Tolerance (default is zero), {\it toler}: ")
    (text . "\newline\tab{2} ")
    (bcStrings (10 "0.0" toler F))
    (text . "\blankline ")
    (text . "\menuitemstyle{} \tab{2} ")
    (text . "\newline \tab{2} ")
    (text . "Ifail value: ")
    (radioButtons ifail
        ("" " -1, Print error messages" minusOne)
        ("" "  1, Suppress error messages" one)))
  htMakeDoneButton('"Continue", 'e02gafSolve)
  htShowPage()

e02gafSolve htPage ==
  m :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'm)
    objValUnwrap htpLabelSpadValue(htPage, 'm)
  n :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'n)
    objValUnwrap htpLabelSpadValue(htPage, 'n)
  la := m+2
--    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'la)
--    objValUnwrap htpLabelSpadValue(htPage, 'la)
  toler := htpLabelInputString(htPage,'toler)
  error := htpButtonValue(htPage,'ifail)
  ifail :=
    error = 'one => '1
    '-1
  ((m = 5 and n = 3) and la = 7) => e02gafDefaultSolve (htPage,toler,ifail) 
  labelList :=
    "append"/[fc(i,n) for i in 1..la] where fc(i,n) ==
      tempList := 
        "append"/[fr(i,j) for j in 1..(n+2)] where fr(i,j) == 
          fnam := INTERN STRCONC ('"f",STRINGIMAGE i, STRINGIMAGE j)
          [['bcStrings,[9, 0.0, fnam, 'F]]]
      prefix := ('"\newline ")
      tempList := [['text,:prefix],:tempList]
  bList := 
    "append"/[fb(i) for i in 1..m] where fb(i) ==
      lnam := INTERN STRCONC ('"l",STRINGIMAGE i)
      [['bcStrings,[9, 0.0, lnam, 'F]]]
  prefix := ('"\blankline \menuitemstyle{} \tab{2} Values of {\it B(m)}: \newline")
  bList := [['text,:prefix],:bList]
  equationPart := [
     '(domainConditions 
        (isDomain P (Polynomial $EmptyMode))
         (isDomain S (String))
          (isDomain F (Float))
           (isDomain PI (PositiveInteger))),
            :labelList,:bList]
  page := htInitPage('"E02GAF - \htbitmap{l1}-approximation by general linear function",nil)
  htSay '"\menuitemstyle{}\tab{2} Values of {\it A(la,n+2)}:"
  htMakePage equationPart
  htSay '"\blankline "
  htMakeDoneButton('"Continue",'e02gafGen)
  htpSetProperty(page,'m,m)
  htpSetProperty(page,'n,n)
  htpSetProperty(page,'la,la)
  htpSetProperty(page,'toler,toler)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()


e02gafDefaultSolve (htPage,toler,ifail) ==
  m := '5
  n := '3
  la := '7
  page := htInitPage('"E02GAF - \htbitmap{l1}-approximation by general linear function",nil)
  htMakePage '(
    (domainConditions 
       (isDomain F (Float)))
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2} Values of {\it A(la,n+2)}:")
    (text . "\newline ")
    (bcStrings (9 "1.0" a11 F))  
    (bcStrings (9 "1.0" a12 F))  
    (bcStrings (9 "1.0" a13 F))  
    (bcStrings (9 "0.0" a14 F))  
    (bcStrings (9 "0.0" a15 F))  
    (text . "\newline ")
    (bcStrings (9 "exp(0.2)" a21 F))  
    (bcStrings (9 "exp(-0.2)" a22 F))  
    (bcStrings (9 "1.0" a23 F))  
    (bcStrings (9 "0.0" a24 F))  
    (bcStrings (9 "0.0" a25 F))  
    (text . "\newline ")
    (bcStrings (9 "exp(0.4)" a31 F))  
    (bcStrings (9 "exp(-0.4)" a32 F))  
    (bcStrings (9 "1.0" a33 F))  
    (bcStrings (9 "0.0" a34 F))  
    (bcStrings (9 "0.0" a35 F))  
    (text . "\newline ")
    (bcStrings (9 "exp(0.6)" a41 F))  
    (bcStrings (9 "exp(-0.6)" a42 F))  
    (bcStrings (9 "1.0" a43 F))  
    (bcStrings (9 "0.0" a44 F))  
    (bcStrings (9 "0.0" a45 F))  
    (text . "\newline ")
    (bcStrings (9 "exp(0.8)" a51 F))  
    (bcStrings (9 "exp(-0.8)" a52 F))  
    (bcStrings (9 "1.0" a53 F))  
    (bcStrings (9 "0.0" a54 F))  
    (bcStrings (9 "0.0" a55 F))  
    (text . "\newline ")
    (bcStrings (9 "0.0" a61 F))  
    (bcStrings (9 "0.0" a62 F))  
    (bcStrings (9 "0.0" a63 F))  
    (bcStrings (9 "0.0" a64 F))  
    (bcStrings (9 "0.0" a65 F))  
    (text . "\newline ")
    (bcStrings (9 "0.0" a71 F))  
    (bcStrings (9 "0.0" a72 F))  
    (bcStrings (9 "0.0" a73 F))  
    (bcStrings (9 "0.0" a74 F))  
    (bcStrings (9 "0.0" a75 F))  
    (text . "\blankline ")
    (text . "\menuitemstyle{}\tab{2} Values of {\it B(m)}:")
    (text . "\newline ")
    (bcStrings (9 "4.501" b1 F))  
    (bcStrings (9 "4.36" b2 F))  
    (bcStrings (9 "4.333" b3 F))  
    (bcStrings (9 "4.418" b4 F))  
    (bcStrings (9 "4.625" b5 F)))  
  htMakeDoneButton('"Continue",'e02gafGen)      
  htpSetProperty(page,'m,m)
  htpSetProperty(page,'n,n)
  htpSetProperty(page,'toler,toler)
  htpSetProperty(page,'la,la)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()

e02gafGen htPage ==
  m := htpProperty(htPage,'m)
  n := htpProperty(htPage,'n)
  la := htpProperty(htPage,'la)
  toler := htpProperty(htPage,'toler)
  ifail := htpProperty(htPage,'ifail)
  nplustwo := n + 2
  alist := htpInputAreaAlist htPage
  y := alist
  for i in 1..m repeat
    right := STRCONC ((first y).1," ")
    y := rest y    
    blist := [right,:blist]    
  bstring := bcwords2liststring blist
  y := REVERSE y
  k := -1
  matform := [[y.(k := k + 1).1 for j in 0..(nplustwo-1)] for i in 0..(la-1)]
  matstring := bcwords2liststring [bcwords2liststring x for x in matform] 
  prefix := STRCONC('"e02gaf(",STRINGIMAGE m,", ",STRINGIMAGE la,", ")
  prefix := STRCONC(prefix,STRINGIMAGE nplustwo,", ",STRINGIMAGE toler,", ")
  prefix := STRCONC(prefix,matstring,",[",bstring,"], ")
  prefix := STRCONC(prefix,STRINGIMAGE ifail,")")
  linkGen prefix


e02daf() == 
  htInitPage('"E02DAF - Least-squares surface fit, bicubic splines",nil)
  htMakePage '(
    (domainConditions 
       (isDomain PI (PositiveInteger)))
    (text . "\windowlink{Manual Page}{manpageXXe02daf} for this routine ")
    (text . "\newline ")
    (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|e02daf| '|NagFittingPackage|)} for this routine")
    (text . "\newline \horizontalline ")
    (text . "Determines a minimal, least squares bicubic B-spline surface fit")
    (text . "\htbitmap{e02daf} to the set of points ")
    (text . "{\em (\htbitmap{xr},\htbitmap{yr},\htbitmap{fr})} with weights ")
    (text . "\htbitmap{wr}, for r = 1,2,...,m. The user must supply internal ")
    (text . "knot sets {\lambda},in the x-direction and {\mu} in the ")
    (text . "y-direction, ")
    (text . "which can be thought of as dividing the data region into panels;")
    (text . "s(x,y) consists of a seperate bicubic polynomial in each panel, ")
    (text . "the polynomial joining together with second derivative ")
    (text . "continuity. Eight additional (external) knots are added to each ")
    (text . "of the knot sets by this routine. The routine minimizes \Sigma, ")
    (text . "the sum of squares of the weighted residuals ")
    (text . "\htbitmap{e02daf1}, for r = 1,2,...,m, subject to the ")
    (text . "given knot sets. \newline ")
    (text . "A call of this routine should be preceded by a call of E02ZAF ")
    (text . "to provide indexing information. ")
    (text . "\blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} \newline ")
    (text . "Number of data points {\it m}: ")
    (text . "\newline\tab{2} ")
    (bcStrings (6 30 m PI))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2}")
    (text . "\newline Knots in x direction {\em px}")
    (text . "\htbitmap{great=} 8: \tab{32} \menuitemstyle{}\tab{34}")
    (text . "Knots in y direction {\em py}:")
    (text . "\newline\tab{2} ")
    (bcStrings (6 8 px PI))
    (text . "\tab{34} ")
    (bcStrings (6 10 py PI))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2}")
    (text . "\newline Rank threshold {\em eps}:")
    (text . "\newline\tab{2} ")
    (bcStrings (10 "0.000001" eps F))
    (text . "\blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} \newline ")
    (text . "Dimension of point {\it npoint}: ")
    (text . "\newline\tab{2} ")
    (bcStrings (6 43 npoint PI))
    -- include a radio button later to allow switching of
    -- x & y if px <= py 
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{} \tab{2} ")
    (text . "\newline \tab{2} ")
    (text . "Ifail value: ")
    (radioButtons ifail
        ("" " -1, Print error messages" minusOne)
        ("" "  1, Suppress error messages" one)))
  htMakeDoneButton('"Continue", 'e02dafSolve)
  htShowPage()

e02dafSolve htPage ==
  m :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'm)
    objValUnwrap htpLabelSpadValue(htPage, 'm)
  px :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'px)
    objValUnwrap htpLabelSpadValue(htPage, 'px)
  py :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'py)
    objValUnwrap htpLabelSpadValue(htPage, 'py)
  npoint :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'npoint)
    objValUnwrap htpLabelSpadValue(htPage, 'npoint)
  nc := (px - 4)*(py - 4)
  nws := (2*nc + 1)*(3*py - 6) -2
  eps := htpLabelInputString(htPage,'eps)
  error := htpButtonValue(htPage,'ifail)
  ifail :=
    error = 'one => '1
    '-1
  ((m = '30 and px = '8) and py = '10) => e02dafDefaultSolve(htPage,eps,nws,npoint,ifail)
  labelList :=
    "append"/[fxy(i) for i in 1..m] where fxy(i) ==
      prefix := ('"\newline \tab{2} ")
      middle := ('"\tab{17} ")
      next := ('"\tab{32} ")
      end := ('"\tab{47} ")
      xnam := INTERN STRCONC ('"x",STRINGIMAGE i)
      ynam := INTERN STRCONC ('"y",STRINGIMAGE i)
      fnam := INTERN STRCONC ('"f",STRINGIMAGE i)
      wnam := INTERN STRCONC ('"w",STRINGIMAGE i)
      [['text,:prefix],['bcStrings,[8, 0.0, xnam, 'F]], 
       ['text,:middle],['bcStrings,[8, 0.0, ynam, 'F]],
         ['text,:next],['bcStrings,[8, 0.0, fnam, 'F]],
           ['text,:end],['bcStrings,[8, 0.0, wnam, 'F]]]
  lamList :=
    "append"/[flam(i) for i in 5..(px-4)] where flam(i) ==
      lnam := INTERN STRCONC ('"l",STRINGIMAGE i)
      [['bcStrings,[8, 0.0, lnam, 'F]]]
  prefix := ('"\blankline \menuitemstyle{} \tab{2} lamda(5) to lamda(px-4): ")
  prefix := STRCONC(prefix,"\newline \tab{2} ")
  postfix := ('"\newline \blankline ")
  lamList := [['text,:prefix],:lamList,['text,:postfix]]
  muList :=
    "append"/[fmu(i) for i in 5..(py-4)] where fmu(i) ==
      mnam := INTERN STRCONC ('"m",STRINGIMAGE i)
      [['bcStrings,[8, 0.0, mnam, 'F]]]
  prefix := ('"\menuitemstyle{} \tab{2} mu(5) to mu(py-4):")
  prefix := STRCONC(prefix,"\newline \tab{2} ")
  muList := [['text,:prefix],:muList]
  pList :=
    "append"/[fp(i) for i in 1..npoint] where fp(i) ==
      prefix := ('"\newline \tab{2} ")
      pnam := INTERN STRCONC ('"p",STRINGIMAGE i)
      [['text,:prefix],['bcStrings,[8, 0.0, pnam, 'F]]]
  prefix := ('"\blankline \menuitemstyle{} \tab{2} Enter values of Point: ")
  pList := [['text,:prefix],:pList]
  equationPart := [
     '(domainConditions 
        (isDomain P (Polynomial $EmptyMode))
         (isDomain S (String))
          (isDomain F (Float))
           (isDomain PI (PositiveInteger))),
            :labelList,:lamList,:muList,:pList]
  page := htInitPage('"E02DAF - Least-squares surface fit, bicubic splines",nil)
  htSay '"\menuitemstyle{}\tab{2} Values of \htbitmap{xr}: "
  htSay '"\tab{15} \menuitemstyle{}\tab{17} Values of \htbitmap{yr}:"
  htSay '"\tab{30} \menuitemstyle{}\tab{32} Values of \htbitmap{fr}:" 
  htSay '"\tab{44} \menuitemstyle{}\tab{46} Values of \htbitmap{wr}:"
  htMakePage equationPart
  htSay '"\blankline "
  htMakeDoneButton('"Continue",'e02dafGen)
  htpSetProperty(page,'m,m)
  htpSetProperty(page,'px,px)
  htpSetProperty(page,'py,py)
  htpSetProperty(page,'nws,nws)
  htpSetProperty(page,'eps,eps)
  htpSetProperty(page,'npoint,npoint)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()


e02dafDefaultSolve (htPage,eps,nws,npoint,ifail) ==
  m := '30
  px := '8
  py := '10
  page := htInitPage('"E02DAF - Least-squares surface fit, bicubic splines",nil)
  htMakePage '(
    (domainConditions 
       (isDomain F (Float)))
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2} Values of \htbitmap{xr}:")
    (text . "\tab{15} \menuitemstyle{} \tab{17} Values of ")
    (text . "\htbitmap{yr}: \tab{30} \menuitemstyle{} \tab{32} ")
    (text . "Values of \htbitmap{fr}: \tab{44} \menuitemstyle{} ")
    (text . "\tab{46} Values of \htbitmap{wr}:")
    (text . "\newline \tab{2} ")
    (bcStrings (8 "-0.52" x1 F))  
    (text . "\tab{17}")
    (bcStrings (8 "0.60" y1 F))    
    (text . "\tab{32}")
    (bcStrings (8 "0.93" f1 F))    
    (text . "\tab{47}")
    (bcStrings (8 "10" w1 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "-0.61" x2 F))  
    (text . "\tab{17}")
    (bcStrings (8 "-0.95" y2 F))    
    (text . "\tab{32}")
    (bcStrings (8 "-1.79" f2 F))    
    (text . "\tab{47}")
    (bcStrings (8 "10" w2 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "0.93" x3 F))  
    (text . "\tab{17}")
    (bcStrings (8 "0.87" y3 F))    
    (text . "\tab{32}")
    (bcStrings (8 "0.36" f3 F))    
    (text . "\tab{47}")
    (bcStrings (8 "10" w3 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "0.09" x4 F))  
    (text . "\tab{17}")
    (bcStrings (8 "0.84" y4 F))    
    (text . "\tab{32}")
    (bcStrings (8 "0.52" f4 F))    
    (text . "\tab{47}")
    (bcStrings (8 "10" w4 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "0.88" x5 F))  
    (text . "\tab{17}")
    (bcStrings (8 "0.17" y5 F))    
    (text . "\tab{32}")
    (bcStrings (8 "0.49" f5 F))    
    (text . "\tab{47}")
    (bcStrings (8 "10" w5 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "-0.70" x6 F))  
    (text . "\tab{17}")
    (bcStrings (8 "-0.87" y6 F))    
    (text . "\tab{32}")
    (bcStrings (8 "-1.76" f6 F))    
    (text . "\tab{47}")
    (bcStrings (8 "10" w6 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "1" x7 F))  
    (text . "\tab{17}")
    (bcStrings (8 "1" y7 F))    
    (text . "\tab{32}")
    (bcStrings (8 "0.33" f7 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w7 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "1" x8 F))  
    (text . "\tab{17}")
    (bcStrings (8 "0.1" y8 F))    
    (text . "\tab{32}")
    (bcStrings (8 "0.48" f8 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w8 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "0.3" x9 F))  
    (text . "\tab{17}")
    (bcStrings (8 "0.24" y9 F))    
    (text . "\tab{32}")
    (bcStrings (8 "0.65" f9 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w9 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "-0.77" x10 F))  
    (text . "\tab{17}")
    (bcStrings (8 "-0.77" y10 F))    
    (text . "\tab{32}")
    (bcStrings (8 "-1.82" f10 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w10 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "-0.23" x11 F))  
    (text . "\tab{17}")
    (bcStrings (8 "0.32" y11 F))    
    (text . "\tab{32}")
    (bcStrings (8 "0.92" f11 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w11 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "-1" x12 F))  
    (text . "\tab{17}")
    (bcStrings (8 "1" y12 F))    
    (text . "\tab{32}")
    (bcStrings (8 "1" f12 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w12 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "-0.26" x13 F))  
    (text . "\tab{17}")
    (bcStrings (8 "-0.63" y13 F))    
    (text . "\tab{32}")
    (bcStrings (8 "8.88" f13 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w13 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "-0.83" x14 F))  
    (text . "\tab{17}")
    (bcStrings (8 "-0.66" y14 F))    
    (text . "\tab{32}")
    (bcStrings (8 "-2.01" f14 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w14 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "0.22" x15 F))  
    (text . "\tab{17}")
    (bcStrings (8 "0.93" y15 F))    
    (text . "\tab{32}")
    (bcStrings (8 "0.47" f15 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w15 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "0.89" x16 F))  
    (text . "\tab{17}")
    (bcStrings (8 "0.15" y16 F))    
    (text . "\tab{32}")
    (bcStrings (8 "0.49" f16 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w16 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "-0.80" x17 F))  
    (text . "\tab{17}")
    (bcStrings (8 "0.99" y17 F))    
    (text . "\tab{32}")
    (bcStrings (8 "0.84" f17 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w17 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "-0.88" x18 F))  
    (text . "\tab{17}")
    (bcStrings (8 "-0.54" y18 F))    
    (text . "\tab{32}")
    (bcStrings (8 "-2.42" f18 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w18 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "0.68" x19 F))  
    (text . "\tab{17}")
    (bcStrings (8 "0.44" y19 F))    
    (text . "\tab{32}")
    (bcStrings (8 "0.47" f19 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w19 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "-0.14" x20 F))  
    (text . "\tab{17}")
    (bcStrings (8 "-0.72" y20 F))    
    (text . "\tab{32}")
    (bcStrings (8 "7.15" f20 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w20 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "0.67" x21 F))  
    (text . "\tab{17}")
    (bcStrings (8 "0.63" y21 F))    
    (text . "\tab{32}")
    (bcStrings (8 "0.44" f21 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w21 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "-0.90" x22 F))  
    (text . "\tab{17}")
    (bcStrings (8 "-0.40" y22 F))    
    (text . "\tab{32}")
    (bcStrings (8 "-3.34" f22 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w22 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "-0.84" x23 F))  
    (text . "\tab{17}")
    (bcStrings (8 "0.20" y23 F))    
    (text . "\tab{32}")
    (bcStrings (8 "2.78" f23 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w23 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "0.84" x24 F))  
    (text . "\tab{17}")
    (bcStrings (8 "0.43" y24 F))    
    (text . "\tab{32}")
    (bcStrings (8 "0.44" f24 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w24 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "0.15" x25 F))  
    (text . "\tab{17}")
    (bcStrings (8 "0.28" y25 F))    
    (text . "\tab{32}")
    (bcStrings (8 "0.70" f25 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w25 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "-0.91" x26 F))  
    (text . "\tab{17}")
    (bcStrings (8 "-0.24" y26 F))    
    (text . "\tab{32}")
    (bcStrings (8 "-6.52" f26 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w26 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "-0.35" x27 F))  
    (text . "\tab{17}")
    (bcStrings (8 "0.86" y27 F))    
    (text . "\tab{32}")
    (bcStrings (8 "0.66" f27 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w27 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "-0.16" x28 F))  
    (text . "\tab{17}")
    (bcStrings (8 "-0.41" y28 F))    
    (text . "\tab{32}")
    (bcStrings (8 "2.32" f28 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w28 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "-0.35" x29 F))  
    (text . "\tab{17}")
    (bcStrings (8 "-0.05" y29 F))    
    (text . "\tab{32}")
    (bcStrings (8 "1.66" f29 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w29 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "-1" x30 F))  
    (text . "\tab{17}")
    (bcStrings (8 "-1" y30 F))    
    (text . "\tab{32}")
    (bcStrings (8 "-1" f30 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w30 F))    
    (text . "\blankline ")
    (text . "\menuitemstyle{}\tab{2} mu(5) to mu(py-4):")
    (text . "\newline \tab{2}")
    (bcStrings (8 "-0.50" mu5 F))
    (bcStrings (8 "0.00" mu6 F))
    (text . "\blankline ")
    (text . "\menuitemstyle{}\tab{2} Enter values for point:")
    (text . "\newline \tab{2}")
    (bcStrings (6 3 p1 PI))     
    (text . "\newline \tab{2}")
    (bcStrings (6 6 p2 PI))     
    (text . "\newline \tab{2}")
    (bcStrings (6 4 p3 PI))     
    (text . "\newline \tab{2}")
    (bcStrings (6 5 p4 PI))     
    (text . "\newline \tab{2}")
    (bcStrings (6 7 p5 PI))     
    (text . "\newline \tab{2}")
    (bcStrings (6 10 p6 PI))    
    (text . "\newline \tab{2}")
    (bcStrings (6 8 p7 PI))     
    (text . "\newline \tab{2}")
    (bcStrings (6 9 p8 PI))     
    (text . "\newline \tab{2}")
    (bcStrings (6 11 p9 PI))    
    (text . "\newline \tab{2}")
    (bcStrings (6 13 p10 PI))   
    (text . "\newline \tab{2}")
    (bcStrings (6 12 p11 PI))   
    (text . "\newline \tab{2}")
    (bcStrings (6 15 p12 PI))   
    (text . "\newline \tab{2}")
    (bcStrings (6 14 p13 PI))   
    (text . "\newline \tab{2}")
    (bcStrings (6 18 p14 PI))   
    (text . "\newline \tab{2}")
    (bcStrings (6 16 p15 PI))   
    (text . "\newline \tab{2}")
    (bcStrings (6 17 p16 PI))   
    (text . "\newline \tab{2}")
    (bcStrings (6 19 p17 PI))   
    (text . "\newline \tab{2}")
    (bcStrings (6 20 p18 PI))   
    (text . "\newline \tab{2}")
    (bcStrings (6 21 p19 PI))   
    (text . "\newline \tab{2}")
    (bcStrings (6 30 p20 PI))   
    (text . "\newline \tab{2}")
    (bcStrings (6 23 p21 PI))   
    (text . "\newline \tab{2}")
    (bcStrings (6 26 p22 PI))   
    (text . "\newline \tab{2}")
    (bcStrings (6 24 p23 PI))   
    (text . "\newline \tab{2}")
    (bcStrings (6 25 p24 PI))   
    (text . "\newline \tab{2}")
    (bcStrings (6 27 p25 PI))   
    (text . "\newline \tab{2}")
    (bcStrings (6 28 p26 PI))   
    (text . "\newline \tab{2}")
    (bcStrings (6 0 p27 PI))    
    (text . "\newline \tab{2}")
    (bcStrings (6 29 p28 PI))   
    (text . "\newline \tab{2}")
    (bcStrings (6 0 p29 PI))    
    (text . "\newline \tab{2}")
    (bcStrings (6 0 p30 PI))    
    (text . "\newline \tab{2}")
    (bcStrings (6 2 p31 PI))    
    (text . "\newline \tab{2}")
    (bcStrings (6 22 p32 PI))   
    (text . "\newline \tab{2}")
    (bcStrings (6 1 p33 PI))    
    (text . "\newline \tab{2}")
    (bcStrings (6 0 p34 PI))    
    (text . "\newline \tab{2}")
    (bcStrings (6 0 p35 PI))    
    (text . "\newline \tab{2}")
    (bcStrings (6 0 p36 PI))    
    (text . "\newline \tab{2}")
    (bcStrings (6 0 p37 PI))    
    (text . "\newline \tab{2}")
    (bcStrings (6 0 p38 PI))    
    (text . "\newline \tab{2}")
    (bcStrings (6 0 p39 PI))    
    (text . "\newline \tab{2}")
    (bcStrings (6 0 p40 PI))    
    (text . "\newline \tab{2}")
    (bcStrings (6 0 p41 PI))    
    (text . "\newline \tab{2}")
    (bcStrings (6 0 p42 PI))    
    (text . "\newline \tab{2}")
    (bcStrings (6 0 p43 PI))    
    (text . "\blankline"))
  htMakeDoneButton('"Continue",'e02dafGen)      
  htpSetProperty(page,'m,m)
  htpSetProperty(page,'px,px)
  htpSetProperty(page,'py,py)
  htpSetProperty(page,'nws,nws)
  htpSetProperty(page,'eps,eps)
  htpSetProperty(page,'npoint,npoint)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()

e02dafGen htPage ==
  m := htpProperty(htPage,'m)
  px := htpProperty(htPage,'px)
  py := htpProperty(htPage,'py)
  nws := htpProperty(htPage,'nws)
  eps := htpProperty(htPage,'eps)
  npoint := htpProperty(htPage,'npoint)
  ifail := htpProperty(htPage,'ifail)
  alist := htpInputAreaAlist htPage
  y := alist
  nadres := (px-7)*(py-7)
  -- point
  for i in 1..npoint repeat
    right := STRCONC ((first y).1," ")
    y := rest y
    pointList := [right,:pointList]
  pstring := bcwords2liststring pointList
  -- mu 
  for i in 1..4 repeat
    muList := ['"0 ",:muList]
  for i in 5..(py-4) repeat
    right := STRCONC ((first y).1," ")
    y := rest y    
    muList := [right,:muList]
  for i in (py-3)..py repeat
    muList := ['"0 ",:muList]
  mustring := bcwords2liststring muList
  -- lamda
  for i in 1..4 repeat
    lamList := ['"0 ",:lamList]
  for i in 5..(px-4) repeat
    right := STRCONC ((first y).1," ")
    y := rest y    
    lamList := [right,:lamList]
  for i in (px-3)..px repeat
    lamList := ['"0 ",:lamList]
  lamstring := bcwords2liststring lamList
  -- x & y
  while y repeat
    one := STRCONC((first y).1," ")
    y := rest y
    two := STRCONC((first y).1," ")
    y := rest y
    three := STRCONC ((first y).1," ")
    y := rest y    
    four :=  STRCONC ((first y).1," ")
    y := rest y
    xlist := [four,:xlist]
    ylist := [three,:ylist]
    flist := [two,:flist]
    wlist := [one,:wlist]
  xstring := bcwords2liststring xlist
  ystring := bcwords2liststring ylist
  fstring := bcwords2liststring flist
  wstring := bcwords2liststring wlist
  nc := (px-4)*(py-4)
  prefix := STRCONC('"e02daf(",STRINGIMAGE m,", ",STRINGIMAGE px,", ")
  prefix := STRCONC(prefix,STRINGIMAGE py,",[",xstring,"],[",ystring,"],[")
  prefix := STRCONC(prefix,fstring,"],[",wstring,"],[",mustring,"],[")
  prefix := STRCONC(prefix,pstring,"], ",STRINGIMAGE npoint,", ")
  prefix := STRCONC(prefix,STRINGIMAGE nc,", ",STRINGIMAGE nws,", ",eps,", [")
  prefix := STRCONC(prefix,lamstring,"], ",STRINGIMAGE ifail,")")
  linkGen prefix


e02dcf() ==
  htInitPage('"E02DCF - Least-squares curve by bicubic splines with automatic knot placement, data on a rectangular grid",nil)
  htMakePage '(
    (domainConditions 
       (isDomain PI (PositiveInteger)))
    (text . "\windowlink{Manual Page}{manpageXXe02dcf} for this routine ")
    (text . "\newline ")
    (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|e02dcf| '|NagFittingPackage|)} for this routine")
    (text . "\newline \horizontalline ")
    (text . "\newline ")
    (text . "Determines a bicubic spline approximation to a set of points ")
    (text . "given on a rectangular grid defined by \htbitmap{mx} ")
    (text . "points \htbitmap{xq} on the x-axis and ")
    (text . "\htbitmap{my} points \htbitmap{yr} on the ")
    (text . "y-axix. The knots \htbitmap{lamdai}, for i = 1,2,...,")
    (text . "\htbitmap{nx} and \htbitmap{mui}, for ")
    (text . "i = 1,2,...,\htbitmap{ny} are chosen for this routine ")
    (text . ", but a single parameter S must be specified to control the ")
    (text . "trade-off between closeness of fit and smoothness of fit.  This ")
    (text . "affects the number of knots required by the spline, which is ")
    (text . "given in the B-spline representation \htbitmap{e02daf}")
    (text . "\blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} \newline ")
    (text . "Grid points on x-axis \htbitmap{mx}: ")    
    (text . "\tab{30} \menuitemstyle{}\tab{32} Grid points on y-axis ")
    (text . "\htbitmap{my}: ")    
    (text . "\newline\tab{2} ")
    (bcStrings (6 11 mx PI))
    (text . "\tab{32} ")
    (bcStrings (6 9 my PI))    
    (text . "\blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} ")
    (text . "Over-estimate of \htbitmap{nx} of knots \lambda of the ") 
    (text . "computed spline, {\it nxest}: \newline\tab{2} ")
    (bcStrings (6 15 nxest PI))
    (text . "\blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} ")
    (text . "Over-estimate of \htbitmap{ny} of knots \mu of the computed ") 
    (text . "spline, {\it nyest}: \newline\tab{2} ")
    (bcStrings (6 13 nyest PI))
    (text . "\blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} \newline ")
    (text . "Smoothing factor {\it s}: ")
    (text . "\newline\tab{2} ")
    (bcStrings (6 "0.1" s F))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{} \tab{2} ")
    (text . "\newline \tab{2} ")
    (text . "Start value: ")
    (radioButtons start
        ("" "  Cold Start - no values needed for {\it nx,ny,lamda,mu} or {\it iwrk}" cold)
        ("" "  Warm Start - uses knots found in a previous call" warm))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{} \tab{2} ")
    (text . "\newline \tab{2} ")
    (text . "Ifail value: ")
    (radioButtons ifail
        ("" " -1, Print error messages" minusOne)
        ("" "  1, Suppress error messages" one)))
  htMakeDoneButton('"Continue", 'e02dcfSolve)
  htShowPage()

e02dcfSolve htPage ==
  mx :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'mx)
    objValUnwrap htpLabelSpadValue(htPage, 'mx)
  my :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'my)
    objValUnwrap htpLabelSpadValue(htPage, 'my)
  nxest :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'nxest)
    objValUnwrap htpLabelSpadValue(htPage, 'nxest)
  nyest :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'nyest)
    objValUnwrap htpLabelSpadValue(htPage, 'nyest)
  wrklist := [my,nxest]
  wrkmax := APPLY ('MAX, wrklist)
  lwrk := 4*(mx + my) +11*(nxest + nyest) + nxest*my + wrkmax +54
  liwrk := 3 + mx + my + nxest + nyest
  s := htpLabelInputString(htPage,'s)
  initial := htpButtonValue(htPage,'start)
  start :=
    initial = 'cold => '1
    '2
  error := htpButtonValue(htPage,'ifail)
  ifail :=
    error = 'one => '1
    '-1
  ((mx = 11 and my = 9) and start = 1) => 
                e02dcfDefaultSolve (htPage,nxest,nyest,lwrk,liwrk,s,ifail)
  start = 1 => e02dcfColdSolve (htPage,mx,my,nxest,nyest,lwrk,liwrk,s,ifail)
  -- warm start not really possible from hyperdoc
  -- as inputing a workspace array of dimension 592 is asking too much
  -- user should use the command line, using the previous calculated 
  -- parameters
  htInitPage('"E02DCF - Least-squares curve by bicubic splines with automatic knot placement, data on a rectangular grid",nil)  
  htMakePage '(
    (domainConditions 
       (isDomain PI (PositiveInteger)))
    (text . "\blankline ")
    (text . "{\center{\em Hyperdoc interface not available for warm starts.}}")
    (text . "\newline ")
    (text . "{\center{\em Please use the command line.}}"))
  htMakeDoneButton('"Continue",'e02dcf)
  htShowPage()



e02dcfColdSolve(htPage,mx,my,nxest,nyest,lwrk,liwrk,s,ifail) ==
  xList :=
    "append"/[f(i) for i in 1..mx] where f(i) ==
      xnam := INTERN STRCONC ('"x",STRINGIMAGE i)
      [['bcStrings,[8, 0.0, xnam, 'F]]] 
  yList :=
    "append"/[g(i) for i in 1..my] where g(i) ==
      ynam := INTERN STRCONC ('"g",STRINGIMAGE i)
      [['bcStrings,[8, 0.0, ynam, 'F]]] 
  prefix:= ('"\blankline \menuitemstyle{}\tab{2} Values of {\it y(my)}: \newline ")
  yList := [['text,:prefix],:yList]
  fList :=
    "append"/[h(i) for i in 1..(mx*my)] where h(i) ==
      fnam := INTERN STRCONC ('"g",STRINGIMAGE i)
      [['bcStrings,[8, 0.0, fnam, 'F]]] 
  prefix:=('"\blankline \menuitemstyle{} \tab{2} Values of {\it f(mx*my)}: \newline ")
  fList := [['text,:prefix],:fList]
  equationPart := [
     '(domainConditions 
        (isDomain P (Polynomial $EmptyMode))
         (isDomain S (String))
          (isDomain F (Float))
           (isDomain PI (PositiveInteger))),
            :xList,:yList,:fList]
  page := htInitPage('"E02DCF - Least-squares curve by bicubic splines with automatic knot placement, data on a rectangular grid",nil)
  htSay '"\menuitemstyle{}\tab{2} Values of {\it x(mx)}: \newline "
  htMakePage equationPart
  htSay '"\blankline "
  htMakeDoneButton('"Continue",'e02dcfColdGen)
  htpSetProperty(page,'mx,mx)
  htpSetProperty(page,'my,my)
  htpSetProperty(page,'nxest,nxest)
  htpSetProperty(page,'nyest,nyest)
  htpSetProperty(page,'lwrk,lwrk)
  htpSetProperty(page,'liwrk,liwrk)
  htpSetProperty(page,'s,s)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()


e02dcfDefaultSolve (htPage,nxest,nyest,lwrk,liwrk,s,ifail) ==
  mx := 11
  my := 9
  page := htInitPage('"E02DCF - Least-squares curve by bicubic splines with automatic knot placement, data on a rectangular grid",nil)
  htMakePage '(
    (domainConditions 
       (isDomain F (Float)))
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2} Values of {\it x(mx)}:")
    (text . "\newline ")
    (bcStrings (8 "0" x1 F))
    (bcStrings (8 "0.5" x2 F))
    (bcStrings (8 "1" x3 F))
    (bcStrings (8 "1.5" x4 F))
    (bcStrings (8 "2" x5 F))
    (bcStrings (8 "2.5" x6 F))
    (bcStrings (8 "3" x7 F))
    (bcStrings (8 "3.5" x8 F))
    (bcStrings (8 "4" x9 F))
    (bcStrings (8 "4.5" x10 F))
    (bcStrings (8 "5" x11 F))
    (text . "\blankline ")
    (text . "\menuitemstyle{}\tab{2} Values of {\it y(my)}:")
    (text . "\newline ")
    (bcStrings (8 "0" y1 F))
    (bcStrings (8 "0.5" y2 F))
    (bcStrings (8 "1" y3 F))
    (bcStrings (8 "1.5" y4 F))
    (bcStrings (8 "2" y5 F))
    (bcStrings (8 "2.5" y6 F))
    (bcStrings (8 "3" y7 F))
    (bcStrings (8 "3.5" y8 F))
    (bcStrings (8 "4" y9 F))
    (text . "\blankline ")
    (text . "\menuitemstyle{}\tab{2} Values of {\it f(mx*my)}:")
    (text . "\newline ")
    (bcStrings (8 "1" f1 F))
    (bcStrings (8 "0.88758" f2 F))
    (bcStrings (8 "0.5403" f3 F))
    (bcStrings (8 "0.070737" f4 F))
    (bcStrings (8 "-0.41515" f5 F))
    (bcStrings (8 "-0.80114" f6 F))
    (bcStrings (8 "-0.97999" f7 F))
    (bcStrings (8 "-0.93446" f8 F))
    (bcStrings (8 "-0.65664" f9 F))
    (bcStrings (8 "1.5" f10 F))
    (bcStrings (8 "1.3564" f11 F))
    (bcStrings (8 "0.82045" f12 F))
    (bcStrings (8 "0.10611" f13 F))
    (bcStrings (8 "-0.62422" f14 F))
    (bcStrings (8 "-1.2317" f15 F))
    (bcStrings (8 "-1.485" f16 F))
    (bcStrings (8 "-1.3047" f17 F))
    (bcStrings (8 "-0.98547" f18 F))
    (bcStrings (8 "2.06" f19 F))
    (bcStrings (8 "1.7552" f20 F))
    (bcStrings (8 "1.0806" f21 F))
    (bcStrings (8 "0.15147" f22 F))
    (bcStrings (8 "-0.83229" f23 F))
    (bcStrings (8 "-1.6023" f24 F))
    (bcStrings (8 "-1.97" f25 F))
    (bcStrings (8 "-1.8729" f26 F))
    (bcStrings (8 "-1.4073" f27 F))
    (bcStrings (8 "2.57" f28 F))
    (bcStrings (8 "2.124" f29 F))
    (bcStrings (8 "1.3508" f30 F))
    (bcStrings (8 "0.17684" f31 F))
    (bcStrings (8 "-1.0404" f32 F))
    (bcStrings (8 "-2.0029" f33 F))
    (bcStrings (8 "-2.475" f34 F))
    (bcStrings (8 "-2.3511" f35 F))
    (bcStrings (8 "-1.6741" f36 F))
    (bcStrings (8 "3" f37 F))
    (bcStrings (8 "2.6427" f38 F))
    (bcStrings (8 "1.6309" f39 F))
    (bcStrings (8 "0.21221" f40 F))
    (bcStrings (8 "-1.2484" f41 F))
    (bcStrings (8 "-2.2034" f42 F))
    (bcStrings (8 "-2.97" f43 F))
    (bcStrings (8 "-2.8094" f44 F))
    (bcStrings (8 "-1.9809" f45 F))
    (bcStrings (8 "3.5" f46 F))
    (bcStrings (8 "3.1715" f47 F))
    (bcStrings (8 "1.8611" f48 F))
    (bcStrings (8 "0.24458" f49 F))
    (bcStrings (8 "-1.4565" f50 F))
    (bcStrings (8 "-2.864" f51 F))
    (bcStrings (8 "-3.265" f52 F))
    (bcStrings (8 "-3.2776" f53 F))
    (bcStrings (8 "-2.2878" f54 F))
    (bcStrings (8 "4.04" f55 F))
    (bcStrings (8 "3.5103" f56 F))
    (bcStrings (8 "2.0612" f57 F))
    (bcStrings (8 "0.28595" f58 F))
    (bcStrings (8 "-1.6946" f59 F))
    (bcStrings (8 "-3.2046" f60 F))
    (bcStrings (8 "-3.96" f61 F))
    (bcStrings (8 "-3.7958" f62 F))
    (bcStrings (8 "-2.6146" f63 F))
    (bcStrings (8 "4.5" f64 F))
    (bcStrings (8 "3.9391" f65 F))
    (bcStrings (8 "2.4314" f66 F))
    (bcStrings (8 "0.31632" f67 F))
    (bcStrings (8 "-1.8627" f68 F))
    (bcStrings (8 "-3.6351" f69 F))
    (bcStrings (8 "-4.455" f70 F))
    (bcStrings (8 "-4.2141" f71 F))
    (bcStrings (8 "-2.9314" f72 F))
    (bcStrings (8 "5.04" f73 F))
    (bcStrings (8 "4.3879" f74 F))
    (bcStrings (8 "2.7515" f75 F))
    (bcStrings (8 "0.35369" f76 F))
    (bcStrings (8 "-2.0707" f77 F))
    (bcStrings (8 "-4.0057" f78 F))
    (bcStrings (8 "-4.97" f79 F))
    (bcStrings (8 "-4.6823" f80 F))
    (bcStrings (8 "-3.2382" f81 F))
    (bcStrings (8 "5.505" f82 F))
    (bcStrings (8 "4.8367" f83 F))
    (bcStrings (8 "2.9717" f84 F))
    (bcStrings (8 "0.38505" f85 F))
    (bcStrings (8 "-2.2888" f86 F))
    (bcStrings (8 "-4.4033" f87 F))
    (bcStrings (8 "-5.445" f88 F))
    (bcStrings (8 "-5.1405" f89 F))
    (bcStrings (8 "-3.595" f90 F))
    (bcStrings (8 "6" f91 F))
    (bcStrings (8 "5.2755" f92 F))
    (bcStrings (8 "3.2418" f93 F))
    (bcStrings (8 "0.42442" f94 F))
    (bcStrings (8 "-2.4769" f95 F))
    (bcStrings (8 "-4.8169" f96 F))
    (bcStrings (8 "-5.93" f97 F))
    (bcStrings (8 "-5.6387" f98 F))
    (bcStrings (8 "-3.9319" f99 F)))
  htMakeDoneButton('"Continue",'e02dcfColdGen)          
  htpSetProperty(page,'mx,mx)
  htpSetProperty(page,'my,my)
  htpSetProperty(page,'nxest,nxest)
  htpSetProperty(page,'nyest,nyest)
  htpSetProperty(page,'lwrk,lwrk)
  htpSetProperty(page,'liwrk,liwrk)
  htpSetProperty(page,'s,s)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()

e02dcfColdGen htPage ==
  mx := htpProperty(htPage,'mx)
  my := htpProperty(htPage,'my)
  nxest := htpProperty(htPage,'nxest)
  nyest := htpProperty(htPage,'nyest)
  lwrk := htpProperty(htPage,'lwrk)
  liwrk := htpProperty(htPage,'liwrk)
  s := htpProperty(htPage,'s)
  cold := '"c"
  ifail := htpProperty(htPage,'ifail)
  alist := htpInputAreaAlist htPage
  y := alist
  for i in 1..(mx*my) repeat
    end := STRCONC((first y).1," ")
    y := rest y
    fList := [end,:fList]
  fstring := bcwords2liststring fList
  for i in 1..my repeat
    mid :=  STRCONC ((first y).1," ")
    y := rest y
    ylist := [mid,:ylist]
  ystring := bcwords2liststring ylist
  while y repeat
    start :=  STRCONC ((first y).1," ")
    y := rest y
    xlist := [start,:xlist]
  xstring := bcwords2liststring xlist
  -- additional entries needed to get it running
  -- but as Start = c they are not used
  prefix := STRCONC('"e02dcf(_"",cold,"_",",STRINGIMAGE mx,", [",xstring,"],")
  prefix := STRCONC(prefix,STRINGIMAGE my,",[",ystring,"],[",fstring,"], ")
  prefix := STRCONC(prefix,STRINGIMAGE s,", ",STRINGIMAGE nxest,", ")
  prefix := STRCONC(prefix,STRINGIMAGE nyest,", ",STRINGIMAGE lwrk,", ")
  prefix := STRCONC(prefix,STRINGIMAGE liwrk,",0,new(1,", STRINGIMAGE nxest,",0.0)$MATRIX DFLOAT,")
  prefix := STRCONC(prefix,"0,new(1,", STRINGIMAGE nyest,",0.0)$MATRIX DFLOAT,")
  end := STRCONC("new(1,", STRINGIMAGE lwrk,",0.0)$MATRIX DFLOAT,[[0 for i in 1..")
  end := STRCONC(end,STRINGIMAGE liwrk,"]]::Matrix Integer,",STRINGIMAGE ifail,")")
  linkGen STRCONC(prefix,end)


e02ddf() ==
  htInitPage('"E02DDF - Least-squares surface fit by bicubic splines with automatic knot placement, scattered data",nil)
  htMakePage '(
    (domainConditions 
       (isDomain PI (PositiveInteger)))
    (text . "\windowlink{Manual Page}{manpageXXe02ddf} for this routine ")
    (text . "\newline ")
    (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|e02ddf| '|NagFittingPackage|)} for this routine")
    (text . "\newline \horizontalline ")
    (text . "\newline ")
    (text . "Determines a bicubic spline approximation to a set of scattered")
    (text . " points ( \htbitmap{xr},\htbitmap{yr}, ")
    (text . "\htbitmap{fr})")
    (text . "with weights \htbitmap{wr}, for r = 1,2,...,m. ")
    (text . "The knots \htbitmap{lamdai}, for i = 1,2,...,")
    (text . "\htbitmap{nx} and \htbitmap{mui}, for ")
    (text . "i = 1,2,...,\htbitmap{ny} are chosen by the routine ")
    (text . ", but a single parameter S must be specified to control the ")
    (text . "trade-off between closeness of fit and smoothness of fit.  This ")
    (text . "affects the number of knots required by the spline, which is ")
    (text . "given in the B-spline representation \htbitmap{e02daf}")
    (text . "\newline \menuitemstyle{}\tab{2} \newline ")
    (text . "Number of data points {\it m}: ")
    (text . "\newline\tab{2} ")
    (bcStrings (6 30 m PI))
    (text . "\blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} ")
    (text . "Over-estimate of \htbitmap{nx} of knots \lambda of the ") 
    (text . "computed spline, {\it nxest}: \newline\tab{2} ")
    (bcStrings (6 14 nxest PI))
    (text . "\blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} ")
    (text . "Over-estimate of \htbitmap{ny} of knots \mu of the computed ") 
    (text . "spline, {\it nyest}: \newline\tab{2} ")
    (bcStrings (6 14 nyest PI))
    (text . "\blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} \newline ")
    (text . "Smoothing factor {\it s}: ")
    (text . "\newline\tab{2} ")
    (bcStrings (6 "10" s F))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{} \tab{2} ")
    (text . "\newline \tab{2} ")
    (text . "Start value: ")
    (radioButtons start
        ("" "  Cold Start - no values needed for {\it nx,ny,lamda,wrk,iwrk}" cold)
        ("" "  Warm Start - uses knots found in a previous call" warm))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{} \tab{2} ")
    (text . "\newline \tab{2} ")
    (text . "Ifail value: ")
    (radioButtons ifail
        ("" " -1, Print error messages" minusOne)
        ("" "  1, Suppress error messages" one)))
  htMakeDoneButton('"Continue", 'e02ddfSolve)
  htShowPage()

e02ddfSolve htPage ==
  m :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'm)
    objValUnwrap htpLabelSpadValue(htPage, 'm)
  nxest :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'nxest)
    objValUnwrap htpLabelSpadValue(htPage, 'nxest)
  nyest :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'nyest)
    objValUnwrap htpLabelSpadValue(htPage, 'nyest)
  u := nxest - 4
  v := nyest - 4
  wlist := [u,v]
  w := APPLY ('MAX, wlist)
  lwrk := (7*u*v + 25*w)*(w + 1) + 2*(u + v + 4*m) + 23*w + 56
  liwrk := m + 2*(nxest - 7)*(nyest - 7)
  s := htpLabelInputString(htPage,'s)
  initial := htpButtonValue(htPage,'start)
  start :=
    initial = 'cold => '1
    '2
  error := htpButtonValue(htPage,'ifail)
  ifail :=
    error = 'one => '1
    '-1
  (m = 30 and start = 1) => e02ddfDefaultSolve (htPage,nxest,nyest,lwrk,liwrk,s,ifail) 
  start = 1 => e02ddfColdSolve (htPage,m,nxest,nyest,lwrk,liwrk,s,ifail) 
  -- need to change as only wrk(1) is required
  labelList :=
    "append"/[f(i) for i in 1..m] where f(i) ==
      prefix := ('"\newline \tab{2} ")
      middle := ('"\tab{17} ")
      post   := ('"\tab{32} ")
      end    := ('"\tab{47} ")
      xnam := INTERN STRCONC ('"x",STRINGIMAGE i)
      ynam := INTERN STRCONC ('"y",STRINGIMAGE i)
      fnam := INTERN STRCONC ('"f",STRINGIMAGE i)
      wnam := INTERN STRCONC ('"w",STRINGIMAGE i)
      [['text,:prefix],['bcStrings,[8, 0.0, xnam, 'F]], 
       ['text,:middle],['bcStrings,[8, 0.0, ynam, 'F]],
        ['text,:post],['bcStrings,[8, 0.0, fnam, 'F]],
         ['text,:end],['bcStrings,[8, 0.0, wnam, 'F]]]
  lamdaList := 
    "append"/[g(i) for i in 1..nxest] where g(i) ==
      lnam := INTERN STRCONC ('"l",STRINGIMAGE i)
      [['bcStrings,[8, 0.0, lnam, 'F]]]
  prefix := ('"\blankline \menuitemstyle{} \tab{2} Values of Lamda: \newline")
  lamdaList := [['text,:prefix],:lamdaList]
  muList := 
    "append"/[h(i) for i in 1..nyest] where h(i) ==
      mnam := INTERN STRCONC ('"m",STRINGIMAGE i)
      [['bcStrings,[8, 0.0, mnam, 'F]]]
  prefix := ('"\blankline \menuitemstyle{} \tab{2} Values of Mu: \newline")
  muList := [['text,:prefix],:muList]
  prefix := ('"\blankline \menuitemstyle{} \tab{2} Value of nx: \newline \tab{2}")
  nxList := [['text,:prefix],['bcStrings,[8, 10, 'nx, 'PI]]]
  prefix := ('"\blankline \menuitemstyle{} \tab{2} Value of ny: \newline \tab{2}")
  nyList := [['text,:prefix],['bcStrings,[8, 9, 'ny, 'PI]]]
  prefix := ('"\blankline \menuitemstyle{} \tab{2} Value of WRK(1): \newline \tab{2}")
  wList := [['text,:prefix],['bcStrings,[8, 0.0, 'wone, 'F]]]
  equationPart := [
     '(domainConditions 
        (isDomain P (Polynomial $EmptyMode))
         (isDomain S (String))
          (isDomain F (Float))
           (isDomain PI (PositiveInteger))),
            :labelList,:lamdaList,:muList,:nxList,:nyList,:wList]
  page := htInitPage('"E02DDF - Least-squares surface fit by bicubic splines with automatic knot placement, scattered data",nil)
  htSay '"\menuitemstyle{}\tab{2} Values of \space{1} "
  htSay '"\htbitmap{xr}: \tab{15} \menuitemstyle{}\tab{17} "
  htSay '"Values of \space{1} \htbitmap{yr}: \tab{30}"
  htSay '"\menuitemstyle{}\tab{32} Values of \space{1} "
  htSay '"\htbitmap{fr}: \tab{45} \menuitemstyle{} "
  htSay '"\tab{47} Values of \htbitmap{wr}:"
  htMakePage equationPart
  htSay '"\blankline "
  htMakeDoneButton('"Continue",'e02ddfWarmGen)
  htpSetProperty(page,'m,m)
  htpSetProperty(page,'nxest,nxest)
  htpSetProperty(page,'nyest,nyest)
  htpSetProperty(page,'lwrk,lwrk)
  htpSetProperty(page,'liwrk,liwrk)
  htpSetProperty(page,'s,s)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()




e02ddfColdSolve(htPage,m,nxest,nyest,lwrk,liwrk,s,ifail)  ==
  labelList :=
    "append"/[f(i) for i in 1..m] where f(i) ==
      prefix := ('"\newline \tab{2} ")
      middle := ('"\tab{17} ")
      post   := ('"\tab{32} ")
      end    := ('"\tab{47} ")
      xnam := INTERN STRCONC ('"x",STRINGIMAGE i)
      ynam := INTERN STRCONC ('"y",STRINGIMAGE i)
      fnam := INTERN STRCONC ('"f",STRINGIMAGE i)
      wnam := INTERN STRCONC ('"w",STRINGIMAGE i)
      [['text,:prefix],['bcStrings,[8, 0.0, xnam, 'F]], 
       ['text,:middle],['bcStrings,[8, 0.0, ynam, 'F]],
        ['text,:post],['bcStrings,[8, 0.0, fnam, 'F]],
         ['text,:end],['bcStrings,[8, 0.0, wnam, 'F]]]
  equationPart := [
     '(domainConditions 
        (isDomain P (Polynomial $EmptyMode))
         (isDomain S (String))
          (isDomain F (Float))
           (isDomain PI (PositiveInteger))),
            :labelList]
  page := htInitPage('"E02DDF - Least-squares surface fit by bicubic splines with automatic knot placement, scattered data",nil)
  htSay '"\menuitemstyle{}\tab{2} Values of \space{1} "
  htSay '"\htbitmap{xr}: \tab{15} \menuitemstyle{}\tab{17} "
  htSay '"Values of \space{1} \htbitmap{yr}: \tab{30}"
  htSay '"\menuitemstyle{}\tab{32} Values of \space{1} "
  htSay '"\htbitmap{fr}: \tab{44} \menuitemstyle{} "
  htSay '"\tab{46} Values of \htbitmap{wr}:"
  htMakePage equationPart
  htSay '"\blankline "
  htMakeDoneButton('"Continue",'e02ddfColdGen)
  htpSetProperty(page,'m,m)
  htpSetProperty(page,'nxest,nxest)
  htpSetProperty(page,'nyest,nyest)
  htpSetProperty(page,'lwrk,lwrk)
  htpSetProperty(page,'liwrk,liwrk)
  htpSetProperty(page,'s,s)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()


e02ddfDefaultSolve (htPage,nxest,nyest,lwrk,liwrk,s,ifail) ==
  m := 30
  page := htInitPage('"E02DDF - Least-squares surface fit by bicubic splines with automatic knot placement, scattered data",nil)
  htMakePage '(
    (domainConditions 
       (isDomain F (Float)))
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2} Values of \htbitmap{xr}:")
    (text . "\tab{15} \menuitemstyle{} \tab{17} Values of ")
    (text . "\htbitmap{yr}: \tab{30} \menuitemstyle{} \tab{32} ")
    (text . "Values of \htbitmap{fr}: \tab{44} \menuitemstyle{} ")
    (text . "\tab{46} Values of \htbitmap{wr}:")
    (text . "\newline \tab{2} ")
    (bcStrings (8 "11.16" x1 F))  
    (text . "\tab{17}")
    (bcStrings (8 "1.24" y1 F))    
    (text . "\tab{32}")
    (bcStrings (8 "22.15" f1 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w1 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "12.85" x2 F))  
    (text . "\tab{17}")
    (bcStrings (8 "3.06" y2 F))    
    (text . "\tab{32}")
    (bcStrings (8 "22.11" f2 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w2 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "19.85" x3 F))  
    (text . "\tab{17}")
    (bcStrings (8 "10.72" y3 F))    
    (text . "\tab{32}")
    (bcStrings (8 "7.97" f3 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w3 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "19.72" x4 F))  
    (text . "\tab{17}")
    (bcStrings (8 "1.39" y4 F))    
    (text . "\tab{32}")
    (bcStrings (8 "16.83" f4 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w4 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "15.91" x5 F))  
    (text . "\tab{17}")
    (bcStrings (8 "7.74" y5 F))    
    (text . "\tab{32}")
    (bcStrings (8 "15.30" f5 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w5 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "0" x6 F))  
    (text . "\tab{17}")
    (bcStrings (8 "20" y6 F))    
    (text . "\tab{32}")
    (bcStrings (8 "34.6" f6 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w6 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "20.87" x7 F))  
    (text . "\tab{17}")
    (bcStrings (8 "20" y7 F))    
    (text . "\tab{32}")
    (bcStrings (8 "5.74" f7 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w7 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "3.45" x8 F))  
    (text . "\tab{17}")
    (bcStrings (8 "12.78" y8 F))    
    (text . "\tab{32}")
    (bcStrings (8 "41.24" f8 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w8 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "14.26" x9 F))  
    (text . "\tab{17}")
    (bcStrings (8 "17.87" y9 F))    
    (text . "\tab{32}")
    (bcStrings (8 "10.74" f9 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w9 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "17.43" x10 F))  
    (text . "\tab{17}")
    (bcStrings (8 "3.46" y10 F))    
    (text . "\tab{32}")
    (bcStrings (8 "18.60" f10 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w10 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "22.8" x11 F))  
    (text . "\tab{17}")
    (bcStrings (8 "12.39" y11 F))    
    (text . "\tab{32}")
    (bcStrings (8 "5.47" f11 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w11 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "7.58" x12 F))  
    (text . "\tab{17}")
    (bcStrings (8 "1.98" y12 F))    
    (text . "\tab{32}")
    (bcStrings (8 "29.87" f12 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w12 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "25" x13 F))  
    (text . "\tab{17}")
    (bcStrings (8 "11.87" y13 F))    
    (text . "\tab{32}")
    (bcStrings (8 "4.4" f13 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w13 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "0" x14 F))  
    (text . "\tab{17}")
    (bcStrings (8 "0" y14 F))    
    (text . "\tab{32}")
    (bcStrings (8 "58.2" f14 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w14 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "9.66" x15 F))  
    (text . "\tab{17}")
    (bcStrings (8 "20" y15 F))    
    (text . "\tab{32}")
    (bcStrings (8 "4.73" f15 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w15 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "5.22" x16 F))  
    (text . "\tab{17}")
    (bcStrings (8 "14.66" y16 F))    
    (text . "\tab{32}")
    (bcStrings (8 "40.36" f16 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w16 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "17.25" x17 F))  
    (text . "\tab{17}")
    (bcStrings (8 "19.57" y17 F))    
    (text . "\tab{32}")
    (bcStrings (8 "6.43" f17 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w17 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "25" x18 F))  
    (text . "\tab{17}")
    (bcStrings (8 "3.87" y18 F))    
    (text . "\tab{32}")
    (bcStrings (8 "8.74" f18 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w18 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "12.13" x19 F))  
    (text . "\tab{17}")
    (bcStrings (8 "10.79" y19 F))    
    (text . "\tab{32}")
    (bcStrings (8 "13.71" f19 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w19 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "22.23" x20 F))  
    (text . "\tab{17}")
    (bcStrings (8 "6.21" y20 F))    
    (text . "\tab{32}")
    (bcStrings (8 "10.25" f20 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w20 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "11.52" x21 F))  
    (text . "\tab{17}")
    (bcStrings (8 "8.53" y21 F))    
    (text . "\tab{32}")
    (bcStrings (8 "15.74" f21 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w21 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "15.2" x22 F))  
    (text . "\tab{17}")
    (bcStrings (8 "0" y22 F))    
    (text . "\tab{32}")
    (bcStrings (8 "21.6" f22 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w22 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "7.54" x23 F))  
    (text . "\tab{17}")
    (bcStrings (8 "10.69" y23 F))    
    (text . "\tab{32}")
    (bcStrings (8 "19.31" f23 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w23 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "17.32" x24 F))  
    (text . "\tab{17}")
    (bcStrings (8 "13.78" y24 F))    
    (text . "\tab{32}")
    (bcStrings (8 "12.11" f24 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w24 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "2.14" x25 F))  
    (text . "\tab{17}")
    (bcStrings (8 "15.03" y25 F))    
    (text . "\tab{32}")
    (bcStrings (8 "53.1" f25 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w25 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "0.51" x26 F))  
    (text . "\tab{17}")
    (bcStrings (8 "8.37" y26 F))    
    (text . "\tab{32}")
    (bcStrings (8 "49.43" f26 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w26 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "22.69" x27 F))  
    (text . "\tab{17}")
    (bcStrings (8 "19.63" y27 F))    
    (text . "\tab{32}")
    (bcStrings (8 "3.25" f27 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w27 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "5.47" x28 F))  
    (text . "\tab{17}")
    (bcStrings (8 "17.13" y28 F))    
    (text . "\tab{32}")
    (bcStrings (8 "28.63" f28 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w28 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "21.67" x29 F))  
    (text . "\tab{17}")
    (bcStrings (8 "14.36" y29 F))    
    (text . "\tab{32}")
    (bcStrings (8 "5.52" f29 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w29 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "3.31" x30 F))  
    (text . "\tab{17}")
    (bcStrings (8 "0.33" y30 F))    
    (text . "\tab{32}")
    (bcStrings (8 "44.08" f30 F))    
    (text . "\tab{47}")
    (bcStrings (8 "1" w30 F))    
    (text . "\blankline"))
  htMakeDoneButton('"Continue",'e02ddfColdGen)          
  htpSetProperty(page,'m,m)
  htpSetProperty(page,'nxest,nxest)
  htpSetProperty(page,'nyest,nyest)
  htpSetProperty(page,'lwrk,lwrk)
  htpSetProperty(page,'liwrk,liwrk)
  htpSetProperty(page,'s,s)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()

e02ddfColdGen htPage ==
  m := htpProperty(htPage,'m)
  nxest := htpProperty(htPage,'nxest)
  nyest := htpProperty(htPage,'nyest)
  lwrk := htpProperty(htPage,'lwrk)
  liwrk := htpProperty(htPage,'liwrk)
  s := htpProperty(htPage,'s)
  cold := '"c"
  ifail := htpProperty(htPage,'ifail)
  alist := htpInputAreaAlist htPage
  y := alist
  while y repeat
    end :=  STRCONC ((first y).1," ")
    y := rest y
    right := STRCONC ((first y).1," ")
    y := rest y    
    mid :=  STRCONC ((first y).1," ")
    y := rest y
    left :=  STRCONC ((first y).1," ")
    y := rest y
    xlist := [left,:xlist]
    ylist := [mid,:ylist]
    flist := [right,:flist]
    wlist := [end,:wlist]
  xstring := bcwords2liststring xlist
  ystring := bcwords2liststring ylist
  fstring := bcwords2liststring flist
  wstring := bcwords2liststring wlist
  -- additional entries nx,ny,lamda,mu,wrk needed to get it running
  -- but they are just set to 0.0
  prefix := STRCONC('"e02ddf(_"",cold,"_",",STRINGIMAGE m,", [",xstring,"],[")
  prefix := STRCONC(prefix,ystring,"],[",fstring,"],[",wstring,"], ")
  prefix := STRCONC(prefix,STRINGIMAGE s,", ",STRINGIMAGE nxest,", ")
  prefix := STRCONC(prefix,STRINGIMAGE nyest,", ",STRINGIMAGE lwrk,", ")
  prefix := STRCONC(prefix,STRINGIMAGE liwrk,", 0,")
  prefix := STRCONC(prefix,"new(1,", STRINGIMAGE nxest,",0.0)$MATRIX DFLOAT,0,")
  prefix := STRCONC(prefix,"new(1,", STRINGIMAGE nyest,",0.0)$MATRIX DFLOAT,")
  prefix := STRCONC(prefix,"new(1,", STRINGIMAGE lwrk,",0.0)$MATRIX DFLOAT,")
--  prefix := STRCONC(prefix,"[[0.0 for i in 1..", STRINGIMAGE nxest,"]],0,")
--  prefix := STRCONC(prefix,"[[0.0 for i in 1..", STRINGIMAGE nyest,"]],")
--  prefix := STRCONC(prefix,"[[0.0 for i in 1..", STRINGIMAGE lwrk,"]],")
  prefix := STRCONC(prefix,STRINGIMAGE ifail,")")
  linkGen prefix

e02ddfWarmGen htPage ==
  m := htpProperty(htPage,'m)
  nxest := htpProperty(htPage,'nxest)
  nyest := htpProperty(htPage,'nyest)
  lwrk := htpProperty(htPage,'lwrk)
  liwrk := htpProperty(htPage,'liwrk)
  s := htpProperty(htPage,'s)
  warm := '"w"
  ifail := htpProperty(htPage,'ifail)
  alist := htpInputAreaAlist htPage
  y := alist
  wrk := (first y).1
  y := rest y
  for i in 1..lwrk repeat
    wrkList := ['"0.0 ",:wrkList]
  wrkList := [wrk,:wrkList]
  wrkstring := bcwords2liststring wrkList
  ny := STRCONC((first y).1," ")
  y := rest y
  nx := STRCONC((first y).1," ")
  y := rest y
  for i in 1..nyest repeat
    mu := STRCONC ((first y).1, " ")
    y := rest y
    muList := [mu,:muList]
  mustring := bcwords2liststring muList
  for i in 1..nxest repeat
    lam := STRCONC ((first y).1, " ")
    y := rest y
    lamList := [lam,:lamList]
  lamstring := bcwords2liststring lamList
  while y repeat
    end :=  STRCONC ((first y).1," ")
    y := rest y
    right := STRCONC ((first y).1," ")
    y := rest y    
    mid :=  STRCONC ((first y).1," ")
    y := rest y
    left :=  STRCONC ((first y).1," ")
    y := rest y
    xlist := [left,:xlist]
    ylist := [mid,:ylist]
    flist := [right,:flist]
    wlist := [end,:wlist]
  xstring := bcwords2liststring xlist
  ystring := bcwords2liststring ylist
  fstring := bcwords2liststring flist
  wstring := bcwords2liststring wlist
  -- additional entries nx,ny,lamda,mu,wrk needed to get it running
  -- but they are just set to 0.0
  prefix := STRCONC('"e02ddf(_"",warm,"_",",STRINGIMAGE m,", [",xstring,"],[")
  prefix := STRCONC(prefix,ystring,"],[",fstring,"],[",wstring,"], ")
  prefix := STRCONC(prefix,STRINGIMAGE s,", ",STRINGIMAGE nxest,", ")
  prefix := STRCONC(prefix,STRINGIMAGE nyest,", ",STRINGIMAGE lwrk,", ")
  prefix := STRCONC(prefix,STRINGIMAGE liwrk,", ",nx,",[",lamstring,"],",ny)
  prefix := STRCONC(prefix,",[",mustring,"],[",wrkstring,"],")
  prefix := STRCONC(prefix,STRINGIMAGE ifail,")")
  linkGen prefix

e02zaf() == 
  htInitPage('"E02ZAF - Sort 2-D sata into panels for fitting bicubic splines",nil)
  htMakePage '(
    (domainConditions 
       (isDomain PI (PositiveInteger)))
    (text . "\windowlink{Manual Page}{manpageXXe02zaf} for this routine ")
    (text . "\newline ")
    (text . "\lispwindowlink{Browser operation page}{(|oPageFrom| '|e02zaf| '|NagFittingPackage|)} for this routine")
    (text . "\newline \horizontalline ")
    (text . "Sorts the set of points {\em (\htbitmap{xr},")
    (text . "\htbitmap{yr})} into panels defined by \space{1}")
    (text . "\htbitmap{px} -8 points \htbitmap{lamdai} ")
    (text . "on the x-axis and \space{1}\htbitmap{py}-8 points ")
    (text . "\htbitmap{muj} on the y axis. The points are ordered ")
    (text . "so that all points in a panel occur before data in succeeding ")
    (text . "panels. Within a panel, the points maintain their original ")
    (text . "order. ") 
    (text . "\blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} \newline ")
    (text . "Number of points to be sorted to be sorted {\it m}: ")
    (text . "\newline\tab{2} ")
    (bcStrings (6 10 m PI))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2}")
    (text . "\newline Intercepts + 8 on x axis {\em px}:")
    (text . "\tab{32} \menuitemstyle{}\tab{34}")
    (text . "Intercepts + 8 on y axis {\em py}:")
    (text . "\newline\tab{2} ")
    (bcStrings (6 9 px PI))
    (text . "\tab{34} ")
    (bcStrings (6 10 py PI))
    (text . "\blankline ")
    (text . "\newline \menuitemstyle{}\tab{2} \newline ")
    (text . "Dimension of point {\it npoint}: ")
    (text . "\newline\tab{2} ")
    (bcStrings (6 45 npoint PI))
    (text . "\blankline ")
    (text . "\newline ")
    (text . "\menuitemstyle{} \tab{2} ")
    (text . "\newline \tab{2} ")
    (text . "Ifail value: ")
    (radioButtons ifail
        ("" " -1, Print error messages" minusOne)
        ("" "  1, Suppress error messages" one)))
  htMakeDoneButton('"Continue", 'e02zafSolve)
  htShowPage()

e02zafSolve htPage ==
  m :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'm)
    objValUnwrap htpLabelSpadValue(htPage, 'm)
  px :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'px)
    objValUnwrap htpLabelSpadValue(htPage, 'px)
  py :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'py)
    objValUnwrap htpLabelSpadValue(htPage, 'py)
  npoint :=
    $bcParseOnly => PARSE_-INTEGER htpLabelInputString(htPage, 'npoint)
    objValUnwrap htpLabelSpadValue(htPage, 'npoint)
  error := htpButtonValue(htPage,'ifail)
  ifail :=
    error = 'one => '1
    '-1
  ((m = '10 and px = '9) and py = '10) => e02zafDefaultSolve(htPage,npoint,ifail)
  labelList :=
    "append"/[fxy(i) for i in 1..m] where fxy(i) ==
      prefix := ('"\newline \tab{2} ")
      middle := ('"\tab{32} ")
      lnam := INTERN STRCONC ('"x",STRINGIMAGE i)
      cnam := INTERN STRCONC ('"y",STRINGIMAGE i)
      [['text,:prefix],['bcStrings,[8, 0.0, lnam, 'F]], 
       ['text,:middle],['bcStrings,[8, 0.0, cnam, 'F]]]
  lamList :=
    "append"/[flam(i) for i in 5..(px-4)] where flam(i) ==
      lnam := INTERN STRCONC ('"l",STRINGIMAGE i)
      [['bcStrings,[8, 0.0, lnam, 'F]]]
  prefix := ('"\blankline \menuitemstyle{} \tab{2} {\it \lambda(5) to ")
  prefix := STRCONC(prefix,"\lambda(px-4)}: \newline \tab{2} ")
  lamList := [['text,:prefix],:lamList]
  muList :=
    "append"/[fmu(i) for i in 5..(py-4)] where fmu(i) ==
      mnam := INTERN STRCONC ('"m",STRINGIMAGE i)
      [['bcStrings,[8, 0.0, mnam, 'F]]]
  prefix := ('"\blankline \menuitemstyle{} \tab{2} {\it \mu(5) to \mu(py-4)}: ")
  prefix := STRCONC(prefix,"\newline \tab{2} ")
  muList := [['text,:prefix],:muList]
  equationPart := [
     '(domainConditions 
        (isDomain P (Polynomial $EmptyMode))
         (isDomain S (String))
          (isDomain F (Float))
           (isDomain PI (PositiveInteger))),
            :labelList,:lamList,:muList]
  page := htInitPage('"E02ZAF - Sort 2-D sata into panels for fitting bicubic splines",nil)
  htSay '"\menuitemstyle{}\tab{2} {\it x(m)}: "
  htSay '"\tab{30} \menuitemstyle{}\tab{32} {\it y(m)}: "
  htMakePage equationPart
  htSay '"\blankline "
  htMakeDoneButton('"Continue",'e02zafGen)
  htpSetProperty(page,'m,m)
  htpSetProperty(page,'px,px)
  htpSetProperty(page,'py,py)
  htpSetProperty(page,'npoint,npoint)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()


e02zafDefaultSolve (htPage,npoint,ifail) ==
  m := '10
  px := '9
  py := '10
  page := htInitPage('"E02ZAF - Sort 2-D sata into panels for fitting bicubic splines",nil)
  htMakePage '(
    (domainConditions 
       (isDomain F (Float)))
    (text . "\newline ")
    (text . "\menuitemstyle{}\tab{2} {\it x(m)}:")
    (text . "\tab{30} \menuitemstyle{} \tab{32} {\it y(m)}:")
    (text . "\newline \tab{2} ")
    (bcStrings (8 "0.00" x1 F))  
    (text . "\tab{32}")
    (bcStrings (8 "0.77" y1 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "0.70" x2 F))  
    (text . "\tab{32}")
    (bcStrings (8 "1.06" y2 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "1.44" x3 F))  
    (text . "\tab{32}")
    (bcStrings (8 "0.33" y3 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "0.21" x4 F))  
    (text . "\tab{32}")
    (bcStrings (8 "0.44" y4 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "1.01" x5 F))  
    (text . "\tab{32}")
    (bcStrings (8 "0.50" y5 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "1.84" x6 F))  
    (text . "\tab{32}")
    (bcStrings (8 "0.02" y6 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "0.71" x7 F))  
    (text . "\tab{32}")
    (bcStrings (8 "1.95" y7 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "1.00" x8 F))  
    (text . "\tab{32}")
    (bcStrings (8 "1.20" y8 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "0.54" x9 F))  
    (text . "\tab{32}")
    (bcStrings (8 "0.04" y9 F))    
    (text . "\newline \tab{2} ")
    (bcStrings (8 "1.531" x10 F))  
    (text . "\tab{32}")
    (bcStrings (8 "0.18" y10 F))    
    (text . "\blankline ")
    (text . "\menuitemstyle{}\tab{2} {\it \lambda(5) to \lambda(px-4)}:")
    (text . "\newline \tab{2}")
    (bcStrings (8 "1.00" l5 F))
    (text . "\blankline ")
    (text . "\menuitemstyle{}\tab{2} {\it \mu(5) to \mu(py-4)}:")
    (text . "\newline \tab{2}")
    (bcStrings (8 "0.80" mu5 F))
    (bcStrings (8 "1.20" mu6 F))
    (text . "\blankline "))
  htMakeDoneButton('"Continue",'e02zafGen)      
  htpSetProperty(page,'m,m)
  htpSetProperty(page,'px,px)
  htpSetProperty(page,'py,py)
  htpSetProperty(page,'npoint,npoint)
  htpSetProperty(page,'ifail,ifail)
  htpSetProperty(page,'inputArea, htpInputAreaAlist htPage)
  htShowPage()

e02zafGen htPage ==
  m := htpProperty(htPage,'m)
  px := htpProperty(htPage,'px)
  py := htpProperty(htPage,'py)
  npoint := htpProperty(htPage,'npoint)
  ifail := htpProperty(htPage,'ifail)
  alist := htpInputAreaAlist htPage
  y := alist
  nadres := (px-7)*(py-7)
  -- mu 
  for i in 1..4 repeat
    muList := ['"0 ",:muList]
  for i in 5..(py-4) repeat
    right := STRCONC ((first y).1," ")
    y := rest y    
    muList := [right,:muList]
  for i in (py-3)..py repeat
    muList := ['"0 ",:muList]
  mustring := bcwords2liststring muList
  -- lamda
  for i in 1..4 repeat
    lamList := ['"0 ",:lamList]
  for i in 5..(px-4) repeat
    right := STRCONC ((first y).1," ")
    y := rest y    
    lamList := [right,:lamList]
  for i in (px-3)..px repeat
    lamList := ['"0 ",:lamList]
  lamstring := bcwords2liststring lamList
  -- x & y
  while y repeat
    right := STRCONC ((first y).1," ")
    y := rest y    
    left :=  STRCONC ((first y).1," ")
    y := rest y
    xlist := [left,:xlist]
    ylist := [right,:ylist]
  xstring := bcwords2liststring xlist
  ystring := bcwords2liststring ylist
  prefix := STRCONC('"e02zaf(",STRINGIMAGE px,", ",STRINGIMAGE py,",[")
  prefix := STRCONC(prefix,lamstring,"],[",mustring,"], ",STRINGIMAGE m,", [")
  prefix := STRCONC(prefix,xstring,"],[",ystring,"], ",STRINGIMAGE npoint,", ")
  prefix := STRCONC(prefix,STRINGIMAGE nadres,", ",STRINGIMAGE ifail,")")
  linkGen prefix