This is the question:

Find the solution [a*, b*, c*, d*] to the nonlinear system of equations: (a) (5 points) 20ac - 8ab + 16c3 - 4bd = 39 12a + 6b + 2c - 2d = 11 4a2 + 2bc - 10c + 2ad- = -7 -3ad - 262 + 7cd = 16. (b) (5 points) 3a - cos(bc) = 0.5 a2 - 81(b + 0.1)2 + sin(c) = -1.06 3 - 107 e-ab + 20c = 3 Use initial iterate x = [1, 1, 1, 1] for (a) and x = [1, 1, 1] for (b), tolerance 10- and maximum iteration of 100. Include in your documentation the solution [a*, b*, c*, d*], function value of the solution, relative error and the number of iteration. Item3Imp. py should contain the code you used to solve the problem using the Newton Method.D 22 J = D np . array ( [[20*c-8*b, -8*a, 20*a-4*d, -4*b], Exception has occurred: ValueError X setting an array element with a sequence. The requested array has an inhomogeneous shape after 1 dimensions. The detected shape was (4, ) + inhomogeneous part. File "C: Item2Imp. py", line 22, in jacobian_a J = np. array( [ [20*c-8*b, -8*a, 20*a-4*d, -4*b], vvvvvvvvvvvvvvvvvvvv vvvvvvvvvvvv File "C: Item2Imp. py", line 42, in newton_raphson J = jacobian(x) vvvvvvvvvvv File "C: Item2Imp. py", line 51, in x_a, func_value_a, iterations_a = newton_raphson (func_a, jacobian_a, x0_a, le-14, 100) vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv ValueError: setting an array element with a sequence. The requested array has an inhomogeneous shape after 1 dimensions. The detected shape was (4, ) + inhomogeneous part.(a) Solution: [nan nan nan nan] (a) Function Value: nan (a) Relative Error: nan (a) Number of Iterations: 100 (b) Solution: [ 5.00000000e-01 -1. 30844385e-18 -5.23598776e-01] (b) Function Value: 1.7763568394002505e-15 (b) Relative Error: 1.0 (b) Number of Iterations: 8\f