《实用数值计算方法》甄西丰 全部源码 C语言

教材《实用数值计算方法》甄西丰著 的全部源码 C语言 含1~13章

#include 
#include 
#include 
#define   MN 51
int       RN=0;
double e = 2.71828;
double   EPS = 0.4E-10;
double   H[MN]DTF[MN]DQF[MN]ERR[MN];
double   X0=2.0H0=0.32;
char     FileName[]=“TB010702.txt“;

int choose = 1;

double   DeltF（double x double h）
{
	return log（（x + h） / x）;
}
double DeltF_ex（double x double h）
{
	return pow（e x + h） - pow（ex）;
}
double DeltF_sin（double x double h）
{
	return sin（x + h） - sin（x）;
}
double DeltF_cos（double x double h）
{
	return cos（x + h） - cos（x）;
}
double DeltF_tan（double x double h）
{
	return tan（x + h） - tan（x）;
}
int Operation（）
{
	int K;
	H[0] = H0;
	if （choose == 1）
		DTF[0] = DeltF（X0 H[0]）;
	if （choose == 2）
		DTF[0] = DeltF_ex（X0 H[0]）;
	if （choose == 3）
		DTF[0] = DeltF_sin（X0 H[0]）;
	if （choose == 4）
		DTF[0] = DeltF_cos（X0 H[0]）;
	if （choose == 5）
		DTF[0] = DeltF_tan（X0 H[0]）;
	DQF[0] = DTF[0] / H[0];
	ERR[0] = 1.0;
	for （K = 1; K < MN; K++）
	{
		H[K] = H[K - 1] / 2.0;
		if （choose == 1）
			DTF[K] = DeltF（X0 H[K]）;
		if （choose == 2）
			DTF[K] = DeltF_ex（X0 H[K]）;
		if （choose == 3）
			DTF[K] = DeltF_sin（X0 H[K]）;
		if （choose == 4）
			DTF[K] = DeltF_cos（X0 H[K]）;
		if （choose == 5）
			DTF[K] = DeltF_tan（X0 H[K]）;
		DQF[K] = DTF[K] / H[K];
		ERR[K] = fabs（DQF[K] - DQF[K - 1]）;
		if （ERR[K] < EPS）break;
	}
	RN = K + 1;
	return 0;
}

