MAT331,

Feb 16th, 2010

If a curve fits, wear it (Chapter 2 of the online book)

restart;

pts;

pts:=[[9.461625032, -15.00000412], [-4.026835626, 5.229076781], [-8.417983819, 17.28331107], [7.335537549, -10.75702150], [.667622007, 1.352242696], [-3.642662167, 8.570158849], [4.108941669, -6.649973580], [-3.592756472, 7.661525316], [0.121773384e-1, -.5948721361], [3.588065808, -7.386147022], [9.371019578, -17.50451275], [6.303804064, -11.16461096], [9.396335873, -15.62444142], [-.807790200, 2.555503241], [.782176902, -4.501842488], [3.640727999, -6.449402329], [1.668897236, -3.014773741], [3.704029489, -3.221704164], [-11.29937335, 17.53901720], [8.475940602, -17.91640911], [-3.666884970, 8.440428354], [-5.617368657, 12.00375535], [-5.626834452, 7.337908700], [-2.138751424, 5.625807468], [-2.268585475, 4.613703098]]:

plot(pts,style=point);

with(plots):

Note: The graph below works only if all the vertical blue segments have length larger than 0.

L:=x->(-x); display(plot(pts, x = 0 .. 5, style = point, symbol = box, color = green), plot([L(x)], x = -12 .. 12, color = red), seq(plot([pts[i][1], t, t = pts[i][2] .. L(pts[i][1])], thickness = 3, color = blue), i = 1 .. nops(pts)));

distancePointLine:=P->(m*P[1]+b-P[2])^2;

distancePointsLine:=pts->sum(distancePointLine(pts[i]),i=1..nops(pts));

DD:=distancePointsLine(pts);

poly:=simplify(DD);

plot3d(poly, m=-10..10,b=-10..10);

eq1:=diff(DD,m); eq2:=diff(DD,b);

sol:=solve({eq1,eq2},{m,b});

foundedLine:=x->subs(sol, m*x+b);

L:=foundedLine; display(plot(pts, x = 0 .. 5, style = point, symbol = box, color = green), plot([L(x)], x = -12 .. 12, color = red), seq(plot([pts[i][1], t, t = pts[i][2] .. L(pts[i][1])], thickness = 3, color = blue), i = 1 .. nops(pts)));

EXERCISE: Consider the points

