краевая задача диф.уравнения методом стрельбы

[Laplas]

добрый вечер!!! нужно численно решить диф.уравнение методом стрельбы:
y"=y+2*x;
краевые условия
y(0)=0;
y(1)=-1;
я сделал так:

#include "stdafx.h"
#include "stdio.h"
#include "math.h"
#include "conio.h"

const int n = 1000;

double func2(double x, double y, double z)
{
 	return 2*x+y;
}
double func1(double x, double y, double z)
{ 
	return z;
}
double RK (double x, double y, double z, double a, double b)
{ 
	double k1, k2, k3, k4,l1,l2,l3,l4, h;
	int i;
		h=(b-a)/n;
	FILE *f;
f=fopen ("1.txt", "w");

		for (i=1; i<n; i++)
		{x=i*h;
		k1=h*func1(x,y,z); 
		l1=h*func2(x,y,z);
		k2=h*func1(x+h/2,y+k1/2,z+l1/2);
		l2=h*func2(x+h/2,y+k1/2,z+l1/2);
		k3=h*func1(x+h/2,y+k2/2,z+l2/2);
		l3=h*func2(x+h/2,y+k2/2,z+l2/2);
		k4=h*func1(x+h,y+k3,z+l3);
		l4=h*func2(x+h,y+k3,z+l3);
		y+=(k1+2*k2+2*k3+k4)*h/6;
		z+=(l1+2*l2+2*l3+l4)*h/6;
		fprintf (f,"%lf\t%lf\t%lf\n",x,y,z);
		}
	fclose (f);
	return y;
}
int _tmain(int argc, _TCHAR* argv[])
{double a=0, b=1, ya=0, yb=-1, z0,z1, Fz0, dFz0, eps=0.001, delta;
 scanf ("%lf",&z0);
	do 
		{
		Fz0=yb-RK(a,ya,z0,a,b);
		dFz0=(RK(a,ya,(z0+eps),a,b)-yb-Fz0)/eps;
		z1=z0+Fz0/dFz0;
		delta=fabs(z1-z0);
		z0=z1;
		}
	while (delta>eps);
    getch();
	return 0;
}

но она почему неправильно считает z=y'(x),
посмотрите пожалуйста, может я ошибся где то????))))
спасибо))

[Sweta]

В методе понимаю мало, но смотря на Вашу программу хочу спросить зачем нужна эта функция?

double func1(double x, double y, double z)
{ 
	return z;
}

Не обижайтесь. Может отвечая на мои глупые вопросы сможете найти ошибку.

[Serebro]

dFz0=(RK(a,ya,(z0+eps),a,b)-yb-Fz0)/eps;

Здесь ошибка.