www.pudn.com > branch.rar > branch.cpp

#include  
#include  
 
double a[10][10]; 
 
void Guass(int n){ 
 int i,j,k; 
 double t; 
 for(k=0;k0.001) 
      {found=0; 
      break; 
      } 
     found=1; 
     } 
    if(diff>0.001) break; 
    } 
   if(diff>0.0001) break; 
   } 
  if(found==0) 
   {if(k%5==0) 
    {printf("  step %d",k); 
    for(i=0;i=100) 
    {printf(" steady_state probability have not been detained in 100"); 
    return; 
    } 
   } 
  else 
   {printf(" step %d",k); 
   for(i=0;i