int ShowTable（）
{ int K;
  printf（“\n ===================================================“）;
  printf（“\n   k     H[k]      DTF[k]      DQF[k]       ERR[k]  “）;
  printf（“\n ---------------------------------------------------“）;
  for（K=0;K  printf（“\n  %2d %10.6f%12.8f%12.8f%12.4e“KH[K]DTF[K]DQF[K]ERR[K]）;
  printf（“\n ---------------------------------------------------“）;
  printf（“\n  ANS=%12.8e“DQF[RN-1]）;
  printf（“\n ===================================================“）;
  return 0;
}
int SaveTable（char*fname）
{ FILE *fp;
  int K;
  fp=fopen（fname“w“）;
  if（!fp）{printf（“\nOpen File Error !“）;return 1;}
  fprintf（fp“\n===================================================“）;
  fprintf（fp“\n  k     H[k]      DTF[k]      DQF[k]       ERR[k]  “）;
  fprintf（fp“\n---------------------------------------------------“）;
  for（K=0;K  fprintf（fp“\n %2d %10.6f%12.8f%12.8f%12.4e“KH[K]DTF[K]DQF[K]ERR[K]）;
  fprintf（fp“\n---------------------------------------------------“）;
  fprintf（fp“\n ANS = %12.8e“DQF[RN-1]）;
  fprintf（fp“\n===================================================“）;
  fclose（fp）;
  printf（“\nThis Table Has Saved In File %s !“fname）;
  return 0;
}
int FormProblem（）
{
	printf（“选择要计算导数的函数：\n“）;
	printf（“1.对数函数f（x）=lnx\n“）;
	printf（“2.指数函数f（x）=e^x\n“）;
	printf（“3.正弦函数f（x）=sin（x）\n“）;
	printf（“4.余弦函数f（x）=cos（x）\n“）;
	printf（“5.正切函数f（x）=cos（x）\n“）;
	scanf（“%d“ &choose）;
	printf（“\nPlease Enter x0:“）;
	scanf（“%lf“ &X0）;
	printf（“\nPlease Enter h0:“）;
	scanf（“%lf“ &H0）;

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

    ..A..H.     33792  2019-04-06 16:53  Program\Ch01_Intro\.vs\P0101\v15\.suo

     文件    1851392  2019-04-06 16:53  Program\Ch01_Intro\.vs\P0101\v15\Browse.VC.db

     文件    3014656  2019-04-06 14:59  Program\Ch01_Intro\.vs\P0101\v15\ipch\AutoPCH\4bf13984959085e8\P0101.ipch

     文件    3014656  2019-04-06 15:12  Program\Ch01_Intro\.vs\P0101\v15\ipch\AutoPCH\971764dbef1a0cd6\P0201.ipch

     文件    2555904  2019-04-06 15:12  Program\Ch01_Intro\.vs\P0101\v15\ipch\AutoPCH\f1098e7dacd0e129\P0206.ipch

     文件     906752  2019-04-06 15:12  Program\Ch01_Intro\Debug\P0101.exe

     文件    2391748  2019-04-06 15:12  Program\Ch01_Intro\Debug\P0101.ilk

     文件       1453  2019-04-06 15:12  Program\Ch01_Intro\Debug\P0101.log

     文件      18519  2019-04-06 14:48  Program\Ch01_Intro\Debug\P0101.obj

     文件    8572928  2019-04-06 15:12  Program\Ch01_Intro\Debug\P0101.pdb

     文件       2318  2019-04-06 15:12  Program\Ch01_Intro\Debug\P0101.tlog\CL.command.1.tlog

     文件      10246  2019-04-06 15:12  Program\Ch01_Intro\Debug\P0101.tlog\CL.read.1.tlog

     文件       1712  2019-04-06 15:12  Program\Ch01_Intro\Debug\P0101.tlog\CL.write.1.tlog

     文件       2054  2019-04-06 15:12  Program\Ch01_Intro\Debug\P0101.tlog\link.command.1.tlog

     文件       3698  2019-04-06 15:12  Program\Ch01_Intro\Debug\P0101.tlog\link.read.1.tlog

     文件        588  2019-04-06 15:12  Program\Ch01_Intro\Debug\P0101.tlog\link.write.1.tlog

     文件        243  2019-04-06 15:12  Program\Ch01_Intro\Debug\P0101.tlog\P0101.lastbuildstate

     文件       5597  2019-04-06 14:57  Program\Ch01_Intro\Debug\P0201.obj

     文件      10892  2019-04-06 15:12  Program\Ch01_Intro\Debug\P0206.obj

     文件      60416  2019-04-06 15:12  Program\Ch01_Intro\Debug\vc141.idb

     文件      77824  2019-04-06 15:12  Program\Ch01_Intro\Debug\vc141.pdb

     文件      33792  2017-02-27 14:44  Program\Ch01_Intro\Debug\vc60.idb

     文件      53248  2017-02-27 14:44  Program\Ch01_Intro\Debug\vc60.pdb

     文件       3233  2019-04-06 14:48  Program\Ch01_Intro\P0101.c

     文件       3387  2017-02-20 11:28  Program\Ch01_Intro\P0101.dsp

     文件        535  2017-02-20 11:36  Program\Ch01_Intro\P0101.dsw

     文件      41984  2017-02-27 14:48  Program\Ch01_Intro\P0101.ncb

     文件      48640  2017-02-27 14:48  Program\Ch01_Intro\P0101.opt

     文件       1173  2017-02-27 14:44  Program\Ch01_Intro\P0101.plg

     文件       1064  2019-04-06 13:31  Program\Ch01_Intro\P0101.sln

............此处省略756个文件信息