CPSC120A
Fundamentals of Computer Science I

Scientific Assignment 3

Use Python to estimate integrals.

Given an arbitrary mathematical function, being able to estimate the area under the curve generated is a really important skill to have. In Physics, this value can tell you a lot about the motion of some projectile. In Chemistry, they can be used to define kinetics for enzymes. They even have use in Accounting, where they can be used to calculate compounded interest on loans. In this assignment, you will write a program that uses the Riemann Sum approximation of the area under a curve.


Details

Using the Turtle module, you are going to write a program that can plot an arbitrary polynomial. To accomplish this, you should define a function called polynomial(x), which takes a floating point value x, and returns the result of calculating the value of the polynomial with the value of x. You should define some polynomial inside of this function. As a starter, use the following polynomial:

$$f(x) = x ^ 3 - 2 \cdot x ^ 2 + 6 \cdot x + 5$$

Using the defined function above, plot the polynomial on the turtle window. In order to make the plotted function easier to view, you should scale the output of your function up, similar to what we did for the sine function activity.

Once you have plotted the function, you are going to then compute and plot a series of rectangles to estimate the area under the curve for some range of values from the domain of the previously defined polynomial. Define a function called estimate_area(start, end, size), which takes integers for the starting x value of the range, the ending x value of the range, and the width of the rectangle used to estimate the area under the curve.

This function should compute a series of rectangles of the size specified, with a height determined by the value of the polynomial at one of the points specified in that range. Your estimate of area under the curve is just the sum of the area of all of the rectangles combined. Print the estimated value of the area, and draw the rectangles you computed on the turtle window as well.


"Hacker" Prompt

Each week, additional exercises related to the assignment will be provided at the end. These exercises are typically more challenging than the regular assignment. Bonus points will be provided to students who complete any of the "Hacker" level assignments.

  1. Even Better Estimates: You likely assumed that the value you computed from the polynomial was one of the corners of the rectangle. This provides a decent estimate, but ultimately will always be an over (or under, depending on which corner, etc) estimate of what the value really is. Alter your program so that shifts where you compute the height of the rectangle so that you use the center of the rectangle.

  2. Convergence: The area under the curve can be computed by taking the integral of the polynomial. This can be best approximated by computing the area where the size of the rectangles approaches 0. While this is "impossible" for you to currently compute in python, we can compute how much the integral changes as we decrease the width of the rectangles. Use a for loop to compute and print a table of ever decreasing widths of rectangles, and their associated estimate of the area.

  3. Square Roots: There are some functions that we can plot, but have some issues that cause Python to not be able to compute values for an area under a curve. For example, \(\sqrt{x}\) cannot be computed, because you cannot compute the square root of negative values, and the fact that it is not, technically, a function (\(\pm y = \sqrt{x}\)). Create a second program that allows you to plot the square root function.


Grading

The assignment will be graded on the following requirements according to the course’s programming assignment rubric.

Functionality (75%): A functional program will:

Style (25%): A program with good style will:

Submission

Submit your program as a .py file on the course Inquire page before class on Wednesday October 7th.