newpts := [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

Problems 11 and 13

Fitting points to a circle

Here 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():

myCirclePts;

plot(myCirclePts,style=point);

d1:=P->((P[2]-b)^2+(P[1]-a)^2-r^2)^2;

dCircle:=simplify(sum(d1(myCirclePts[i]),i=1..nops(myCirclePts) ));

eq1:=diff(dCircle,a); eq2:=diff(dCircle,b); eq3:=diff(dCircle,r);

solCircle:=solve({eq1,eq2,eq3},{a,b,r});

Nik8JS9JImFHNiJeJCQhK2VIR28qKSEjNSQiKzwvJHlIIyEiKi9JImJHRiZeJCQhKyk+XD0wKEYqJCEra1A3XUNGKi9JInJHRiYkIiIhRjg8JS9GJV4kJCErYE13SjVGLSQiK15YIik+TEYqL0YvXiQkISs9SDkiNCJGLSQiKzdYLGlGRi1GNTwlL0YlJCErJz5lN0wpRiovRi8kIStMPmVtJSpGKkY1PCUvRiVeJEY8JCErXlgiKT5MRiovRi9eJEZCJCErN1gsaUZGLUY1PCUvRiVeJEYoJCErPC8keUgjRi0vRi9eJEYxJCIra1A3XUNGKkY1PCUvRiUkIStRPFZ5c0YqL0YvJCErJHBHX3gpRiovRjYkIiskbzBeIj1GLTwlRmpuRl1vL0Y2JCErJG8wXiI9Ri0=

solCircle[6];

PCUvSSJhRzYiJCErUTxWeXMhIzUvSSJiR0YlJCErJHBHX3gpRigvSSJyR0YlJCIrJG8wXiI9ISIq

theCircle:=subs(solCircle[6],(x-a)^2+(y-b)^2=r^2);

LywmKiQsJkkieEc2IiIiIiQiK1E8VnlzISM1RigiIiNGKCokLCZJInlHRidGKCQiKyRwR194KUYrRihGLEYoJCIrUyczWUgkISIq

graph1:=implicitplot(theCircle,x=-3..3,y=-3..3,color=blue,thickness=3);

LSUlUExPVEc2JS1JJ0NVUlZFU0c2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkc2IjYkWCxJKWFueXRoaW5nR0YoRio2Iy9JLnNvdXJjZV9ydGFibGVHRipYLEYsRipGKltnbCchJSIhISNleCJocSIjNDA1QTQwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDBDMDA0MjhGNUMyOEY1QzI5QkZGMjUxMERBNzJBNUFGNkMwMDNFM0MyMzU4QzIzNjhCRkY1MUVCODUxRUI4NTFGQzAwM0MwQUU4MjA0MDA3M0JGRjVFRjQ2RDMwMjNDOEJDMDAzM0U2QzA4QzI1QkY5QkZGOEY1QzI4RjVDMjhGNkMwMDI2NjAxRTYzODBGQzRCRkZDN0JBQTQ4MEFDMUMwQzAwMjUyODlGNEZGRDBCRkJGRkNDQ0NDQ0NDQ0NDQ0RDMDAyM0Q3MEEzRDcwQTNEQkZGRDEwN0M3N0Y4NjMzRkMwMDBGMDRCMzIyQjYzNkNDMDAwNTFFQjg1MUVCODUyQzAwMDUxRUI4NTFFQjg1MkMwMDEwMkM3MUFERjc5MDhCRkZERkVENjI4RjZFM0ZFQzAwMjNENzBBM0Q3MEEzREJGRkNDQ0NDQ0NDQ0NDQ0RDMDAyQjU0MjUyMkU5QTZFQkZGQTREMTU0MTc0MDdDOEMwMDM3RDRDNjk4MzZDQzBCRkY4RjVDMjhGNUMyOEY2QzAwM0U4QTdEMjZENUQxQ0JGRjdFMkQwOTc0OTFBRkZDMDA0MjhGNUMyOEY1QzI5QkZGNTFFQjg1MUVCODUxRkMwMDRCODA0QkEzOUU5MUZCRkYzOTdFODEyOTAxODgwQzAwNEVDNURFMjNEMTI3OEJGRjE0N0FFMTQ3QUUxNDhDMDA1M0JCMTkzRUUxRTBEQkZFRDlEQjYzNjE2QUI5MEMwMDU2NTVGM0Y0NzIxRTlCRkVBRTE0N0FFMTQ3QUUxQzAwNTdDNzAwNzIzQjBBNkJGRTU5OTYyQjc1Mjc2NzFDMDA1N0FFRjAwM0ZERDQ1QkZFMzMzMzMzMzMzMzMzM0MwMDU3QTQwMTNEQUEwRUJCRkREMjhCRENBRTY0QkI5QzAwNTUwQUFENUZGNUY3RUJGRDcwQTNENzBBM0Q3MEFDMDA1MzUyMUJBMTJFRUREQkZEMDc5RjU0M0Y0MzkxREMwMDRGQUZFQzgyNTRGRTZCRkJFQjg1MUVCODUxRUI4QzAwNEFEMTRGOUNDOUE3QUJGQjMwODFBMkYzMjM1MDRDMDA0ODY3NzgwNzFGMzc3M0ZBMjhCMjU2MEM0OTM2QkMwMDQyOEY1QzI4RjVDMjkzRkI3Qzg2NEQyMDcwRDVFQzAwM0YxNzY1OUMzNkI5RTNGQkVCODUxRUI4NTFFQjhDMDAzRDZFRTQ4QkQwRkU5M0ZDRjAxQUZGMzdGODVFQ0MwMDMzN0M5MTNCMkQ5QTYzRkQ3MEEzRDcwQTNENzBBQzAwMjlFNzg1ODAyQzk3NTNGRDhCQjE2MTE3OUFDOTlDMDAyNzM4QkI3RjFDNEVGM0ZEQURDQUY4MDQ4MjJGOEMwMDIzRDcwQTNENzBBM0QzRkUwNjZFNkRENDREMjgzQzAwMThBNUQ4RTVCNzIxMTNGRTMzMzMzMzMzMzMzMzNDMDAwRTFCNjUxNERDQzlDM0ZFNDQxNURDQTQ5RDY2OUMwMDA5NTc2MkFFNDYxMUYzRkU1MzBBQ0YxREVEMDJGQzAwMDUxRUI4NTFFQjg1MjNGRTc5MjZENUI0ODhDNzdCRkZFRkM2OUUwRDc3OTZGM0ZFQUI4NEI0MTU0Rjc5QkJGRkNDQ0NDQ0NDQ0NDQ0QzRkVBQzZFOEVBMDEyRkY2QkZGQ0JGOUQ2QUMzMjc1NzNGRUFFMTQ3QUUxNDdBRTFCRkZDQTdEMzdFMTQwMDRDM0ZFRDQyRkY3NkVCMEY1OEJGRkEyNjlFNzNDNzczMzEzRkVFNjRFMkFCQTc0RkFEQkZGOEY1QzI4RjVDMjhGNjNGRUY2Rjc0OEQ5QkJCNDJCRkY3NjVDRUMxQUYyNTRGM0ZGMDdBMEFGMzFCRDU2N0JGRjUxRUI4NTFFQjg1MUYzRkYwOUIwNzdBMzc3NTY2QkZGNDcyMTFCN0E4MTkzRDNGRjEzNURDNzI3MEE1NDNCRkYxNDdBRTE0N0FFMTQ4M0ZGMTM2QjBGN0Q2QjdFN0JGRjEzNkIwRjdENkI3RTczRkYxNDdBRTE0N0FFMTQ4QkZFRkI2MEUzMkE1OTY2RjNGRjE2MTBDRUUwOTNCOTBCRkVCMTQwNTYxMzEyRjcyM0ZGMTYyMjJDQjAwQUJGMkJGRUFFMTQ3QUUxNDdBRTEzRkYxNDdBRTE0N0FFMTQ4QkZFODU5QjgwNjJFQ0JBMzNGRjEwQTI3MTc0MDJCRTFCRkUzMzMzMzMzMzMzMzMzM0ZGMEY3MzI4Mzk1MTYzREJGRTI5MjNDMTE2NzlEMUMzRkYwMjJBMDNDQkFFMkEzQkZENzBBM0Q3MEEzRDcwQTNGRUVFMTdGMDk2MTlBNDBCRkNGNUQwNjYyRjcwQ0Q3M0ZFRDVFQTI4ODg0Nzc2M0JGQkVCODUxRUI4NTFFQjgzRkVBRTE0N0FFMTQ3QUUxM0Y5RTVCODM3QUU4QzVCNjNGRTkyMEFBMTU0MkQwQ0IzRkJFQjg1MUVCODUxRUI4M0ZFMzMzMzMzMzMzMzMzMzNGRDcwQTE0QjRCNTFGQUMzRkUzMzMyMDVCMUZFNzJGM0ZENzBBM0Q3MEEzRDcwQTNGRDcwQTNENzBBM0Q3MEEzRkUyNUE1QjBCRENGNzFEM0ZENEE5RUZDMjk4Q0M3QTNGRTMzMzMzMzMzMzMzMzMzRkJFQjg1MUVCODUxRUI4M0ZFNzIyMTAwMkZBQ0Y4NEJGQkIxM0NEODIwNENGRjAzRkVBNkNCNzIwRTQ3MTA4QkZCRUI4NTFFQjg1MUVCODNGRUFBMUQyQjQwQ0U3RDNCRkMyOUI1QzVGOEM2QTE4M0ZFQUUxNDdBRTE0N0FFMUJGRDcwQTNENzBBM0Q3MEEzRkVDQjQxNUMwMTlCNUU3QkZEQjQxNjgyNTlEMzMwMTNGRUNGQ0REMDg5MTI4RERCRkUzMzMzMzMzMzMzMzMzM0ZFREJENDQ4REM3NjU5RkJGRTYxMjUzNTc4RTgzQTAzRkVEQzA2N0QyNkZDQjRGQkZFQUUxNDdBRTE0N0FFMTNGRURDNUE4RTE0RDAzQkFCRkVENzI5QzZBMDA5NTkwM0ZFRDcyOUM2QTAwOTU5MEJGRjE0N0FFMTQ3QUUxNDgzRkVDQ0Q0MkJBQUE5MDM3QkZGMjBCNjIzQTA5MzA5MDNGRUM2OEFGRjkzMTE5NzJCRkY1MDUwNDJBQUJFRkNEM0ZFQUUxNDdBRTE0N0FFMUJGRjUxOERCMkZFNTFDREUzRkVBRDU4RDZBMDdBQTVGQkZGNTFFQjg1MUVCODUxRjNGRUFEMjE2MTNFNERENDJCRkY3QzE2NjJEODgwRDlGM0ZFODc4OEVFQTZDNDQzNEJGRjhGNUMyOEY1QzI4RjYzRkU3NjVBMTg5MjkyNzFGQkZGQTQ2NUI1N0MwRTFCNDNGRTVENDY0QzNGQ0E0QUZCRkZDN0IxOENCRkFDRDk5M0ZFMzMzMzMzMzMzMzMzM0JGRkNBQjFFNDE3RDE2NzEzRkUyRUZENjFDOTNDNjdCQkZGQ0NDQ0NDQ0NDQ0NDRDNGRTJDMDk2ODdCMTc1OEVCRkZFQjlBRkZCQTQxNzdCM0ZERUJEQ0EyQzAxMDFDMkMwMDA1MUVCODUxRUI4NTIzRkQ4MDQwMkE1MzRFQ0U4QzAwMDYxMUNFREU0OUYwMDNGRDc4M0M4QjZEMzBDN0NDMDAwNkY4M0E2NUJGMjg5M0ZENzBBM0Q3MEEzRDcwQUMwMDEzOUY5REFDQkEwNzUzRkNEREQwRTUwOTExMTkzQzAwMUVGMEQxMUQ1MEYxOTNGQkVCODUxRUI4NTFFQjhDMDAyMENENTY1NjRENUFDM0ZCOEE0RUExRDNFOENBNEMwMDIzRDcwQTNENzBBM0QzRkFENzREODZGQjMwODVDQzAwMkJBOTM0MkNEMkYwMUJGQUUyN0ZDMTk4MTBDNzhDMDAyRkMzRThDRTk4RkM5QkZCRUI4NTFFQjg1MUVCOEMwMDM0RjNGQTAxMDlCQjhCRkNDRjc4QjFEQUU5NjYzQzAwM0I3QTNFMkZCOUI0OUJGRDcwQTNENzBBM0Q3MEFDMDAzQ0Q2N0NFNEUyNDRCQkZEOUU2QUQxMkFEOTVGQUMwMDQyOEY1QzI4RjVDMjlCRkUyRjdEMDE5NUJCRkJGQzAwNDJCODNFQUY3NzdBNUJGRTMyOEZBOTE5MkM1NDNDMDA0MkMwQkU5OEIzQjBCQkZFMzMzMzMzMzMzMzMzM0MwMDQ1MTM4OTlDRDJBNjVCRkVBNDAzQzUxMUQ0MUYxQzAwNDU0ODlBQ0NGOTdEMkJGRUFFMTQ3QUUxNDdBRTFDMDA0NDA1ODgyN0EyOEQ1QkZGMTE4RTg5NEE1NDdFRkMwMDQzRjU2NDg3MThCODZCRkYxNDdBRTE0N0FFMTQ4QzAwNDI4RjVDMjhGNUMyOUJGRjI1MTBEQTcyQTVBRjYwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDBbZ2wnISUiISEjXXgiZHEiI0MwMDQyOEY1QzI4RjVDMjlCRkYyNTEwREE3MkE1QUY2QzAwM0UzQzIzNThDMjM2OEJGRjUxRUI4NTFFQjg1MUZDMDAzQzBBRTgyMDQwMDczQkZGNUVGNDZEMzAyM0M4QkMwMDMzRTZDMDhDMjVCRjlCRkY4RjVDMjhGNUMyOEY2QzAwMjY2MDFFNjM4MEZDNEJGRkM3QkFBNDgwQUMxQzBDMDAyNTI4OUY0RkZEMEJGQkZGQ0NDQ0NDQ0NDQ0NDREMwMDIzRDcwQTNENzBBM0RCRkZEMTA3Qzc3Rjg2MzNGQzAwMEYwNEIzMjJCNjM2Q0MwMDA1MUVCODUxRUI4NTJDMDAwNTFFQjg1MUVCODUyQzAwMTAyQzcxQURGNzkwOEJGRkRGRUQ2MjhGNkUzRkVDMDAyM0Q3MEEzRDcwQTNEQkZGQ0NDQ0NDQ0NDQ0NDREMwMDJCNTQyNTIyRTlBNkVCRkZBNEQxNTQxNzQwN0M4QzAwMzdENEM2OTgzNkNDMEJGRjhGNUMyOEY1QzI4RjZDMDAzRThBN0QyNkQ1RDFDQkZGN0UyRDA5NzQ5MUFGRkMwMDQyOEY1QzI4RjVDMjlCRkY1MUVCODUxRUI4NTFGQzAwNEI4MDRCQTM5RTkxRkJGRjM5N0U4MTI5MDE4ODBDMDA0RUM1REUyM0QxMjc4QkZGMTQ3QUUxNDdBRTE0OEMwMDUzQkIxOTNFRTFFMERCRkVEOURCNjM2MTZBQjkwQzAwNTY1NUYzRjQ3MjFFOUJGRUFFMTQ3QUUxNDdBRTFDMDA1N0M3MDA3MjNCMEE2QkZFNTk5NjJCNzUyNzY3MUMwMDU3QUVGMDAzRkRENDVCRkUzMzMzMzMzMzMzMzMzQzAwNTdBNDAxM0RBQTBFQkJGREQyOEJEQ0FFNjRCQjlDMDA1NTBBQUQ1RkY1RjdFQkZENzBBM0Q3MEEzRDcwQUMwMDUzNTIxQkExMkVFRERCRkQwNzlGNTQzRjQzOTFEQzAwNEZBRkVDODI1NEZFNkJGQkVCODUxRUI4NTFFQjhDMDA0QUQxNEY5Q0M5QTdBQkZCMzA4MUEyRjMyMzUwNEMwMDQ4Njc3ODA3MUYzNzczRkEyOEIyNTYwQzQ5MzZCQzAwNDI4RjVDMjhGNUMyOTNGQjdDODY0RDIwNzBENUVDMDAzRjE3NjU5QzM2QjlFM0ZCRUI4NTFFQjg1MUVCOEMwMDNENkVFNDhCRDBGRTkzRkNGMDFBRkYzN0Y4NUVDQzAwMzM3QzkxM0IyRDlBNjNGRDcwQTNENzBBM0Q3MEFDMDAyOUU3ODU4MDJDOTc1M0ZEOEJCMTYxMTc5QUM5OUMwMDI3MzhCQjdGMUM0RUYzRkRBRENBRjgwNDgyMkY4QzAwMjNENzBBM0Q3MEEzRDNGRTA2NkU2REQ0NEQyODNDMDAxOEE1RDhFNUI3MjExM0ZFMzMzMzMzMzMzMzMzM0MwMDBFMUI2NTE0RENDOUMzRkU0NDE1RENBNDlENjY5QzAwMDk1NzYyQUU0NjExRjNGRTUzMEFDRjFERUQwMkZDMDAwNTFFQjg1MUVCODUyM0ZFNzkyNkQ1QjQ4OEM3N0JGRkVGQzY5RTBENzc5NkYzRkVBQjg0QjQxNTRGNzlCQkZGQ0NDQ0NDQ0NDQ0NDRDNGRUFDNkU4RUEwMTJGRjZCRkZDQkY5RDZBQzMyNzU3M0ZFQUUxNDdBRTE0N0FFMUJGRkNBN0QzN0UxNDAwNEMzRkVENDJGRjc2RUIwRjU4QkZGQTI2OUU3M0M3NzMzMTNGRUU2NEUyQUJBNzRGQURCRkY4RjVDMjhGNUMyOEY2M0ZFRjZGNzQ4RDlCQkI0MkJGRjc2NUNFQzFBRjI1NEYzRkYwN0EwQUYzMUJENTY3QkZGNTFFQjg1MUVCODUxRjNGRjA5QjA3N0EzNzc1NjZCRkY0NzIxMUI3QTgxOTNEM0ZGMTM1REM3MjcwQTU0M0JGRjE0N0FFMTQ3QUUxNDgzRkYxMzZCMEY3RDZCN0U3QkZGMTM2QjBGN0Q2QjdFNzNGRjE0N0FFMTQ3QUUxNDhCRkVGQjYwRTMyQTU5NjZGM0ZGMTYxMENFRTA5M0I5MEJGRUIxNDA1NjEzMTJGNzIzRkYxNjIyMkNCMDBBQkYyQkZFQUUxNDdBRTE0N0FFMTNGRjE0N0FFMTQ3QUUxNDhCRkU4NTlCODA2MkVDQkEzM0ZGMTBBMjcxNzQwMkJFMUJGRTMzMzMzMzMzMzMzMzMzRkYwRjczMjgzOTUxNjNEQkZFMjkyM0MxMTY3OUQxQzNGRjAyMkEwM0NCQUUyQTNCRkQ3MEEzRDcwQTNENzBBM0ZFRUUxN0YwOTYxOUE0MEJGQ0Y1RDA2NjJGNzBDRDczRkVENUVBMjg4ODQ3NzYzQkZCRUI4NTFFQjg1MUVCODNGRUFFMTQ3QUUxNDdBRTEzRjlFNUI4MzdBRThDNUI2M0ZFOTIwQUExNTQyRDBDQjNGQkVCODUxRUI4NTFFQjgzRkUzMzMzMzMzMzMzMzMzM0ZENzBBMTRCNEI1MUZBQzNGRTMzMzIwNUIxRkU3MkYzRkQ3MEEzRDcwQTNENzBBM0ZENzBBM0Q3MEEzRDcwQTNGRTI1QTVCMEJEQ0Y3MUQzRkQ0QTlFRkMyOThDQzdBM0ZFMzMzMzMzMzMzMzMzMzNGQkVCODUxRUI4NTFFQjgzRkU3MjIxMDAyRkFDRjg0QkZCQjEzQ0Q4MjA0Q0ZGMDNGRUE2Q0I3MjBFNDcxMDhCRkJFQjg1MUVCODUxRUI4M0ZFQUExRDJCNDBDRTdEM0JGQzI5QjVDNUY4QzZBMTgzRkVBRTE0N0FFMTQ3QUUxQkZENzBBM0Q3MEEzRDcwQTNGRUNCNDE1QzAxOUI1RTdCRkRCNDE2ODI1OUQzMzAxM0ZFQ0ZDREQwODkxMjhEREJGRTMzMzMzMzMzMzMzMzMzRkVEQkQ0NDhEQzc2NTlGQkZFNjEyNTM1NzhFODNBMDNGRURDMDY3RDI2RkNCNEZCRkVBRTE0N0FFMTQ3QUUxM0ZFREM1QThFMTREMDNCQUJGRUQ3MjlDNkEwMDk1OTAzRkVENzI5QzZBMDA5NTkwQkZGMTQ3QUUxNDdBRTE0ODNGRUNDRDQyQkFBQTkwMzdCRkYyMEI2MjNBMDkzMDkwM0ZFQzY4QUZGOTMxMTk3MkJGRjUwNTA0MkFBQkVGQ0QzRkVBRTE0N0FFMTQ3QUUxQkZGNTE4REIyRkU1MUNERTNGRUFENThENkEwN0FBNUZCRkY1MUVCODUxRUI4NTFGM0ZFQUQyMTYxM0U0REQ0MkJGRjdDMTY2MkQ4ODBEOUYzRkU4Nzg4RUVBNkM0NDM0QkZGOEY1QzI4RjVDMjhGNjNGRTc2NUExODkyOTI3MUZCRkZBNDY1QjU3QzBFMUI0M0ZFNUQ0NjRDM0ZDQTRBRkJGRkM3QjE4Q0JGQUNEOTkzRkUzMzMzMzMzMzMzMzMzQkZGQ0FCMUU0MTdEMTY3MTNGRTJFRkQ2MUM5M0M2N0JCRkZDQ0NDQ0NDQ0NDQ0NEM0ZFMkMwOTY4N0IxNzU4RUJGRkVCOUFGRkJBNDE3N0IzRkRFQkRDQTJDMDEwMUMyQzAwMDUxRUI4NTFFQjg1MjNGRDgwNDAyQTUzNEVDRThDMDAwNjExQ0VERTQ5RjAwM0ZENzgzQzhCNkQzMEM3Q0MwMDA2RjgzQTY1QkYyODkzRkQ3MEEzRDcwQTNENzBBQzAwMTM5RjlEQUNCQTA3NTNGQ0RERDBFNTA5MTExOTNDMDAxRUYwRDExRDUwRjE5M0ZCRUI4NTFFQjg1MUVCOEMwMDIwQ0Q1NjU2NEQ1QUMzRkI4QTRFQTFEM0U4Q0E0QzAwMjNENzBBM0Q3MEEzRDNGQUQ3NEQ4NkZCMzA4NUNDMDAyQkE5MzQyQ0QyRjAxQkZBRTI3RkMxOTgxMEM3OEMwMDJGQzNFOENFOThGQzlCRkJFQjg1MUVCODUxRUI4QzAwMzRGM0ZBMDEwOUJCOEJGQ0NGNzhCMURBRTk2NjNDMDAzQjdBM0UyRkI5QjQ5QkZENzBBM0Q3MEEzRDcwQUMwMDNDRDY3Q0U0RTI0NEJCRkQ5RTZBRDEyQUQ5NUZBQzAwNDI4RjVDMjhGNUMyOUJGRTJGN0QwMTk1QkJGQkZDMDA0MkI4M0VBRjc3N0E1QkZFMzI4RkE5MTkyQzU0M0MwMDQyQzBCRTk4QjNCMEJCRkUzMzMzMzMzMzMzMzMzQzAwNDUxMzg5OUNEMkE2NUJGRUE0MDNDNTExRDQxRjFDMDA0NTQ4OUFDQ0Y5N0QyQkZFQUUxNDdBRTE0N0FFMUMwMDQ0MDU4ODI3QTI4RDVCRkYxMThFODk0QTU0N0VGQzAwNDNGNTY0ODcxOEI4NkJGRjE0N0FFMTQ3QUUxNDhDMDA0MjhGNUMyOEY1QzI5QkZGMjUxMERBNzJBNUFGNi1JJ0NPTE9VUkdGJzYmSSRSR0JHRickIiIhRjVGNCQiKisrKysiISIpLUkqVEhJQ0tORVNTR0YnNiMiIiQtSStBWEVTTEFCRUxTR0YnNiRJInhHRipJInlHRio=

display(graph1,plot(myCirclePts,style=point));

LSUlUExPVEc2KC0lJ0NVUlZFU0c2JTdlcTckJCEkXyMhIiMkITFXNDg2JSl5VzYhIzo3JCQhMikpNFckby1AJ1sjISM7JCEkSyIhIiM3JCQhMk0wYXQlRzNwQyEjOyQhMm4lZmtfciI0UCIhIzs3JCQhMnhlOGtcemFTIyEjOyQhJGMiISIjNyQkITF4UDY1JDMpKkgjISM6JCEySkEnKSkqbyI+IXkiISM7NyQkITIlcCZcWT4tLkgjISM7JCEjPSEiIjckJCEkRyMhIiMkITJrSFJEJlxfOz0hIzs3JCQhMiZlI1JcIzRMPEAhIzskISQvIyEiIzckJCEkLyMhIiMkITFScDREak5FQCEjOjckJCEyUjBldyZmcnU9ISM7JCEkRyMhIiM3JCQhIz0hIiIkITJ1cWpOUzAmUUIhIzs3JCQhMiV5ejYqND5RayIhIzskITI8LSMpMyE0PU9DISM7NyQkISRjIiEiIyQhMTNsLVg4ZylbIyEjOjckJCExcy80NFooR1wiISM6JCEkXyMhIiM3JCQhJEsiISIjJCEyTSV6KW93XykqZSMhIzs3JCQhLyZvZmInZUM3ISM4JCEvOi5XTVQ6RSEjODckJCEkMyIhIiMkITJPW0JtSFpUbCMhIzs3JCQhMSN5STxDPl1EKiEjOyQhMUFwI2UyKVx1RSEjOjckJCEjJSkhIiMkITFjJ1I3YmdkbyMhIzo3JCQhMGghPSRIUShcbiEjOiQhMVI+b3FoLSZvIyEjOjckJCEiJyEiIiQhMjFZTzJgI3AlbyMhIzs3JCQhMTpVRWpsNmNYISM7JCEyI3lOblYkKVFrRSEjOzckJCEjTyEiIyQhMiUpKlE2TkslNGwjISM7NyQkITIlM0FYKEhQV2QjISM8JCExelpEcWliQUUhIzo3JCQhIzchIiMkITIjbz5Qa0VeJWUjISM7NyQkITB2Jz4tIlFVVighIzskITJGLnkqPXdsbEQhIzs3JCQiMU1ONXpceUBPISM8JCEkXyMhIiM3JCQiMW1FeTw8OiFIKiEjPCQhMm0jeTw8OiFIXCMhIzs3JCQiIzchIiMkITFYKT0leW0lKnpDISM6NyQkIjJWXzlmIyoqUUFDISM8JCEyQVgiZiMqKlFBUyMhIzs3JCQiI08hIiMkITIyOW9dKHpQRkIhIzs3JCQiMUFpLWQkKT1rUSEjOyQhMUFFcU4pPWtJIyEjOjckJCIydGohKTR4Jj4oPiUhIzwkISRHIyEiIzckJCIxSHJtbERoRF4hIzskITJGcm1sRGhEPiMhIzs3JCQiIichIiIkITJtWGlXLTYtNiMhIzs3JCQiMT14ZkhIekhqISM7JCEyPHhmSEh6SDIjISM7NyQkIjE5KCk0OiU9PmknISM7JCEkLyMhIiM3JCQiMVIiUSE0VENtdCEjOyQhMTlRITRUQ20kPiEjOjckJCIxeCx4PiNvKlwkKSEjOyQhIz0hIiI3JCQiMXJJalIpNHlPKSEjOyQhMTJMJ1IpNHkneiIhIzo3JCQiIyUpISIjJCEyWFJMJHpKKDR6IiEjOzckJCIxYnYtY1hHVyIqISM7JCEyYnYtY1hHV2oiISM7NyQkIjFSVXhoNjopXCohIzskISRjIiEiIzckJCIxSyJSYyJRYkIpKiEjOyQhMksiUmMiUWJCWSIhIzs3JCQiMilRUGspZiZ6SDUhIzskISRLIiEiIzckJCIwbSgpeSMqW3kuIiEjOSQhMG0oKXkjKlt5RiIhIzk3JCQiMTdweihvXGMyIiEjOiQhJDMiISIjNyQkIjJaJjRxakImZTIiISM7JCEyWiY0cWpCJmUyIiEjOzckJCIkMyIhIiMkITE0R0twY3Q0KiohIzs3JCQiMjRDcTkoUj4nMyIhIzskITE0Q3E5KFI+WSkhIzs3JCQiMkNZP0MnKmVrMyIhIzskISMlKSEiIzckJCIkMyIhIiMkITBpZk8yPyY0dyEjOjckJCIyTiRROydveVwxIiEjOyQhIichIiI3JCQiMipbTEYzNE5nNSEjOyQhMShbTEYzNE4hZSEjOzckJCIxI2V4YGZgJTM1ISM6JCEjTyEiIzckJCIxLGIwaFJFXScqISM7JCEyM11iNSdSRV1DISM8NyQkIjEpKT15NzYteSIqISM7JCEjNyEiIzckJCIjJSkhIiMkIjInPlx2bHJma0ghIz03JCQiMSxmPT5OUF95ISM7JCIjNyEiIzckJCIiJyEiIiQiMShwNUIpRyEqKmYkISM7NyQkIjFrUmZXLCIqKipmISM7JCIjTyEiIzckJCIjTyEiIyQiMXNBLkN1SE5kISM7NyQkIjFMNnpZNXNHSyEjOyQiIichIiI3JCQiIzchIiMkIjEqMyhwYi0zSHMhIzs3JCQhMiZvcztpIjR4MCIhIzwkIjFvcztpIjR4RCkhIzs3JCQhIzchIiMkIjEpXGRbInpgQSQpISM7NyQkITJLW3NOT2lPWCIhIzwkIiMlKSEiIzckJCEjTyEiIyQiMSlmX2tGSSlwKikhIzs3JCQhMSQqSD5kNm5lVSEjOyQiMSQqSD5kNm5lISohIzs3JCQhIichIiIkIjF4LXQheVJOSCohIzs3JCQhMWU9Jz4xcXQqbyEjOyQiMWY9Jz4xcXRIKiEjOzckJCEjJSkhIiMkIjE5ZmcpeiR5LiQqISM7NyQkITFqSm0nKmZTLSMqISM7JCIxakptJypmUy0jKiEjOzckJCEkMyIhIiMkIjEzJnorTGkwKyohIzs3JCQhMkdkZUpBeng3IiEjOyQiMUdkZUpBengpKSEjOzckJCEyb1kneXdZczg4ISM7JCIjJSkhIiM3JCQhMmwnejQzJW8mPTghIzskIjFrJ3o0MyVvJlEpISM7NyQkISRLIiEiIyQiMV8lejQkSFgiUSkhIzs3JCQhMjs8Ql1lO1pbIiEjOyQiMTs8Ql1lO1p3ISM7NyQkISRjIiEiIyQiMUJ5RXk3YzZ0ISM7NyQkITJrJSk+XSZwPFU7ISM7JCIxaiUpPl0mcDwjbyEjOzckJCEwPlVcJUgwIXkiISM5JCIiJyEiIjckJCEybWlvNSFweCJ6IiEjOyQiMW1pbzUhcHgiZiEjOzckJCEjPSEiIiQiMTcsQSR5IzRnZSEjOzckJCEyWzNtZyVRTD8+ISM7JCIydyUzbWclUUwhWyEjPDckJCEkLyMhIiMkIjJiQ0hLclpDdiQhIzw3JCQhMilwSVhmJT11LyMhIzskIjIkKXBJWGYlPXVPISM8NyQkITEqSCQqM0tdVzAjISM6JCIjTyEiIzckJCEwWTc3ZjNMOiMhIzkkIjIuZ0NAImYzTEIhIzw3JCQhMUhIOVpUc1RBISM6JCIjNyEiIzckJCExL0BxJlJtaUQjISM6JCIxVi9AcSZSbWkqISM8NyQkISRHIyEiIyQiMTR3MjImM0t2JiEjPDckJCExcyplLz4sNk0jISM6JCExMEc1YTQpKSopZSEjPDckJCExZFNcIT1tSlAjISM6JCEjNyEiIzckJCExRixIdWJwOEMhIzokITFFKCk0ZFUvakEhIzs3JCQhMih5NXRkIW9ZWSMhIzskISNPISIjNyQkITJ2K2AtbCZIdkMhIzskITJhIypwdVxWcS8lISM8NyQkISRfIyEiIyQhMXA0XWxmXUZmISM7NyQkITJpK2ZxcVo3XyMhIzskITFQKjQlSEhfKClmISM7NyQkITEmUXNlNDI6XyMhIzokISInISIiNyQkITJRUCRvUShlJ1JEISM7JCExaWk7OEVULiMpISM7NyQkITE2Vkk+I3k3YSMhIzokISMlKSEiIzckJCEyWFo1KT0pPTlgIyEjOyQhMmFfKj0iPSJlbzUhIzs3JCQhMig+KXkvSEU0YCMhIzskISQzIiEiIzckJCEkXyMhIiMkITFXNDg2JSl5VzYhIzotJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIjNSEiIi0lKlRISUNLTkVTU0c2IyIiJC0lJ0NVUlZFU0c2JTdUNyQkIStCeS5QNSEiKiQhKzMhKSoqZUQhIio3JCQhKyh6I3pBOyEiKiQiKyUqcGQicCghIzU3JCQhK0pueHk1ISIqJCErVzI0QEYhIio3JCQhKi8yRUsqISIqJCIqS3QrLSIhIik3JCQhK00vKnAvIiEiKiQhK0UhKUhdRSEiKjckJCEqTSZScHghIiokIStvO2ZxRSEiKjckJCEoZnVfJiEiKiQhKyM+a0hYIyEiKjckJCEpRWEjPikhIikkIitPMVMicCohIzU3JCQhKldeLkwqISIqJCErLTw4cUUhIio3JCQhKzwlPisoPiEiKiQhKyRbVjI/IyEiKjckJCErOTsmSE8jISIqJCEreS1DZDwhIzU3JCQhKkdHKFJMISIqJCIrWGlPOCkpISM1NyQkIitZVVRMNSEiKiQhK14oKltINyEiKjckJCErWXQtJFsjISIqJCErLWZQKSp6ISM1NyQkIit1cXkmNCIhIiokISsmKUdhbnchIzU3JCQiKiZHSCM9KSEiKiQhKlZxcUwpISM1NyQkISspWycqKnk1ISIqJCEqYUk1biMhIik3JCQhK09eKnAtIyEiKiQhKyxyZ2g/ISIqNyQkISg5K2UiISInJCErVWs9VUMhIio3JCQhK0JqZSNHIyEiKiQiKyJHaVtNIyEjNjckJCEqTlNWYiYhIiokISooUi8hbyMhIik3JCQiKjRPXWckISIqJCIrWE53O1shIzU3JCQhKzl4eF8/ISIqJCErUDJQN0AhIio3JCQhKnAiUkAnKiEiKiQiK1tMXTQ1ISIqNyQkIikiKlxycSEiKSQiKy5cdy43ISM1NyQkISpwaXdLIiEiKSQiK25EXjwhKSEjNTckJCErXydSJXk+ISIqJCErSERbR0IhIio3JCQhKyp6JCpwayIhIiokIStvQlhJQiEiKjckJCErKSllU21DISIqJCErT1hmM1MhIzU3JCQhKiVbNCYzJyEiKiQiK0Y3LkYqKSEjNTckJCErVVs6WT8hIiokIStLMkJ5QCEiKjckJCIqKVtqKCp6ISIqJCEqUnAmUTshIzU3JCQhKyIqNDBWQiEiKiQhK2haOTY6ISIqNyQkISpmYGtTJSEiKiQhKypcRF95IyEiKjckJCIrJVttJTM1ISIqJCErdHg7U1chIzU3JCQhKk1xLypSISIqJCIrJ0c3WWEqISM1NyQkISouIz0rYiEiKiQiKzdhbnA1ISIqNyQkIipWKnA9ViEiKiQhK3VlRWRCISIqNyQkIilQWlUnKSEiKiQhKnBFW2YjISIpNyQkISptXjs8IiEiKiQhK2Z2UDhEISIqNyQkIStIbyhRLiMhIiokIitMNjlKSCEjNTckJCEqYjYyTyohIiokIStGKlsnW0YhIio3JCQhK0xjLEY7ISIqJCIrZUkjUlEnISM1NyQkISpgLHInZSEiKiQhKyhRXHZzIyEiKjckJCEqaTNBdiMhIiokIitFMCdcLCkhIzU3JCQhKndUPylcISIqJCErXiRRKVJGISIqNyQkISosZG9AKCEiKiQiK0d5YCZSKiEjNTckJCErN3MjNFEiISIqJCIrdmwkPVsqISM1NyQkISt5STIsOiEiKiQhK09sPyNlIyEiKjckJCErImYtNzkjISIqJCEqRUxAJD4hIiktJSZDT0xPUkc2JiUkUkdCRyQiIzUhIiIkIiIhISIiJCIiISEiIi0lJlNUWUxFRzYjJSZQT0lOVEctJiUmX0FYSVNHNiMiIiI2Ji0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiIhISIiLSUqTElORVNUWUxFRzYjIiIhLSUqVEhJQ0tORVNTRzYjIiIiLSUtVFJBTlNQQVJFTkNZRzYjJCIiISEiIi0mJSZfQVhJU0c2IyIiIzYmLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiEtJSpUSElDS05FU1NHNiMiIiItJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSUrQVhFU0xBQkVMU0c2JS1JI21pRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjVRIng2Ii8lJ2ZhbWlseUdRKERFRkFVTFQ2Ii8lJXNpemVHUSMxMDYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSV0cnVlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lLG1hdGh2YXJpYW50R1EnaXRhbGljNiItSSNtaUc2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY1USJ5NiIvJSdmYW1pbHlHUShERUZBVUxUNiIvJSVzaXplR1EjMTA2Ii8lJWJvbGRHUSZmYWxzZTYiLyUnaXRhbGljR1EldHJ1ZTYiLyUqdW5kZXJsaW5lR1EmZmFsc2U2Ii8lKnN1YnNjcmlwdEdRJmZhbHNlNiIvJSxzdXBlcnNjcmlwdEdRJmZhbHNlNiIvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXTYiLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV02Ii8lJ29wYXF1ZUdRJmZhbHNlNiIvJStleGVjdXRhYmxlR1EmZmFsc2U2Ii8lKXJlYWRvbmx5R1EmZmFsc2U2Ii8lKWNvbXBvc2VkR1EmZmFsc2U2Ii8lKmNvbnZlcnRlZEdRJmZhbHNlNiIvJStpbXNlbGVjdGVkR1EmZmFsc2U2Ii8lLHBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lNnNlbGVjdGlvbi1wbGFjZWhvbGRlckdRJmZhbHNlNiIvJSxtYXRodmFyaWFudEdRJ2l0YWxpYzYiLSUlRk9OVEc2JSUoREVGQVVMVEclKERFRkFVTFRHIiM1LSUlUk9PVEc2Jy0lKUJPVU5EU19YRzYjJCIjISohIiItJSlCT1VORFNfWUc2IyQiJCsiISIiLSUtQk9VTkRTX1dJRFRIRzYjJCIlNVEhIiItJS5CT1VORFNfSEVJR0hURzYjJCIlK1EhIiItJSlDSElMRFJFTkc2Ig==

JSFH