>> 201
あ。これy_exactの算出が間違ってるな
データプロットしたらそれっぽくなったんで満足してしまった　すまん
下記が修正パッチ

--- a.c 2015-06-07 19:14:32.946483097 +0900
+++ b.c 2015-06-07 19:15:01.269484468 +0900
@@ -98,3 +98,3 @@

- INIT_FUNC[M](N1-1, x1, y1);
+ INIT_FUNC[M](N1-1, x1, y_exact);

@@ -138,3 +138,2 @@
for (i = 0; i < N1; ++i) {
- y_exact[i] = sin(x1[i]);
fprintf(fp, "%f, %f, %f\n", x1[i], y1[i], y_exact[i]);