#include?“stdio.h“???
#include?“stdlib.h“???
#include????

int?Func（double?y[]?double?d[]）
{
	d[0]?=?-?y[1]?-?y[2]?;??????/*y0‘=y1*/
	d[1]?=?y[0]?+?0.2?*?y[1];?????/*y1‘=y0*/
	d[2]?=?0.2?+?（y[1]?-?5.7）?*?y[2];?????/*y2‘=y2*/
	return（1）;
}

void?Euler1（double?t?double?y[]?int?n?double?h?int?k?double?z[3][11]）
/*int?n;??????/*整型變量，微分方程組中方程的個(gè)數(shù)，也是未知函數(shù)的個(gè)數(shù)
int?k;??????/*整型變量。積分步數(shù)（包括起始點(diǎn)這一步）
double?t;???????/*雙精度實(shí)型變量對(duì)微分方程進(jìn)行積分的起始點(diǎn)t0
double?h;???????/*雙精度實(shí)型變量。積分步長(zhǎng)
double?y[];?/*雙精度實(shí)型一維數(shù)組，長(zhǎng)度為n。存放n個(gè)未知函數(shù)yi在起始點(diǎn)t0處的函數(shù)值
double?z[];?/*雙精度實(shí)型二維數(shù)組，體積為nxk。返回k個(gè)積分點(diǎn)（包括起始點(diǎn)）上的未知函數(shù)值
*/
{
	int?i?j;
	double?*d;
	d?=?（double?*）malloc（n*sizeof（double））;
	if?（d?==?NULL）
	{
		printf（“內(nèi)存分配失敗\n“）;
		exit（1）;
	}
	/*將方程組的初值賦給數(shù)組z[i*k]*/
	for?（i?=?0;?i?<=?n?-?1;?i++）
		z[i][0]?=?y[i];
	for?（j?=?1;?j?<=?k?-?1;?j++）
	{
		Func（y?d）;??????????/*求出f（x）*/
		for?（i?=?0;?i?<=?n?-?1;?i++）
			y[i]?=?z[i][j?-?1]?+?h*d[i];

		Func（y?d）;
		for?（i?=?0;?i?<=?n?-?1;?i++）
			d[i]?=?z[i][j?-?1]?+?h*d[i];
		for?（i?=?0;?i?<=?n?-?1;?i++）
		{
			y[i]?=?（y[i]?+?d[i]）?/?2.0;
			z[i][j]?=?y[i];
		}
	}
	free（d）;
	return;
}
int?Savedata（double?xx[11]?double?Ynnn[3][11]）
{
	FILE?*fp;//文件指針??
	int?i?j;

	/*文件的打開(kāi)*/
	fp?=?fopen（“Euler_data.txt“?“w“）;//fopen打開(kāi)文件，這個(gè)文件可以是當(dāng)前不存在的。“w”以寫(xiě)入的形式打開(kāi)，“r”以讀的形式打開(kāi)??
	if?（fp?==?NULL）?//判斷如果文件指針為空??
	{
		printf（“File?cannot?open!?“）;
		exit（0）;//在以0的形式退出，必須在文件開(kāi)頭有#include?stdlib?頭文件即standard?library標(biāo)準(zhǔn)庫(kù)頭文件??
	}

	//寫(xiě)入東西??
	for?（j?=?0;?j?<=?9;?j++）
		fprintf（fp?“%f\t“?xx[j]）;//寫(xiě)入指針fp，寫(xiě)入的東西就是剛才的用戶(hù)輸入的d注意這里的fp和d沒(méi)有引號(hào)??
	fprintf（fp?“%f\n“?xx[10]）;
	for?（i?=?0;?i?<=?2;?i++）
	{
		for?（j?=?0;?j?<=?9;?j++）
			fprintf（fp?“%f\t“?Ynnn[i][j]）;//寫(xiě)入指針fp，寫(xiě)入的東西就是剛才的用戶(hù)輸入的d注意這里的fp和d沒(méi)有引號(hào)??
		fprintf（fp?“%f\n“?Ynnn[i][10]）;
	}
	//關(guān)閉文件??
	fclose（fp）;
	return?0;
}

