sl_ex.cpp

矩阵计算库

// This is an example of the use of solution to find the cube root of 
// the integers -10 to 10

// you will need to compile and link solution.cpp and except.cpp

#define WANT_STREAM
#define WANT_MATH

#include "include.h"
#include "solution.h"

// the cube class

class Cube : public R1_R1
{ Real operator()() { return x*x*x; } };

main()
{
   // construct the Cube object
   Cube cube;
   // and then the solve object
   OneDimSolve cube_root(cube);
   // Now do the solves
   for (int i=-10; i<=10; i++)
      cout << i << "   "  << cube_root.Solve(i,0,1.5) << endl;
   return 0;
}