#include<stdio.h>
#include<math.h>

/******************* 図6.5の数値解を求めるプログラム *******************/

#define a  3.4		//α
#define z0 0.070588	//z_0
#define imax 10		//z_imax(iの最大値)

double f(double zi){
	double z;
	z = a*zi*(1-zi);
	return z;
}

int main(void){
	int i;
	double z1, z2;

	z1 = z0;
	printf("       z[i]           z[i+1]\n");
	for(i=0; i<imax+1; i++){
		z2 = f(z1);
		printf("z[%3d]=%9.5lf, z[%3d]=%9.5lf\n", i, z1, i+1, z2);
		if(fabs(z2-z1) > pow(10, 15)){
			printf("overflow error!!\n");
			return 1;
		}
		z1 = z2;
	}

	return 0;
}