MAT331, Feb 4th, 2010Announcements: REU program - Wikis will be used for second project.Finish problem 4 Problem 53D plotsp:=2 *sin( x) - x^3 - 0.5:
plot(p,x=-5..5);
plot(p,x=-2..2);
fsolve(p,x);If a curve fits, wear it (Chapter 2 of the online book)EXERCISE: Suppose we measure the height of a baby, and we get the following results:When the baby is 1 day old, the height is 20.1 in.When the baby is 15 days old, the height is 20.6 in.When the baby is 20 days old, the height is 21.2 in.When the baby is 29 days old, the height is 21.9 in.JSFH1. Make a graph of the points determined by the above data2. What is the height of the baby at 25 days old?3. Find a polynomial passing through all the points.4. Find a line passing through all the points.5.EXERCISE: Given the points, [1,-1/2], [2,1], [2.4,1] and [3,1/3] find a piecewise linear function with domain the interval [1,3] passing through all of them.EXERCISE: Is it possible to find a linear function passing through [1,-1/2], [2,1] and [3,1/3]? Why?EXERCISE: Find the linear function passing closest to the points [1,-1/2], [2,1] and [3,1/3]EXERCISE: Find the quadratic function passing closest to the points [1,-1/2], [2,1] and [3,1/3]EXERCISE: Find the quadratic function passing closest to the points [1,2], [2.4,1.5] , [3,1] and [-0.4,3]EXERCISE: Is it possible to find a linear function passing through the following points? Justify your answer.EXERCISE: Find the linear function passing closest to the points [1,-1/2], [2,1], [2.4,1] and [3,1/3] . You can use the graph below to study the distance. pts:=[ [1,-1/2], [2,1], [2.4,1] , [3,1/3] ];Note: The graph below works only if all the vertical blue segments have length larger than 0. L:=x->x+3;with(plots):display(plot(pts, x = 0 .. 5, style = point, symbol = box, color = green), plot([L(x)], x = 0 .. 5, color = red), seq(plot([pts[i][1], t, t = pts[i][2] .. L(pts[i][1])], thickness = 3, color = blue), i = 1 .. 3));EXERCISE: Consider the pointsnewpts := [seq([i, 1.5*i+1], i = -1 .. 5), [8, 1.5*8+.5], seq([i, 1.5*i+1], i = 9 .. 11)];a. Find a polynomial passing through those points.b. Fit the points with a line using the least squares method. Use the graph of the previous exercise to visualize the line you found, and see how "close" the line and the points are.c. Fit the points with a quadratic function using the least squares method.d. Find a piecewise linear function passing through all the points.e. Plot the points and all the functions together.f. Compare the results. In particular, compute the value of the function you found in each case when x=12.3 Fitting points to a circleHere is a program that produces a list of points approximately in a circle. You do not have to worry about the details of this program (but you are wellcome to do it if you want), just run it to obtain the points. Seed:=readlib(randomize)():
with(stats):
myrand := (a,b) -> evalf(a+(b-a)*(rand()/10^12)):
points_wanted := 50:
circlePts := proc()
local fuzz, a, b, r,s,x,i;
a:= myrand(-3,3);
b:= myrand(-3,3);
r:= myrand(1,2);
fuzz := [ random[normald[0,0.1]](points_wanted)];
x := [ random[uniform[a-r,a+r]](points_wanted)];
for i from 1 to points_wanted do if (myrand(-10,10) > 0) then s[i]:=1 else s[i]:=-1 fi:od:
return [seq([x[i],b+s[i]*sqrt(r^2-(x[i]-a)^2)+fuzz[i]],i=1..points_wanted)];
end:myCirclePts:=circlePts():