//?example1_2_1.cpp?:?定義控制臺應用程序的入口點。
//

#include?“stdafx.h“
#include
#include
#include

void?main（?）
{
??double?*x*a*b*c*d*rhs*ans*u;
??double?ABhPialphabetalambdamu;
??int?ijN;
??double?*?f（double?*?x?int?N）;
??double?q（double?x）;
??double?*chase_algorithm（double?*a?double?*b?double?*c?int?n?double?*d）;
??m=640;
??Pi=3.14159265359;
??A=0;
??B=Pi;
??h=（B-A）/m;
??alpha=1.0;
??beta=-exp（Pi）;
??lambda=mu=0;
?
??x=（double?*）malloc（sizeof（double）*（m+1））;
??for（i=0;i<=m;i++）
????x[i]=A+i*h;
?
??rhs=f（xm+1）;
??for（i=1;i	??rhs[i]=rhs[i]*h*h;
??rhs[0]=h*alpha;
??rhs[m]=-h*beta;

??a=（double?*）malloc（sizeof（double）*（m+1））;
??b=（double?*）malloc（sizeof（double）*（m+1））;
??c=（double?*）malloc（sizeof（double）*（m+1））;
??d=（double?*）malloc（sizeof（double）*（m+1））;
??for（i=0;i<=m;i++）
??{
	??d[i]=rhs[i];
	??a[i]=1.0;
	??b[i]=h*h*q（x[i]）-2;
	??c[i]=a[i];
??}
??b[0]=lambda*h-1.0;
??b[m]=-（mu*h+1.0）;
??free（rhs）;
??ans=（double?*）malloc（sizeof（double）*（m+1））;
??ans=chase_algorithm（abcm+1d）;
??free（a）;
??free（b）;
??free（c）;
??free（d）;
??
??u=（double?*）malloc（sizeof（double）*（m+1））;???//u為數值解
??for（i=0;i<=m;i++）
	??u[i]=ans[i];
??free（ans）;

??i=m/5;
??for（j=1;j<=4;j++）
	?printf（“x=%f==%d*Pi/5?u=%f\n“A+j*i*hju[i*j]）;
??free（u）;
}


double?*?f（double??*?x?int?N）
{
	int?i;
	double?*ans;
	ans=（double?*）malloc（sizeof（double）*（N+1））;
	for（i=0;i<=N;i++）
		ans[i]=exp（x[i]）*（2*cos（x[i]）-sin（x[i]））;
	return?ans;
}

double?q（double?x）
{
	double?z;
	z=-1.0;
	return?z;
}

double?*?chase_algorithm（double?*a?double?*b?double?*c?int?n?double?*d）
{
??double?*?ans*g*wp;
??int?i;
??ans=（double?*）malloc（sizeof（double）*n）;
??g=（double?*）malloc（sizeof（double）*n）;
??w=（double?*）malloc（sizeof（double）*n）;
??g[0]=d[0]/b[0];
??w[0]=c[0]/b[0];
??
??for（i=1;i??{
	??p=b[i]-a[i]*w[i-1];
	??g[i]=（d[i]-a[i]*g[i-1]）/p;
	??w[i]=c[i]/p;
??}
??ans[n-1]=g[n-1];
??i=n-2;
???do
	?{
		ans[i]=g[i]-w[i]*ans[i+1];
		i=i-1;
		//printf（“i=%d?Now?it‘s?OK!\n“i）;
	?}
??while（i>=0）;
??free（g）;
??free（w）;
??return?ans;
}

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????173437??2014-12-11?13:54??微分方程數值解法的課件\第10講（拋物型方程導數邊界條件）\拋物型方程導數邊界.pptx

?????文件?????615179??2014-12-18?11:21??微分方程數值解法的課件\第11講（對流方程若干解法）\對流方程差分法.pptx

?????文件?????279817??2014-12-21?21:30??微分方程數值解法的課件\第12講（二階雙曲顯差分、隱差分）\二階雙曲方程顯、隱差分法.pptx

