小编Ans*_*nha的帖子

我编写了一个1 Dimmension cfd问题,但我的数值解决方案与分析解决方案相同(最多6位小数).

我使用TDMA方法进行数值求解,对于解析解,我直接用函数T(x)中的x值代替.

分析解T(x)出来了 T(x) = -(x^2)/2 +11/21(x);

E. g.然后是4格点;

  x0  = 0.000000,    x1  = 0.333333 ,   x2  = 0.666666 ,   x3  = 0.999999 .
T(x0) = 0.000000 , T(x1) = 0.119048 , T(x2) = 0.126984 , T(x3) = 0.023810.

对于数值解决方案,我使用过TDMA技术,请参考下面的代码.

输入n = 4表示结果.

#include<stdio.h>

void temp_matrix(int n, double *a, double *b, double *c, double *d, double *T);

int main() {
  int Bi = 20.0;
  int n;
  printf("%s ", "Enter the Number of total Grid Points");
  scanf("%d", &n);
  float …