哭求:fortran语言矩阵求逆的程序答:do k=1,N-1 elmax=dabs(Ab(k,k))id_max=k do i=k+1,n if (dabs(Ab(i,k))>elmax) then elmax=Ab(i,k)id_max=i end if end do vtemp1=Ab(k,:)vtemp2=Ab(id_max,:)Ab(k,:)=vtemp2 Ab(id_max,:)=vtemp1 !do i=k+1,N temp=Ab(i,k)/Ab(k,k)Ab(i,:)=Ab...
怎么把一个多项式分解成两个因式的积?答:!求解au=b,u !n表示为方程维数 subroutine lu(a,b,n,u)implicit real(8) (a-h,o-z)real(8)::a(n,n),b(n),u(n),a_bak(n,n),b1(n),aL(n,n),aU(n,n),y(n)!exchange rows do i=1,n tmpMax=0.d0 do ic=i,n if(tmpMax<dabs(a(ic,i))) then tmpMax=dabs(a(...
如何用matlab三角矩阵分解?答:下面是LU分解的Fortran子程序 希望可以有所帮助 !求解au=b,u !n表示为方程维数 subroutine lu(a,b,n,u)implicit real(8) (a-h,o-z)real(8)::a(n,n),b(n),u(n),a_bak(n,n),b1(n),aL(n,n),aU(n,n),y(n)!exchange rows do i=1,n tmpMax=0.d0 do ic=i,n if(tmp...
如何将矩阵A分解成两个矩阵相乘的形式?答:!求解au=b,u !n表示为方程维数 subroutine lu(a,b,n,u)implicit real(8) (a-h,o-z)real(8)::a(n,n),b(n),u(n),a_bak(n,n),b1(n),aL(n,n),aU(n,n),y(n)!exchange rows do i=1,n tmpMax=0.d0 do ic=i,n if(tmpMax<dabs(a(ic,i))) then tmpMax=dabs(a(...
如何将矩阵A分解成一个矩阵和一个向量?答:下面是LU分解的Fortran子程序 希望可以有所帮助 !求解au=b,u !n表示为方程维数 subroutine lu(a,b,n,u)implicit real(8) (a-h,o-z)real(8)::a(n,n),b(n),u(n),a_bak(n,n),b1(n),aL(n,n),aU(n,n),y(n)!exchange rows do i=1,n tmpMax=0.d0 do ic=i,n if(tmp...