用高斯消去法求方程组解.m

来自「本程序是用MATLAB编写」· M 代码 · 共 41 行

function gjy12
clc
clear
 a=[0.012 0.01 0.167 0.6781;1 0.8334 5.91 12.1;3200 1200 4.2 981]
 %a=[2 -3 1 -3;-3 5 -7 4;1 -2 3 -4]
m=3;n=4;
   for j=1:1:m
        if a(j,j)==0
         
         warndlg('Error')
       
        end 
     for i=j+1:1:m
      g=-a(i,j)/a(j,j)
	  for k=j:1:n
          
	      a(i,k)=g*a(j,k)+a(i,k)
       
      end 
      
  end
     
end
%-----------------------------------------------------------
for j=m:-1:1
    for i=j-1:-1:1
      l=-a(i,j)/a(j,j);
       for k=j:1:n
             
           a(i,k)=l*a(j,k)+a(i,k)
       end
   end
end
%-------------------------------------
disp('kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk')
    for i=1:1:m
        t=a(i,i);
        for j=i:1:n
            a(i,j)=a(i,j)/t
         end
    end