lescape.c

来自「包括线性插值」· C语言 代码 · 共 25 行

void Lescape(double a[37],double b[13],int n)/*将三对角矩阵的元素采用一微数组存储,n微矩阵，共存3n-2个数据*/
{
  int i,j;
  double p[13],q[12],y[13];
  p[0]=a[0];
  j=1;
  for(i=0;i<=n-2;i++)
  {  	
  	q[i]=a[2*i+j]/p[i];
  	p[i+1]=a[2*(i+1)+j]-a[2*(i+1)+(j-1)]*q[i];
  	j=j+1;  	
  }
  
  y[0]=b[0]/p[0];
  j=1;
  for(i=1;i<=n-1;i++)
  {  	
  	y[i]=(b[i]-a[(2*i)+(j-1)]*y[i-1])/p[i];
  	j=j+1;
  } 
  b[n-1]=y[n-1];
  for(i=n-2;i>=0;i--)
  	b[i]=y[i]-q[i]*b[i+1];
}