?????文件??????10916??2014-11-25?21:08??微分方程數值解法的課件\第12講（二階雙曲顯差分、隱差分）\追趕法子程序.docx

?????文件?????326182??2014-12-25?20:20??微分方程數值解法的課件\第13講（橢圓型方程五點差分格式）\橢圓方程五點差分.pptx

?????文件????????811??2010-01-19?15:27??微分方程數值解法的課件\第1講（歐拉法）\example1.txt

?????文件???????1277??2010-03-10?10:47??微分方程數值解法的課件\第1講（歐拉法）\example2.txt

?????文件???????1003??2010-03-23?21:23??微分方程數值解法的課件\第1講（歐拉法）\Exercises\改進的歐拉方法.txt

?????文件???????1286??2010-03-07?22:27??微分方程數值解法的課件\第1講（歐拉法）\Exercises\梯形方法.txt

?????文件????????787??2010-02-17?14:05??微分方程數值解法的課件\第1講（歐拉法）\Exercises\歐拉方法.txt

?????文件?????487524??2014-10-23?09:04??微分方程數值解法的課件\第1講（歐拉法）\section1.pptx

?????文件????????968??2010-03-09?21:41??微分方程數值解法的課件\第2講（歐拉法續、局部截斷誤差相容性等）\example3.txt

?????文件????????983??2010-03-09?21:50??微分方程數值解法的課件\第2講（歐拉法續、局部截斷誤差相容性等）\example4.txt

?????文件?????402791??2014-10-23?11:04??微分方程數值解法的課件\第2講（歐拉法續、局部截斷誤差相容性等）\section2.pptx

?????文件????????963??2010-02-19?23:02??微分方程數值解法的課件\第3講（龍格－庫塔方法）\example2.txt

?????文件????????966??2010-02-22?14:13??微分方程數值解法的課件\第3講（龍格－庫塔方法）\Exercises\problem1.txt

?????文件?????386592??2014-10-23?11:32??微分方程數值解法的課件\第3講（龍格－庫塔方法）\section3.pptx

?????文件????????947??2010-03-14?17:37??微分方程數值解法的課件\第3講（龍格－庫塔方法）\流程圖解例2.txt

?????文件???????1452??2010-03-22?11:25??微分方程數值解法的課件\第4講（線性多步法）\example3_顯格式.txt

?????文件???????1051??2010-03-22?11:05??微分方程數值解法的課件\第4講（線性多步法）\example3_隱格式.txt

?????文件???????1448??2010-03-25?23:42??微分方程數值解法的課件\第4講（線性多步法）\example4四階adams顯（初值龍四提供）.txt

?????文件???????1379??2010-03-23?16:23??微分方程數值解法的課件\第4講（線性多步法）\example4預估-校正法（初值龍四提供）.txt

?????文件?????855040??2011-03-03?23:53??微分方程數值解法的課件\第4講（線性多步法）\section?4.ppt

?????文件???????1729??2010-04-06?23:01??微分方程數值解法的課件\第5講（一階方程組及高階方程初值問題）\example2.txt

?????文件?????229099??2014-11-13?08:35??微分方程數值解法的課件\第5講（一階方程組及高階方程初值問題）\section5.pptx

?????文件???????1877??2014-11-16?16:03??微分方程數值解法的課件\第6講（常微兩點邊值問題）\example1.txt

?????文件???????2195??2014-11-15?10:27??微分方程數值解法的課件\第6講（常微兩點邊值問題）\example2.cpp

?????文件???????2241??2014-11-15?11:03??微分方程數值解法的課件\第6講（常微兩點邊值問題）\example3.cpp

?????文件????1148098??2014-11-16?16:01??微分方程數值解法的課件\第6講（常微兩點邊值問題）\section6.pptx

?????文件???????1960??2014-11-15?18:17??微分方程數值解法的課件\第7講（外推及緊差分）\Dirichlet.txt

............此處省略25個文件信息