您所在位置：主题：[原创]请各位高手帮看看这个程序什么地方错了
艾瑞达a
[专家分：0] 发布于 2011-10-20 12:22:00
此程序是求线性方程组 ，能运行，就是结果不正确
program ex123
implicit none
integer i,j,k,A(3,3),B(3),x(3)
read(*,*) ((A(i,j),i=1,3),j=1,3)
read(*,*) (B(i),i=1,3)
Do 30,k = 1,2
Do 20,i = k+1,3
Do 10,j = k,3
A(i,j)=A(i,j)-A(k,j)*(A(i,k)/A(k,k))
B(i)=B(i)-B(k)*(A(i,k)/A(k,k))
10 continue
20 continue
30 continue
x(3)=B(3)/A(3,3)
x(2)=(B(2)-A(2,3)*x(3))/A(2,2)
x(1)=(B(1)-A(1,2)*x(2)-A(1,3)*x(3))/A(1,1)
Do 40,i = 1,3
print*,x(i)
40 continue
end
program ex123
implicit none
integer i,j,k,A(3,3),B(3),x(3)
read(*,*) ((A(i,j),i=1,3),j=1,3)
read(*,*) (B(i),i=1,3)
Do 30,k = 1,2
Do 20,i = k+1,3
Do 10,j = k,3
A(i,j)=A(i,j)-A(k,j)*(A(i,k)/A(k,k))
B(i)=B(i)-B(k)*(A(i,k)/A(k,k))
10 continue
20 continue
30 continue
x(3)=B(3)/A(3,3)
x(2)=(B(2)-A(2,3)*x(3))/A(2,2)
x(1)=(B(1)-A(1,2)*x(2)-A(1,3)*x(3))/A(1,1)
Do 40,i = 1,3
print*,x(i)
40 continue
end