quadratic_coordinate_transformation.cpp

来自「算法的一些集合」· C++ 代码 · 共 16 行

#include "vs.h"
static double f(double x) { return (x-1.0)*(x-1.0); }
int main() {
	double X[3] = {1.0, 1.5, 2.0};
   Quadrature qp(1,4);
   H0 z(qp),
   	N(3, (double*)0, qp);
   N[0] = (z-1.0)*z/2.0;
   N[1] = (1.0-z)*(1.0+z);
   N[2] = z*(z+1.0)/2.0;
   H0 y = N[0]*f(X[0])+N[1]*f(X[1])+N[2]*f(X[2]);
   C0 area = y | J(0.5);
   cout << area << endl;
	return 0;
}