小唐？
typedef CArray<double,double>CDoubleArray;
void ZG(CDoubleArray *A,CDoubleArray *B,CDoubleArray *C,CDoubleArray *G,int &LOGI)
{
//追赶法
register long I;
int N;
N=A->GetSize();
if(LOGI==0)
{
(*C)[0]=(*C)[0]/(*B)[0];
for(I=1;I<N;I++)//DO 2 I=2,N
{
(*B)=(*B)-(*A)*(*C)[I-1];
(*C)=(*C)/(*B);
}
(*A)[0]=0.;
(*C)[N-1]=0.;
LOGI=1;
}
(*G)[0]=(*G)[0]/(*B)[0];
for(I=1;I<N;I++)//DO 20 I=2,N
{
(*G)=((*G)-(*A)*(*G)[I-1])/(*B);
}
for(I=N-2;I>-1;I--)//DO 30 I=N-1,1,-1
{
(*G)=(*G)-(*C)*(*G)[I+1];
}
return;
}
