题目内容

已知三元一次方程组

请画出解三元一次联列方程组的框图和用Scilab语言编写程序.

答案:略
解析:

解三元一次联列方程组的框图:

解三元一次联立方程组的程序.

write(%io(2),“please input the data at this format!”);

write(%io(2),“——>a11 a12 a13”);

write(%io(2),“——>a21 a22 a23”);

write(%io(2),“——>a31 a32 a33”);

write(%io(2),“——>b1 b2 b3”);

[a]=read(%io(1),3,3);

[b]=read(%io(1),1,3);

for k=1∶2

  t=abs(a,(k,k));

  m=k;

for i=k+1∶3

  if abs(a(i,k))>t then t=abs(a(i,k));

   m=i;end;

end;

if t<le-6 then print(%io(2),“no solution or ununique solution”,t);abort;end;

if m<>k

for i=k∶3

  c=a(m,i);a(m,i)=a(k,i);a(k,i)=c;

end

c=b(m);b(m)=b(k);b(k)=c;

end

for i=k+1∶3

  c=-a(i,k)a(k,k);

  for j=(k+1)∶3

   a(i,j)=a(i,j)+c*a(k,j);

  end

  b(i)=b(i)-c*b(k);

 end

end

if a(3,3)< le-6 then print(%io(2),“no solution or ununique solution”);

abort;end;

x(3)=b(3)/a(3,3);

for k=2∶?-1∶1

  s=0;

  for i=(k+1)∶3

    s=s+a(k,i)*x(i);

  end

x(k)=(b(k)-s)/a(k,k);

end

print(%io(2),x);


练习册系列答案
相关题目

违法和不良信息举报电话:027-86699610 举报邮箱:58377363@163.com

精英家教网