This question is to develop a program that solves systems of linear equations. Assume the number of equations is always the same as the number of unknowns. The general form of such a system is as follows:

a11x1 + a12x2 + … + a1nxn = b1,
a21x1 + a22x2 + … + a2nxn = b2,

an1x1 + an2x2 + … + annxn = bn.

Your program first asks the user to enter the number n of unknowns. Then ask the user to enter the coefficients and the constant term of each equation, one equation in one line. For instance, if n is 3, the user may enter the following:

120 -32.1 20.4 50
1.6 -3.5 2.7 1
-0.2 0.15 0.09 0

This input represents the following system of equations

120x1 -32.1x2 + 20.4x3 = 50,
1.6 x1 -3.5 x2 + 2.7 x3 = 1,
-0.2x1 + 0.15x2 + 0.09 x3 = 0.

A 2-dimensional array is built dynamically to store the input data (so-called the augmented matrix of the system of equations). Then use an algorithm, called scaled Gauss elimination with backward substitution, to solve this system of equations. The algorithm is as follows:

First, divide each of the equations by the coefficient that has the largest absolute value (not including the constant term) to convert this system to a so-called normalized system of equations. Then use the Gauss elimination algorithm:

Suppose the normalized system is as follows:

a11x1 + a12x2 + … + a1nxn = b1,
a21x1 + a22x2 + … + a2nxn = b2,

an1x1 + an2x2 + … + annxn = bn.

Find the coefficient in the first column (a11, a21, …, an1) with the largest absolute value. If ai1 is the largest, swap equation 1 with equation i. Then, we may assume a11 has the largest absolute value. If this value is zero, this system does not have a unique solution. Output an error message and terminate the program. Otherwise, denote the i-th equation by Ei. Replace Ei by Ei - E1 * (ai1 / a11), i = 2, 3, …, n. After that, in the first column, all coefficients under a11 become 0. The system has the following form:

a11x1 + a12x2 + … + a1nxn = b1,
a22x2 + … + a2nxn = b2,

an2x2 + … + annxn = bn.

Next, find the coefficient among a22, a32, …, an2 that has the largest absolute value. If aj2 has the largest absolute value, interchange E2 and Ej. Then we may assume a22 has the largest absolute value. If this value is zero, this system does not have a unique solution. Output an error message and terminate the program. Otherwise, Replace Ei, i = 3, 4, …, n, by Ei - E2 * (ai2 / a22). Then all coefficients under a22 become zero.

Keep doing this. Finally, we have a system that has the so-called echelon form:

a11x1 + a12x2 + a13x3 + … + a1nxn = b1,
a22x2 + a23x3 + … + a2nxn = b2,

annxn = bn,

where aii are not zero, i = 1, 2, …, n.

In the last step, find the values backwards: Form the last equation, we have

xn = bn / ann.

Substituting this value into the second last equation, we have

xn-1 = (bn-1- an-1, nxn-1) / an-1, n-1.

In this way, we can find xn-2, xn-3, … . Finally,

x1 = (b1 - a12x2 - a13x3 -- a1nxn) / a11.

To make your calculation more accurate, all coefficients should be of the type long double. In this algorithm, we need to check whether a coefficient is zero. However, it is not good to compare a double value with zero because a real number may not be expressed by a binary number accurately. You should consider a value being zero if its absolute value is less than, say, 10-10.

At the end of the program, don't forget to free the memory used by the dynamic array.

2. This question is to review how to manipulate c-strings and linked lists. (Do not use C++ class string in this question).

You know that there is a method called tokenize in Java, which separates a string into parts according to a set of specified characters called delimiters. Unfortunately, there is no such built-in function in C++. This question is to write such a function tokenize( ). This function accepts two strings. The first string is to be tokenized, and the second string contains all delimiters. For instance, if the first string is "sin(x)+cos(y)-tan(z)", and the second string is "+_()", then this function should separate the first string into the following tokens:

"sin", "(", "x", ")", "+", "cos", "(", "y", ")", "-", "tan", "(", "z", ")"

Each delimiter is also regarded as a token.

Note that, if you have two delimiters side by side (as ")+" and ")-" in the previous example, you should not have an empty token between them.

Your function should return a linked list that stores the tokens of the input string in the same order as they are in the input string. The node of the linked list has the type

struct Node
{
char * token;