三对角方程直接解法--追赶法+TDMA程序

TnnnnnnfffbacbacbacbAfAx,,,111122211voidZG(intn){for(i=1;i=n-1;i++)U[i]=c[i];L[1]=b[1];for(i=2;i=n;i++){M[i]=a[i]/L[i-1];L[i]=b[i]-M[i]*U[i-1];}Y[1]=f[1];for(i=2;i=n;i++)Y[i]=f[i]-M[i]*Y[i-1];X[n]=Y[n]/L[n];for(i=n-1;i=1;i--)X[i]=(Y[i]-U[i]*X[i+1])/L[i];}0;0;111iiiiiiiiiBniCiDTCTBTAvoidTDMA(intn){P[1]=B[1]/A[1];Q[1]=D[1]/A[1];for(i=2;i=n;i++){P[i]=B[i]/(A[i]-C[i]*P[i-1]);Q[i]=(D[i]+C[i]*Q[i-1])/(A[i]-C[i]*P[i-1]);}T[n]=Q[n];for(i=n-1;i=1;i--)T[i]=P[i]*T[i+1]+Q[i];}