#include
double?y（int?n）;
double?x（int?n）;
double?X[41]Y[41];

double?yp（int?n）
{
???return?Y[n-1]+0.1*（1/（1+X[n-1]*X[n-1]）-2*Y[n-1]*Y[n-1]）;????	
}

double?yc（int?n）
{
????return??Y[n-1]+0.1*（1/（1+X[n]*X[n]）-2*yp（n）*yp（n））;
}

double?y（int?n）
{
	if（n==0）