（ [微分方程的数值解法与程序实现][华冬英李祥贵][电子课件

网上有文字版图书，这是课后习题以及程序实现，用c写的，包括讲义，共三部分

// 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个文件信息