常微分方程数值解Euler法c语言

用c语言实现Euler法，附带详细使用说明文档

#include
#include
#define N 20

main（）
{
	float h;//存放步长
	float ab;//存放区间端点值
	int n;
	float t0u0;//存放初值
	float u[N];//存放结果

	void input（float *afloat *bfloat *t0float *u0float *h）;//从文件输入abt0u0h
	int get_n（float afloat bfloat h）;//得到n
	void Euler（float t0float u0float hint nfloat u[]）;//返回结果，并存放于r[N]
	void output（float t0float hfloat u[]int n）;//输出结果

	input（&a&b&t0&u0&h）;

	printf（“%f%f%f%f%f“abt0u0h）;//测试输入
	printf（“\n“）;

	//n=get_n（abh）;
	n=10;
	Euler（t0u0hnu）;
	output（t0hun）;
}

//从文件输入abt0u0h
void input（float *pafloat *pbfloat *pt0float *pu0float *ph）
{
	FILE *f_a_b*f_t0_u0*f_h;
	f_a_b=fopen（“a_b.dat““r“）;
	f_t0_u0=fopen（“t0_u0.dat““r“）;
	f_h=fopen（“h.dat““r“）;

	fscanf（f_a_b“%

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        121  2009-04-02 10:08  readme.txt

     文件          3  2009-03-05 22:14  t0_u0.dat

     文件          3  2009-03-05 22:13  a_b.dat

     文件       1631  2009-03-06 21:34  Euler.c

     文件          3  2009-03-05 22:14  h.dat

----------- ---------  ---------- -----  ----

                 1761                    5