double incriment(double l, double e, double i, double y4, double y3, double y2, double y1, double y0) {
	double yStar;
	double x = pow((l / 554), 4);
	yStar = ((2 * 16 * (x)) / (e * i)) + 3 * y4 + 2 * y3 - 2 * y2 + 3 * y1 + y0;
	return yStar;
}

int main()
{
	double e, q, l, y0, y1, y2, y3, y4, y5;

	e = 5; q = 10; l = 213; y0 = 0; y1 = 1; y2 = 2; y3 = 3; y4 = 4;

	//y5 = incriment(l, e, 1, y4, y3, y2, y1, y0);
	for (int i = 1; i < 551; i++) {

		y5 = incriment(l, e, i, y4, y3, y2, y1, y0);
		std::cout <<"Y "<<i+4<<" = "<< y5 << "\n";
		y0=y1;
		y1=y2;
		y2=y3;
		y4=y5;
	}

}