#include<iostream>
#include<iomanip>
using namespace std;
const float h=0.01; //步长
const float Tmax=10;//最大时间
void gauss(float arr[][3],int n);
int main(){
    //变量说明
    float I3[1001],U3[1001],I5[1001],U5[1001];
    float YN[2][3];//结点方程的增广阵
    float t; //时间变量
    int m=2;//阶数
    float Gl,Gc;
    int k;
    //输入电路参数
    float Is=1,R2=20,R4=2,L=2,C=0.2;
    //计算Gl,Gc
    Gl=h/(2*L);
    Gc=2*C/h;
    //初始值计算
    I3[0]=0;U5[0]=1;
    I5[0]=(-1-R2)/(R2+R4);
    U3[0]=(1-I5[0])*R2;
    t=0.01;k=0;
    while(t<=Tmax){
        //增广阵
        YN[0][0]=1/R2+1/R4+Gl;
             YN[0][1]=-1/R4;
        YN[1][0]=-1/R4;
        YN[1][1]=1/R4+Gc;
       
        YN[0][2]=Is-I3[k]-Gl*U3[k];
        YN[1][2]=I5[k]+Gc*U5[k];
        //调用gauss函数，解方程
                  gauss(YN,m);
        U3[k+1]=YN[0][2];
        U5[k+1]=YN[1][2];
        //计算下一次初值
        I3[k+1]=Gl*U3[k+1]+I3[k]+Gl*U3[k];
        I5[k+1]=Gc*(U5[k+1]-U5[k])-I5[k];
        //时间变量增加一个步长，k加1
        t=t+h;
        k++;
    }
    //每隔20个时间点输出
    cout<<"T(s)"<<setw(15)<<"IL3(A)"<<setw(15)<<"UC5(V)"<<endl;
    for(k=0;k<=1000;k++){
        if(k%20==0){
            cout<<k*h<<setw(15)<<I3[k]<<setw(15)<<U5[k]<<endl;
        }
    }
    return 0;
}
//gaus function
void gauss(float arr[][3],int n){
    int m;
    m=n+1;
    int i,j,k;
    
    //消元
    float a;
    for(k=0;k<=n-2;k++){
        for(i=k+1;i<=n-1;i++){
          a=arr[i][k]/arr[k][k];
            for(j=0;j<=n;j++){
               arr[i][j]=arr[i][j]-a*arr[k][j];
            }
        }
    }
    //回代求解
    arr[n-1][n]=arr[n-1][n]/arr[n-1][n-1];
    for(j=n-1;j>=1;j--){
        for(i=j-1;i>=0;i--){
        arr[i][n]=arr[i][n]-arr[i][j]*arr[j][n];
        }
        arr[j-1][n]=arr[j-1][n]/arr[j-1][j-1];
    }
}
这是一个解电路计算的程序，计算最后一次值U3[1000],U5[1000]时，结果不对，不知道怎么回事，请大侠帮忙看看