point_collocation.cpp

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

#include "vs.h"
int main() {
	double x[2] = {0.25, 0.75};
   C0 M(2, 2, (double*)0), b(2, x);
   for(int i = 0; i < 2; i++) {
   	C2 z(x[i]),
      	phi = VECTOR_OF_TANGENT_OF_TANGENT_BUNDLE("int, int", 2, 1);
      phi[0] = z*(1-z); phi[1] = z.pow(2)*(1-z);
      M[i] = (+dd(phi))(0) +((C0)phi);
   }
   C0 c = -b / M;
   cout << c << endl;
	return 0;
}