Program (in a Python function; hint: extend the function from Assignment 1) a (mxm) LU factorization algorithm (using partial pivoting) for a matrix A and compute the PLU factors. The factorization is obtained by elimination steps k = 1,...,m-1 so that A-₁AV -k+1,...,m and where the multipliers mix are given by mi When km-1, A, is upper triangular, that is, U₁ = A). The lower triangular matrix is obtained using the multipliers m.,. that is Lij = mij Vi>j. Li,i = 1, and Lij = 0 Vi

(Numerical analysis) Here’s a challenging problem for those who know a little calculus. The Newton-Raphson method can be used to find the roots of any equation y(x)=0. In this method, the (i+1)stapproximation,xi+1,toarootofy(x)=0 is given in terms of the ith approximation, xi, by the following formula, where y’ denotes the derivative of y(x) with respect to x: xi+1=xiy(xi)/y(xi) For example, if y(x)=3x2+2x2,theny(x)=6x+2 , and the roots are found by making a reasonable guess for a first approximation x1 and iterating by using this equation: xi+1=xi(3xi2+2xi2)/(6xi+2) a. Using the Newton-Raphson method, find the two roots of the equation 3x2+2x2=0. (Hint: There’s one positive root and one negative root.) b. Extend the program written for Exercise 6a so that it finds the roots of any function y(x)=0, when the function for y(x) and the derivative of y(x) are placed in the code.

(Thermodynamics) The work, W, performed by a single piston in an engine can be determined by this formula: W=Fd F is the force provided by the piston in Newtons. d is the distance the piston moves in meters. a. Determine the units of W by calculating the units resulting from the right side of the formula. Check that your answer corresponds to the units for work listed in Table 1.1. b. Determine the work performed by a piston that provides a force of 1000 N over a distance of 15 centimeters.

(Statics) An annulus is a cylindrical rod with a hollow center, as shown in Figure 6.7. Its second moment of inertia is given by this formula: I4(r24r14) I is the second moment of inertia (m4). r2 is the outer radius (m). r1 is the inner radius (m). a. Using this formula, write a function called annulusMoment ( ) that accepts two double-precision numbers as parameters (one for the outer radius and one for the inner radius), calculates the corresponding second moment of inertia, and displays the result. b. Include the function written in Exercise 5a in a working program. Make sure your function is called from main(). Test the function by passing various data to it.