#include double _j1(double x) { return x; } double _y1(double x) { return x; }