KM2016 / Elliptic / fem4

#include <ctype.h>
#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#include <string.h>

 static double funD( double xi, double bi []);
 static double funC( double xi, double bi []);
 static double funA( double xi, double bi []);
 static double fun_q( double xi, double bi []);

  void K_ij(int N, double K [], double x [], double bi [], double a, 
	  double (*funD)( double xi, double bi []), 
	  double (*funC)( double xi, double bi []), 
	  double (*funA)( double xi, double bi []));
 void q_i(int , double [], double [], double [], double (*)( double, double *));
 double domik(double x, int j, int mp, double *xx);
 double intgD(int, int, double *, double *, double (*)(double, double *));
 double intgCq(int, int, int , double *, double *, double (*)(double, double *), double (*)(double, int, int, double *));
 double intgA(int, int ,int, int , int , double *, double *, double (*)(double, double *), double (*)(double, int, int, double *));
 int i_max(int N, int p, double V[]);
 int lufact(int, double [], double [], double []);
 void SaveToFile(int, double*, const char* ); 
 int output(int test, int N, double* x, double *b, double *solution, const char* filename);