工程电磁场迭代法程序

#include

Void main ()

{

Double a[5][5],t,d,b;

Int i,j,n=100000;

For (i=0;i<5;i++)

{

For (j=0;j<5;j++)

If(i==0)

A[i][j]=100;

Else

A[i][j]=0;

}

For (i=0;i<5;i++)

{

For (i=0;i<5;i++)

Cout<Cout<}

For (d=1;i<2;d+=0.001)

{

Int flag =0;

Double p[5][5];

For (i=0;i<5;i++)

For (j=0;j<5;j++)

P[i][j]=a[i][j];

Int k=0;

While(flag<9)

{

For (i=1;i<4;i++)

For (j=1;j<4;j++)

T=p[i][j];

P[i][j]=p[i][j]+d*(p[i+1][j]+p[i][j+1]p[i-1][j]+p[i][j-1]-4*p[i][j]/4);

If (fabs(p[i][j]-t)<1e-5)

Flag++;

}

K++;

}

If(n>=k)

{

N=k;

B=d;

}

}

}

Cout<<\"最佳收敛因子:\"<Cout<<\"各点电压最终值:\"<Double q[5][5];

For(i=0;i<5;i++)

For(j=0;j<5;j++)

Q[i][j]=a[i][j];

While(n)

{

For(i=0;i<5;i++)

For(j=0;j<5;j++)

Q[i][j]=q[i][j]+b*(q[i+1][j]+q[i][j+1]+q[i-1][j]+q[i][j-1]-4*q[i][j])/4;

N--;

}

For(i=0;i<5;i++)

{

For(j=0;j<5;j++)

Cout<Cout<}

}