C语言实现Crout分解

最新推荐文章于 2024-08-23 07:32:51 发布

平林寒山98

最新推荐文章于 2024-08-23 07:32:51 发布

阅读量1.2k

点赞数

分类专栏：数值分析

本文链接：https://blog.youkuaiyun.com/qq_40755405/article/details/106205667

版权

在这里插入图片描述
PS：3、4公式要同时进行，第4公式的j要从k开始。

#include<stdio.h>
#define N 4
#define col 4//列数
#define row 4//行数
void print(double a[row][col]) {
   

	for (int i = 0; i < row; i++)//找到最大主元，记录行号和列号
	{
   
		for (int j = 0; j < col; j++)
		{
   
			printf("%f  ", a[i][j]);
		}
		printf("\n");
	}
	printf("\n");
}

void printXY(double XY[N]) {
   

	for (int i = 0; i < N; i++)
	{
   
		printf("%f  ", XY[i]);
		
	}
	printf("\n\n");
}

void step1(double a[row][col], double L[row][col], double  U[row][col])
{
   
	for (int i = 1; i <= N; i++)
	{
   
		L[i - 1][0] = a