///修改方程系数
void CTouchDlg::calculate()
{
//UpdateData(TRUE);
int a2=atoi(cs6)；//左上
int b2=atoi(cs7);
int c2=atoi(cs10)；//左下
int d2=atoi(cs11);
int a1=atoi(cs8)；//右上
int b1=atoi(cs9);
int c1=atoi(cs12)；//右下
int d1=atoi(cs13);
int e=1；//第一个系数
//开始
int hang=三、//初始化行列
int lie=4;
long double a[100][101]
long double b[100][101];
//初始化
for(int i=0;i<100;i )
{
for(int j=0;j<101;j )
{
a[i][j]=0;
}
}
a[0][0]=1; //第1组
a[0][1]=a2;
a[0][2]=b2;
a[0][3]=0;
a[1][0]=1; //第2组
a[1][1]=a1;
a[1][2]=b1;
a[1][3]=e0;
a[2][0]=1; //第3组
a[2][1]=c2;
a[2][2]=d2;
a[2][3]=0;
int z=3;
//jisuan
{
for(int k3=0;k3<100;k3 )
{
for(int h=0;h<101;h )
{
b[k3][h]=a[k3][h];
}
}
for(int zi=0;zi
{
//换列(zi 1);i=zi 1!!!111
//...换列...
int jilu=-1;
for(int v=0;v
{
if(b[zi v][zi]!=0)//changge!!!
{
jilu=zi 1 v;
break;
}
}
if(jilu==-1)
{
}
else
{
long double kkk=1;
//exchange(i,jilu);
//int j=jilu;///记录
//...交换...11111111111111111111
{
long double temp=0;
for(int t3=0;t3
{
temp=b[zi][t3];
b[zi][t3]=b[jilu-1][t3];
b[jilu-1][t3]=temp;
}
}
//...end 交换...11111111111111111111
//chengk(i,(1.0/b[i-1][i-1]));
//...检查...11111111111111
{
long double k=(1.0/b[zi][zi]);/注意！！
for(int c=0;c
{
b[zi][c]=b[zi][c]*k;
}
}
//... end 检查...
{
for(int w=0;w
{
if(b[zi][zi]!=0)
{
kkk=0.0-(b[zi 1 w][zi]/b[zi][zi]);
//chengplus(zi 1,kkk,zi 1 1 w);
//...相乘...
for(int t=0;t
{
b[zi 1 w][t] =b[zi][t]*kkk;
}
//... end 相乘...
}
}
}//555555
}
//...end 换列...11111111111111111
}
for(int zii=0;zii
{
//Uhualie(z-zii);i=z-zii;!!!
//...22222222222222222
if(b[z-zii-1][z-zii-1]!=0)
{
long double kkk1=1;
for(int w=0;w
{
if(b[z-zii-1][z-zii-1]!=0)
{
kkk1=0.0-(b[z-zii-2-w][z-zii-1]/b[z-zii-1][z-zii-1]);
//chengplus(i,kkk1,z-zii-1-w);
//...55
for(int t=0;t
{
b[z-zii-1-w-1][t] =b[z-zii-1][t]*kkk1;
}
//... end 55 ...
}
}
}
//...end 5...
}//end 2222222222222222222
//... 6666...3333333333333333
for(int zck=0;zck
{
if(b[zck][zck]!=0)
{
//chengk(zck 1,(1.0/b[zck][zck]));
long double k2=(1.0/b[zck][zck]);
//... ...
for(int c1=0;c1
{
b[zck][c1]=b[zck][c1]*k2;
}
//... end ......................................
}
}
//...........................end .......................................333333333333
}//end jisuan
{
//........................... 打印..................................4444444444444444
int testa=0;
int testb=0;
for(int m=0;m
{
for(int n=0;n
{
if(b[m][n]!=0){testb+=1;}
}
if(testb==0&&b[m][lie-1]!=0){testa+=1;}
testb=0;
}
if(testa!=0)
{
//wujie
}
else
{
//UpdateData(TRUE);
//cout <
//m_edt3.Format("%d",atoi(m_edt2) + atoi(m_edt1));xianshi
long double d1=b[0][3];
long double d2=b[1][3];
long double d3=b[2][3];
long double d4=b[2][3];
cs0.Format("%lf",d1);
cs1.Format("%lf",d2);
cs2.Format("%lf",d3);
UpdateData(FALSE);
}
//...........................end printb.................................4444444444444444
}
}