追赶法求解方程组（C语言）

最新推荐文章于 2024-08-18 22:08:24 发布

landcruiser007

最新推荐文章于 2024-08-18 22:08:24 发布

阅读量1.1w

点赞数 9

分类专栏：计算方法文章标签：计算方法追赶法

本文链接：https://blog.csdn.net/landcruiser007/article/details/79303596

版权

计算方法专栏收录该内容

16 篇文章 11 订阅

订阅专栏

编写用追赶法解三对角线方程组的程序，并解下列方程组：

这里写图片描述

#include<stdio.h>
#include<math.h>
void main()
{
    int i;
    int n=4;

    float x[4];
    float c[4][5]={2,-1,0,0,5,
                   -1,2,-1,0,-12,
                   0,-1,2,-1,11,
                   0,0,-1,2,-1};

    float y[10];
    float d[10][11]={-4,1,0,0,0,0,0,0,0,0,-27,
                     1,-4,1,0,0,0,0,0,0,0,-15,
                     0,1,-4,1,0,0,0,0,0,0,-15,
                     0,0,1,-4,1,0,0,0,0,0,-15,
                     0,0,0,1,-4,1,0,0,0,0,-15,
                     0,0,0,0,1,-4,1,0,0,0,-15,
                     0,0,0,0,0,1,-4,1,0,0,-15,
                     0,0,0,0,0,0,1,-4,1,0,-15,
                     0,0,0,0,0,0,0,1,-4,1,-15,
                     0,0,0,0,0,0,0,0,1,-4,-15};

    void zhui(float *,int,float[]);
    zhui(c[0],4,x); 
        zhui(d[0],10,y);    
    for(i=0;i<=3;i++)printf("x[%d]=%f\n",i+1,x[i]);
    printf("\n");
    for(i=0;i<=9;i++)printf("x[%d]=%f\n",i+1,y[i]); 
}
//追赶法
void zhui(float *u,int n,float x[])
{
    int i,j;
    float p;

for(i=1;i<=n-1;i++)
    {
        p=*(u+i*(n+1)+i-1)/(*(u+(i-1)*(n+1)+i-1));
        *(u+i*(n+1)+i-1)=0;
        *(u+i*(n+1)+i)-=p*(*(u+(i-1)*(n+1)+i));
        *(u+i*(n+1)+n)-=p*(*(u+(i-1)*(n+1)+n));
    }

    x[n-1]=(*(u+(n-1)*(n+1)+n))/(*(u+(n-1)*(n+1)+n-1));
    for(j=n-2;j>=0;j--)
    {

        x[j]=((*(u+j*(n+1)+n))-x[j+1]*(*(u+j*(n+1)+j+1)))/(*(u+j*(n+1)+j));
    }
}