Can you code that in python with using NumPy? and explain step by step.

Q2. A linear equation system is given below. 5x1x2+2x3=14x1+12x26x3=4x1+0.5x2+3.5x3=16 Solve the system with Gauss-Siedel iteration method using x1(0)=0.0000,x2(0)=0.0004,x3(0)= 0.0000 as initial values by using 4 decimal points (continue to iteration until results are same to 4 decimal points). Your codes should print estimated values for each iteration. Instructions - Full solution (step by step) should be provided in your report. - The Python codes for Q2 should be also designed and delivered (hint: use NumPy package). You should upload .py file along with your reports