// Test_poly.scala    test  polyval  and  polyfit

object Test_poly
{
  var P = Poly

  def main(args: Array[String])
  {
    println("Test_poly.scala running")
    var pwr = 3
    var pn = pwr+1
    var n = 10
    var x =  Array.ofDim[Double](n)
    var y =  Array.ofDim[Double](n)
    var yc = Array.ofDim[Double](n)
    var pc:Array[Double] = Array(1.2, 2.4, 3.6, 4.8)
    var p =  Array.ofDim[Double](pn)
    var err = 0.0

    for(i <- 0 until n)
    {
      x(i) = (i+1)*0.25
      println("i="+i+", x(i)="+x(i))
    }

    for(i <- 0 until n)
    {
      y(i) = P.polyval(pwr, pc, x(i))
      yc(i) = y(i)
      println("i="+i+", x(i)="+x(i)+", y(i)="+y(i))
    }

    println("computed values")

    p = P.polyfit(pwr, n, x, y)
    for(j <- 0 until pn) println("j="+j+", p(j)="+p(j))

    for(i <- 0 until n)
    {
      y(i) = P.polyval(pwr, p, x(i))
      err = (y(i)-yc(i)).abs
      println("i="+i+", x(i)="+x(i)+", y(i)="+y(i)+", err="+err)
    }
    println("Test_poly.scala finished")
  }
} // end Test_poly.scala