User built-in commands or intrinsic functions

Question

For this assignment, you will need to write some short codes to complete the following problems.  You are free to write in any language you prefer, including MATLAB, provided you do not use any built-in commands or intrinsic functions to shortcut the methods you are asked to investigate (though you may use them to verify your answers).  For instance, you may not use MATLAB's fzero command to find the roots of an equation.  The only restrictions on your codesare that they must write your output to a file, and numbers in your results must be properly formatted with an appropriate number of significant figures.  Use exponential notation where appropriate.Please turn in both hard and electronic copies of both your codes and your results.  Your codes and results may be combined into a single file for electronic submission, provided that it is clear which results correspond to which code.

1.  Root finding has its place in several applications in the areas of vibration (modal analysis of wing flutter), acoustics (microphone and speaker placement for acoustic signature measurements), and computational fluid dynamics (matching governing equations to proper numerical schemes for aerodynamic calculations).

Determine the highest real root, out to exactly three decimal points, of the polynomial

                                                                          f(x) = 1.5x³ – 8.1x² + 18.6x - 10

                a) graphically (no code required; simply plot; you may use any function of any software

                                                to do this)

                b) using fixed-point iteration (five iterations, x_o = 1)

                c) using Newton-Raphson (five iterations, x_o = 1)

                d) using Secant method (five iterations, x_-1 = 0, x­_o = 1)

 

Compute the approximate percent relative errors out to exactly four decimal points for each iteration of your solutions as described in the notes.

2.  There are many applications where numerical integration is required in aerospace engineering.  Examples include integrating a pressure distribution about an airfoil to find lift and drag, integrating a time-dependent drag profile on a rocket to determine time of flight and burnout altitude, and integrating a probability distribution function to find the average fuel drop diameter in an aircraft engine.  In many instances, the function to be integrated may be a set of discrete data points rather than being analytical, but the technique is similar.

Integrate the following function both analytically and numerically:

                                                                                       I(f) =  ₀∫³x² e^-xdx

 

Use both the trapezoidal and Simpson's 1/3 rules to numerically integrate the function out to exactly seven decimal points.  For both cases, start with n = 1, then double the number of segments until your last computation includes 512 segments.  Compute the percent relative errors, out to exactly four decimal points, for the numerical results (as compared to the analytic answer).

%Errot = (Numerical Answer - Analytic Answer)/Analytic Answer