数值计算方法程序代码

数值计算方法：最小二乘法的曲线拟合、四阶龙格库塔法求解常微分方程、牛顿法求解非线性方程、列主元素消去法求解线性方程组、拉格朗日插值、改进欧拉方法求解常微分方程、变步长梯形积分、LU分解法求解线性方程组

#include
#include

void main（）
{
	float a[20][20]u[20][20]l[20][20]b[20]y[20]x[20]c;
	int ijrnk;
	FILE*fp1=fopen（“输入.txt““r“）;
	fscanf（fp1“%d“&n）;
	for（i=1;i<=n;i++）
	{for（j=1;j<=n;j++）
	fscanf（fp1“%f“&a[i][j]）;}
	for（i=1;i<=n;i++）
	fscanf（fp1“%f“&b[i]）;
	fclose（fp1）;




    for（j=1;j<=n;j++）
    {u[1][j]=a[1][j];}
	for（i=1;i<=n;i++）
	l[i][i]=1;
    for（i=2;i<=n;i++）
    l[i][1]=a[i][1]/u[1][1];



	    for（r=2;r<=n;r++）
		{  
			for（i=r;i<=n;i++）
			{ c=0;
				for（k=1;k<=r-1;k++）
			{c=c+l[r][k]*u[k][i];
				
				u[r][i]=a[r][i]-c;}}
			for（i=r+1;i<=n;i++）
			{ c=0;
				for（k=1;k<=r-1;k++）
				{c=c+l[i][k]*u[k][r];}
			l[i][r]=（a[i][r]-c）/u[r][r];}
		}


 
          y[1]=b[1];
		for（i=2;i<=n;i++）
		{	c=0;
 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2016-10-24 23:54  数值计算方法201501000414马子儒\
     目录           0  2016-10-24 23:32  数值计算方法201501000414马子儒\LU分解法求解线性方程组\
     文件        1167  2016-10-24 23:31  数值计算方法201501000414马子儒\LU分解法求解线性方程组\LU分解法求解线性方程组.cpp
     文件          32  2016-10-24 23:29  数值计算方法201501000414马子儒\LU分解法求解线性方程组\输入.txt
     文件          31  2016-10-24 23:31  数值计算方法201501000414马子儒\LU分解法求解线性方程组\输出.txt
     目录           0  2016-10-24 23:28  数值计算方法201501000414马子儒\列主元素消去法求解线性方程组\
     文件         876  2016-10-24 23:27  数值计算方法201501000414马子儒\列主元素消去法求解线性方程组\列主元素消去法求解线性方程组.cpp
     文件          36  2016-10-24 23:25  数值计算方法201501000414马子儒\列主元素消去法求解线性方程组\输入.txt
     文件          30  2016-10-24 23:27  数值计算方法201501000414马子儒\列主元素消去法求解线性方程组\输出.txt
     目录           0  2016-10-24 23:56  数值计算方法201501000414马子儒\变步长梯形积分\
     文件         721  2016-10-19 17:27  数值计算方法201501000414马子儒\变步长梯形积分\变步长梯形求积法.cpp
     文件          17  2016-10-19 17:19  数值计算方法201501000414马子儒\变步长梯形积分\输入.txt
     文件           9  2016-10-19 17:27  数值计算方法201501000414马子儒\变步长梯形积分\输出.txt
     目录           0  2016-10-24 23:44  数值计算方法201501000414马子儒\四阶龙格－库塔法求解常微分方程的初值问题\
     文件         602  2016-10-18 19:56  数值计算方法201501000414马子儒\四阶龙格－库塔法求解常微分方程的初值问题\四阶龙格－库塔法求解常微分方程的初值问题.cpp
     文件          13  2016-10-18 19:52  数值计算方法201501000414马子儒\四阶龙格－库塔法求解常微分方程的初值问题\输入.txt
     文件         135  2016-10-18 19:56  数值计算方法201501000414马子儒\四阶龙格－库塔法求解常微分方程的初值问题\输出.txt
     目录           0  2016-10-24 23:41  数值计算方法201501000414马子儒\拉格朗日插值\
     文件         666  2016-10-24 23:41  数值计算方法201501000414马子儒\拉格朗日插值\拉格朗日插值.cpp
     文件          23  2016-10-24 23:40  数值计算方法201501000414马子儒\拉格朗日插值\输入1.txt
     文件           2  2016-10-24 23:40  数值计算方法201501000414马子儒\拉格朗日插值\输入2.txt
     文件          10  2016-10-24 23:41  数值计算方法201501000414马子儒\拉格朗日插值\输出.txt
     目录           0  2016-10-24 23:44  数值计算方法201501000414马子儒\改进欧拉方法求解常微分方程的初值问题\
     文件         505  2016-10-18 19:58  数值计算方法201501000414马子儒\改进欧拉方法求解常微分方程的初值问题\改进欧拉方法求解常微分方程的初值问题.cpp
     文件          14  2016-10-18 19:18  数值计算方法201501000414马子儒\改进欧拉方法求解常微分方程的初值问题\输入.txt
     文件          81  2016-10-18 19:58  数值计算方法201501000414马子儒\改进欧拉方法求解常微分方程的初值问题\输出.txt
     目录           0  2016-10-24 23:42  数值计算方法201501000414马子儒\最小二乘法的曲线拟合\
     文件         974  2016-10-18 20:50  数值计算方法201501000414马子儒\最小二乘法的曲线拟合\最小二乘法的曲线拟合.cpp
     文件          34  2016-10-18 20:40  数值计算方法201501000414马子儒\最小二乘法的曲线拟合\输入.txt
     文件          22  2016-10-18 20:50  数值计算方法201501000414马子儒\最小二乘法的曲线拟合\输出.txt
     目录           0  2016-10-24 23:24  数值计算方法201501000414马子儒\牛顿法求解非线性方程\
............此处省略3个文件信息