数值分析实验之改进的欧拉公式程序代码

#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）