
68 Programming for Chemical Engineers
38 printf(“\tValue at right side = “);
39 scanf(“%lf”,&simtab[i+1][1]);
40 }
41 printf(“\n____________________\n\n”);
42 for(j=1; j<=ne; j++) simtab[0][j+1] = j;
43 for(i=ne+1; i<=ne+nc; i++) simtab[i-ne+1][0] = i;
44 do{
45 xmax = 0.0;
46 for(j=2; j<=ne+1; j++)
47 {
48 if (simtab[1][j] > 0.0 && simtab[1][j] > xmax)
49 {
50 xmax = simtab[1][j];
51 pivot2 = j;
52 }
53 }
54 flag = 1000000.0;
55 for (i=2; i<=nc+1; i++)
56 {
57 if (simtab[i][pivot2] < 0.0)
58 {
59 aux = fabs(simtab[i][1] / simtab[i][pivot2]);
60 if (aux < flag)
61 {
62 flag = aux;
63 pivot1 = i;
64 }
65 }
66 }
67 aux = simtab[0][pivot2];
68 simtab[0][pivot2] = ...