int?main（）
{
	int?i?j;
	double?y[3]?z[3][11]?t?h?x[11];
	y[0]?=?0;??????????/*初值y0（0）=0*/
	y[1]?=?2;???????????/*初值y1（0）=2*/
	y[2]?=?2;???????????/*初值y2（0）=2*/
	t?=?0.0;??????????????/*起始點(diǎn)t=0*/
	h?=?0.1;?????????????/*步長(zhǎng)為0.01*/
	Euler1（t?y?3?h?11?z）;
	printf（“定步長(zhǎng)歐拉法結(jié)果：\n“）;
	for?（i?=?0;?i?<=?10;?i++）
	{
		x[i]?=?i*h;
		printf（“t=%f\t???“?x[i]）;
		for?（j?=?0;?j?<=?2;?j++）
			printf（“y（%d）=%f\t??“?j?z[j][i]）;
		printf（“\n“）;
	}
	Savedata（x?z）;
	getchar（）;
	return?0;
	
}

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----
?????目錄???????????0??2018-06-11?11:09??Euler\Euler\
?????目錄???????????0??2018-06-11?11:08??Euler\Euler\Debug\
?????文件???????33280??2018-05-16?11:53??Euler\Euler\Debug\Euler.exe
?????文件??????244784??2018-05-16?11:53??Euler\Euler\Debug\Euler.ilk
?????文件??????593920??2018-05-16?11:53??Euler\Euler\Debug\Euler.pdb
?????文件?????3473408??2018-05-16?12:07??Euler\Euler\Euler.sdf
?????文件?????????961??2018-05-15?15:04??Euler\Euler\Euler.sln
?????文件???????24064??2018-05-16?12:07??Euler\Euler\Euler.v12.suo
?????目錄???????????0??2018-06-11?11:09??Euler\Euler\Euler\
?????目錄???????????0??2018-06-11?11:09??Euler\Euler\Euler\Debug\
?????文件????????1765??2018-05-16?11:53??Euler\Euler\Euler\Debug\Euler.log
?????文件???????12784??2018-05-16?11:53??Euler\Euler\Euler\Debug\Euler.obj
?????目錄???????????0??2018-06-11?11:09??Euler\Euler\Euler\Debug\Euler.tlog\
?????文件????????2586??2018-05-16?11:53??Euler\Euler\Euler\Debug\Euler.tlog\cl.command.1.tlog
?????文件????????1900??2018-05-16?11:53??Euler\Euler\Euler\Debug\Euler.tlog\CL.read.1.tlog
?????文件????????1548??2018-05-16?11:53??Euler\Euler\Euler\Debug\Euler.tlog\CL.write.1.tlog
?????文件?????????188??2018-05-16?11:53??Euler\Euler\Euler\Debug\Euler.tlog\Euler.lastbuildstate
?????文件????????2626??2018-05-16?11:53??Euler\Euler\Euler\Debug\Euler.tlog\link.command.1.tlog
?????文件????????2524??2018-05-16?11:53??Euler\Euler\Euler\Debug\Euler.tlog\link.read.1.tlog
?????文件?????????628??2018-05-16?11:53??Euler\Euler\Euler\Debug\Euler.tlog\link.write.1.tlog
?????文件????????3459??2018-05-16?11:30??Euler\Euler\Euler\Debug\stdafx.obj
?????文件??????109568??2018-05-16?11:53??Euler\Euler\Euler\Debug\vc120.idb
?????文件???????86016??2018-05-16?11:53??Euler\Euler\Euler\Debug\vc120.pdb
?????文件????????2692??2018-05-16?11:53??Euler\Euler\Euler\Euler.cpp
?????文件????????4233??2018-05-15?15:04??Euler\Euler\Euler\Euler.vcxproj
?????文件????????1312??2018-05-15?15:04??Euler\Euler\Euler\Euler.vcxproj.filters
?????文件?????????410??2018-05-16?12:04??Euler\Euler\Euler\Euler_data.txt
?????文件????????1496??2018-05-15?15:04??Euler\Euler\Euler\ReadMe.txt
?????文件?????????211??2018-05-15?15:04??Euler\Euler\Euler\stdafx.cpp
?????文件?????????234??2018-05-15?15:04??Euler\Euler\Euler\stdafx.h
?????文件?????????236??2018-05-15?15:04??Euler\Euler\Euler\targetver.h
............此處省略0個(gè)文件信息