zhouwen.prg

202用列主元消取法解线性方程 ***********★*******★********★************ 一．功能 当线性方程组有唯一解时求其解 。 二．算法简介消元过程

function  zhouwenliang
para n
Public b(10,10)
 n==4

set proc to myproc
for x=1 to n
s=0
for y=1 to m
*set proc to myproc
n=val(thisform.qn1.text)
m=val(thisform.qm1.text)
ph=alltrim(thisform.in1.value)
 for i=1 to n
    for j=1 to m
    s=at(",",ph)
    a(i,j)=val(subs(ph,1,s-1)) 
     if(j==m)
      ph=alltrim(subs(ph,s+3))
     else
      ph=alltrim(subs(ph,s+1))     
    endif
   endfor
 endfor
* ? a(1,2)   
a(1,y)=0
a(2,y)=0
a(3,y)=0
a(4,y)=0
for j1=1 to n
for j2=1 to n
for j3=1 to n
if(j1#j2.and.j1#j3.and.j2#j3)
if (x==1)
s=s+qupi(j1,j2,j3)*a(2,j1)*a(3,j2)*a(4,j3)
endif
if (x==2)
s=s+qupi(j1,j2,j3)*a(1,j1)*a(3,j2)*a(4,j3)
endif
if (x==3)
s=s+qupi(j1,j2,j3)*a(1,j1)*a(2,j2)*a(4,j3)
endif
if (x==4)
s=s+qupi(j1,j2,j3)*a(1,j1)*a(2,j2)*a(3,j3)
endif
endif
endfor
endfor
endfor
if(mod((x+y),2)==0)
b(x,y)=(s-2)
endif
if (mod((x+y),2)==1)
b(x,y)=(-1)*(s-2)
endif
thisform.edit1.value=""
for i=1 to n
p=""
for j=1 to n
p=p+str(d(i,j),n)
endfor
thisform.edit1.value=thisform.edit1.value+p+chr(13)
endfor