(*^ ::[ Information = "This is a Mathematica Notebook file. It contains ASCII text, and can be transferred by email, ftp, or other text-file transfer utility. It should be read or edited using a copy of Mathematica or MathReader. If you received this as email, use your mail application or copy/paste to save everything from the line containing (*^ down to the line containing ^*) into a plain text file. On some systems you may have to give the file a name ending with ".ma" to allow Mathematica to recognize it as a Notebook. The line below identifies what version of Mathematica created this file, but it can be opened using any other version as well."; FrontEndVersion = "NeXT Mathematica Notebook Front End Version 2.2"; NeXTStandardFontEncoding; fontset = title, inactive, noPageBreakBelow, noPageBreakInGroup, nohscroll, preserveAspect, groupLikeTitle, center, M7, bold, L1, e8, 24, "Times"; ; fontset = subtitle, inactive, noPageBreakBelow, noPageBreakInGroup, nohscroll, preserveAspect, groupLikeTitle, center, M7, bold, L1, e6, 18, "Times"; ; fontset = subsubtitle, inactive, noPageBreakBelow, noPageBreakInGroup, nohscroll, preserveAspect, groupLikeTitle, center, M7, italic, L1, e6, 14, "Times"; ; fontset = section, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeSection, grayBox, M22, bold, L1, a20, 18, "Times"; ; fontset = subsection, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeSection, blackBox, M19, bold, L1, a15, 14, "Times"; ; fontset = subsubsection, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeSection, whiteBox, M18, bold, L1, a12, 12, "Times"; ; fontset = text, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = smalltext, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 10, "Times"; ; fontset = input, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeInput, M42, N23, bold, L1, 12, "Courier"; ; fontset = output, output, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeOutput, M42, N23, L-5, 12, "Courier"; ; fontset = message, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeOutput, M42, N23, L1, 12, "Courier"; ; fontset = print, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeOutput, M42, N23, L1, 12, "Courier"; ; fontset = info, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeOutput, M42, N23, L1, 12, "Courier"; ; fontset = postscript, PostScript, formatAsPostScript, output, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeGraphics, M7, l34, w282, h287, L1, 12, "Courier"; ; fontset = name, inactive, noPageBreakInGroup, nohscroll, preserveAspect, M7, italic, B65535, L1, 10, "Times"; ; fontset = header, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, italic, L1, 12, "Times"; ; fontset = leftheader, L0, 12; fontset = footer, inactive, nohscroll, noKeepOnOnePage, preserveAspect, center, M7, italic, L1, 12, "Times"; ; fontset = leftfooter, L0, 12; fontset = help, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = clipboard, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = completions, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12, "Courier"; ; fontset = special1, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = special2, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = special3, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = special4, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = special5, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; currentKernel; ] :[font = title; inactive; preserveAspect; startGroup] GEOMASS :[font = section; inactive; preserveAspect; startGroup] BRIEF ABSTRACT :[font = subsection; inactive; preserveAspect; endGroup] Determine the gravitational force at a point from a nearby mass and determine the underlying mass distribution given force information. :[font = section; inactive; Cclosed; preserveAspect; startGroup] GENERAL INFORMATION :[font = subsection; inactive; preserveAspect; endGroup] FileName: GEOMASS Full title: Determining the mass distribution of a bar from observed forces using physical principles. Last Revision Date: 29 May 1996. Developer: Brian J. Winkel, Department of Mathematical Sciences, United States Military Academy, West Point NY 10996 USA. Phone: 914-938-3200. Email: ab3646@usma2.usma.edu. FAX: 914-938-2409. Contact: Brian J. Winkel, Department of Mathematical Sciences, United States Military Academy, West Point NY 10996 USA. Phone: 914-938-3200. Email: ab3646@usma2.usma.edu. FAX: 914-938-2409. Aaron D. Klebanoff, Department of Mathematics, Rose-Hulman Institute of Technology, Terre Haute IN 47803 USA. Phone: 812-877-8151. Email: Aaron.Klebanoff@Rose-Hulman.Edu. FAX: 812-877-3198. Support: The production of this material is supported by the National Science Foundation under Division of Undergraduate Education grant DUE-9352849: Development Site for Complex, Technology-Based Problems in Calculus with Applications in Science and Engineering and the Arvin Foundation of Columbus IN. :[font = section; inactive; Cclosed; preserveAspect; startGroup] STATEMENT OF PROBLEM :[font = subsection; inactive; Cclosed; preserveAspect; startGroup] Introduction :[font = subsubsection; inactive; preserveAspect] "The general problem of geological prospecting is to determine the location, shape, and constitution of subterranean bodies from measurements at the earth's surface. We will treat a much simplified one-dimensional model. Namely, we suppose that on a horizontal line (of length 1 m) measurements are made of the vertical component of gravitational force due to a variable distribution of mass along a parallel line one unit (m) below the surface." [1, p. 9] :[font = subsubsection; inactive; preserveAspect] From Newton's law of gravitation it is known that there is a force of attraction between two masses, say m1 and m2 (in kg), at a distance of r m apart equal to F = G (m1 m2)/r^2 and the direction of the force is along the line connecting the masses. :[font = subsubsection; inactive; preserveAspect; endGroup] G is the gravitational constant in N m^2/kg^2 or m^3/(kg s^2). Here G = 6.67 x 10^(-11) m^3/(kg s^2). :[font = subsection; inactive; Cclosed; preserveAspect; startGroup] Forward Problem - Problem 1 :[font = subsubsection; inactive; preserveAspect; startGroup] We have the situation (see sketch below) in which there is a bar (1 m long) of mass density x(t) = 100 + 90 t^3 kg/m running along the x-axis, beginning at the origin (0, 0) and going to (1, 0), i.e., t goes from 0 to 1. We consider another unit point mass (1 kg) at a point along the line y = 1, say (s, 1). Distances in meters. Determine the vertical component of the gravitational force between the bar and the unit point mass (1 kg) particle at each point (s, 1) (0 <= s <= 1). :[font = postscript; PostScript; formatAsPostScript; output; inactive; preserveAspect; pictureLeft = 59; pictureTop = 2; pictureWidth = 290; pictureHeight = 218; endGroup; endGroup] %! %%Creator: Mathematica %%AspectRatio: 0.752212 MathPictureStart /Courier findfont 10 scalefont setfont % Scaling calculations 0 1 0 1 [ [ 74 362 0 0 ] [ 526 702 0 0 ] ] MathScale % Start of Graphics %!PS-Adobe-2.0 EPSF-2.0 %%Title: %%Creator: Draw %%CreationDate: Thu Jan 11 18:21:20 1996 %%For: me %%DocumentFonts: (atend) %%Pages: 0 0 %%BoundingBox: 74 362 526 702 %%NXNextStepVersion: 3.0 %%EndComments %%BeginProcSet: /usr/lib/NextStep/printPackage.ps 3.0 %! % NeXT Printing Package % Version: 3.1 % Copyright: 1988, NeXT, Inc. /__NXdef{1 index where{pop pop pop}{def}ifelse}bind def /__NXbdef{1 index where{pop pop pop}{bind def}ifelse}bind def /UserObjects 10 array __NXdef /defineuserobject{ exch dup 1 add dup UserObjects length gt{ array dup 0 UserObjects putinterval /UserObjects exch def }{pop}ifelse UserObjects exch 3 -1 roll put }__NXbdef /undefineuserobject{UserObjects exch null put}__NXbdef /execuserobject{UserObjects exch get exec}__NXbdef /__NXRectPath{4 2 roll moveto 1 index 0 rlineto 0 exch rlineto neg 0 rlineto closepath}__NXbdef /__NXProcessRectArgs{ 1 index type /arraytype eq{ exch 0 4 2 index length 1 sub{ dup 3 add 1 exch{1 index exch get exch}for 5 1 roll 5 index exec }for pop pop }{exec}ifelse }__NXbdef /rectfill{gsave newpath {__NXRectPath fill} __NXProcessRectArgs grestore}__NXbdef /rectclip{newpath {__NXRectPath} __NXProcessRectArgs clip newpath}__NXbdef /rectstroke{ gsave newpath dup type /arraytype eq{dup length 6 eq}{false}ifelse{ {gsave __NXRectPath null concat stroke grestore} dup length array cvx copy dup 2 4 -1 roll put __NXProcessRectArgs }{{__NXRectPath stroke} __NXProcessRectArgs}ifelse grestore }__NXbdef /_NXLevel2 systemdict /languagelevel known {languagelevel 2 ge}{false}ifelse __NXdef /xyshow{ 0 1 3 index length 1 sub{ currentpoint 4 index 3 index 1 getinterval show 3 index 3 index 2 mul 1 add get add exch 3 index 3 index 2 mul get add exch moveto pop }for pop pop }__NXbdef /xshow{ 0 1 3 index length 1 sub{ currentpoint 4 index 3 index 1 getinterval show exch 3 index 3 index get add exch moveto pop }for pop pop }__NXbdef /yshow{ 0 1 3 index length 1 sub{ currentpoint 4 index 3 index 1 getinterval show 3 index 3 index get add moveto pop }for pop pop }__NXbdef /arct{arcto pop pop pop pop}__NXbdef /setbbox{pop pop pop pop}__NXbdef /ucache{}__NXbdef /ucachestatus{mark 0 0 0 0 0}__NXbdef /setucacheparams{cleartomark}__NXbdef /uappend{systemdict begin cvx exec end}__NXbdef /ueofill{gsave newpath uappend eofill grestore}__NXbdef /ufill{gsave newpath uappend fill grestore}__NXbdef /ustroke{ gsave newpath dup length 6 eq {exch uappend concat}{uappend}ifelse stroke grestore }__NXbdef /__NXustrokepathMatrix dup where {pop pop}{matrix def}ifelse /ustrokepath{ newpath dup length 6 eq{ exch uappend __NXustrokepathMatrix currentmatrix exch concat strokepath setmatrix }{uappend strokepath}ifelse } __NXbdef /upath{ [exch {/ucache cvx}if pathbbox /setbbox cvx {/moveto cvx}{/lineto cvx}{/curveto cvx}{/closepath cvx}pathforall]cvx } __NXbdef /setstrokeadjust{pop}__NXbdef /currentstrokeadjust{false}__NXbdef /selectfont{exch findfont exch dup type /arraytype eq {makefont}{scalefont}ifelse setfont}__NXbdef /_NXCombineArrays{ counttomark dup 2 add index dup length 3 -1 roll { 2 index length sub dup 4 1 roll 1 index exch 4 -1 roll putinterval exch }repeat pop pop pop }__NXbdef /flushgraphics{}def /setwindowtype{pop pop}def /currentwindowtype{pop 0}def /setalpha{pop}def /currentalpha{1.0}def /hidecursor{}def /obscurecursor{}def /revealcursor{}def /setcursor{4 {pop}repeat}bind def /showcursor{}def /NextStepEncoding where not{ /NextStepEncoding StandardEncoding 256 array copy def 0 [129/Agrave/Aacute/Acircumflex/Atilde/Adieresis/Aring/Ccedilla/Egrave /Eacute/Ecircumflex/Edieresis/Igrave/Iacute/Icircumflex/Idieresis /Eth/Ntilde/Ograve/Oacute/Ocircumflex/Otilde/Odieresis/Ugrave/Uacute /Ucircumflex/Udieresis/Yacute/Thorn/mu/multiply/divide/copyright 176/registered 181/brokenbar 190/logicalnot 192/onesuperior 201/twosuperior 204/threesuperior 209/plusminus/onequarter/onehalf/threequarters/agrave /aacute/acircumflex/atilde/adieresis/aring/ccedilla/egrave/eacute /ecircumflex/edieresis/igrave 226/iacute 228/icircumflex/idieresis/eth /ntilde 236/ograve/oacute/ocircumflex/otilde/odieresis 242/ugrave/uacute /ucircumflex 246/udieresis/yacute 252/thorn/ydieresis] {dup type /nametype eq {NextStepEncoding 2 index 2 index put pop 1 add}{exch pop}ifelse }forall pop /NextStepEncoding NextStepEncoding readonly def /_NXfstr 128 string dup 0 (_NX) putinterval def /_NXfindfont /findfont load def /findfont{ % Because we can never let NextStepEncoding get into % SharedFontDirectory, we cannot reencode a font to NextStepEncoding % if we are in shared mode. So if currentshared is true, % we call the normal findfont and return that /currentshared where {pop currentshared} {false} ifelse {_NXfindfont} {dup _NXfstr 3 125 getinterval cvs length 3 add _NXfstr 0 3 -1 roll getinterval cvn exch FontDirectory 2 index known {pop FontDirectory exch get} {_NXfindfont dup /Encoding get StandardEncoding eq { dup length dict exch {1 index /FID ne {2 index 3 1 roll put}{pop pop}ifelse}forall dup /Encoding NextStepEncoding put definefont }{exch pop} ifelse }ifelse }ifelse }bind def }{pop}ifelse /_NXImageString {/__NXImageString where{pop}{/__NXImageString 4000 string __NXdef}ifelse __NXImageString}__NXbdef /_NXDoImageOp{ 3 dict begin /parr 5 array def 1 index{dup}{1}ifelse /chans exch def chans 2 add 2 roll parr 0 chans getinterval astore pop 5 index 4 index mul 2 index{1 sub 8 idiv 1 add mul}{mul 1 sub 8 idiv 1 add}ifelse 4 index mul /totbytes exch def pop exch pop gsave matrix invertmatrix concat 0.5 setgray 0 0 4 2 roll rectfill grestore {0 1 chans 1 sub{parr exch get exec length totbytes exch sub /totbytes exch def}for totbytes 0 le{exit}if}loop end }__NXbdef /alphaimage{1 add _NXDoImageOp}def _NXLevel2{ /NXCalibratedRGBColorSpace where{pop}{ /NXCalibratedRGBColorSpace {mark /NXCalibratedRGB /ColorSpace findresource exch pop}stopped {cleartomark /NXCalibratedRGB[/CIEBasedABC 2 dict dup begin /MatrixLMN[.4124 .2126 .0193 .3576 .7152 .1192 .1805 .0722 .9505]def /WhitePoint[.9505 1 1.089] def end] /ColorSpace defineresource}if def}ifelse /nxsetrgbcolor{NXCalibratedRGBColorSpace setcolorspace setcolor}__NXbdef /nxsetgray{dup dup nxsetrgbcolor}__NXbdef /_NXCalibratedImage{exch{array astore dup length true}{false}ifelse 8 -1 roll{NXCalibratedRGBColorSpace setcolorspace}if 8 dict dup 9 1 roll begin /ImageType 1 def /MultipleDataSources exch def currentcolorspace 0 get /Indexed eq{pop /Decode[0 2 6 index exp 1 sub]def} {2 mul dup array /Decode exch def 1 sub 0 1 3 -1 roll{Decode exch dup 2 mod put}for}ifelse /DataSource exch def /ImageMatrix exch def /BitsPerComponent exch def /Height exch def /Width exch def end image}__NXbdef } { /setcmykcolor{ 1.0 exch sub dup dup 6 -1 roll sub dup 0 lt{pop 0}if 5 1 roll 4 -1 roll sub dup 0 lt{pop 0}if 3 1 roll exch sub dup 0 lt{pop 0}if setrgbcolor}__NXbdef /currentcmykcolor{currentrgbcolor 3{1.0 exch sub 3 1 roll}repeat 0}__NXbdef /colorimage{2 copy 3 ne or{_NXDoImageOp}{4 index dup 8 ne exch 4 ne and{_NXDoImageOp}{ pop pop save 6 1 roll 12 dict begin/Proc exch def/Res 0 string def /Alloc{2 index length mul 2 add dup 2 index load length gt{1.2 mul round cvi string def}{pop pop}ifelse}def 1 index 8 eq{/Unpack{.34 Alloc}def}{ /Wid 4 index 3 mul def exch pop 8 exch/Str1 0 string def/Lim Wid def /Unpack{.67 Alloc/Str1 2 Alloc 0 exch Lim exch {dup -4 bitshift 17 mul Str1 exch 4 index exch put 15 and 17 mul Str1 exch 3 index 1 add exch put 2 sub dup 0 le{0 lt Wid exch{exch 1 sub exch}if}if exch 2 add exch }forall/Lim exch def Str1 exch 0 exch getinterval }def }ifelse /Ops[{.3 mul add 1}{.59 mul add 2}{.11 mul add round cvi Res exch 2 index exch put 1 add 0.0 0}]def/Val 0.0 def/Phase 0 def {0 Val Phase Proc/Res Unpack{exch Ops exch get exec}forall/Phase exch def/Val exch def Res exch 0 exch getinterval} image end restore}ifelse}ifelse }__NXbdef /nxsetrgbcolor{setrgbcolor}__NXbdef /nxsetgray{setgray}__NXbdef /setpattern{pop .5 setgray}__NXbdef /_NXCalibratedImage{dup 1 eq {pop pop image}{colorimage}ifelse pop}__NXbdef } ifelse /_NXSetCMYKOrRGB where{pop}{ mark{systemdict /currentwindow get exec}stopped {{pop pop pop setcmykcolor}}{{nxsetrgbcolor pop pop pop pop}}ifelse /_NXSetCMYKOrRGB exch def cleartomark }ifelse %%EndProcSet gsave /__NXbasematrix matrix currentmatrix def grestore %%EndProlog %%BeginSetup /oval { translate scale newpath 0.5 0.5 0.5 0 360 arc closepath } def /line { moveto rlineto stroke } def /setup { setlinewidth setlinecap setlinejoin gsave } def /arrow { newpath moveto dup rotate -13 6 rlineto 4 -6 rlineto -4 -6 rlineto closepath gsave 0 setlinejoin stroke grestore fill neg rotate } def %%EndSetup 0 0 526 702 rectclip 74 362 452 340 rectclip 0 0 4.883721 setup 0 0 0 setup 0 nxsetgray 126 0 243 567 line grestore 0 0 4.883721 setup 0 nxsetgray 126 0 243 459 line grestore 0 0 0 setup 0 0 0 setup 0 nxsetgray 0 18 243 432 line grestore 0 0 0 setup 0 nxsetgray 45 0 243 441 line 180 243 441 arrow grestore 0 0 0 setup 0 nxsetgray 36 0 315 441 line 0 351 441 arrow grestore 0 0 0 setup gsave /Times-Bold findfont 16 scalefont [1 0 0 -1 0 0] makefont 69 exch defineuserobject 69 execuserobject setfont 0 nxsetgray [1 0 0 -1 0 881] concat 69 execuserobject setfont 0 nxsetgray 297 446 moveto (t) show grestore grestore grestore 0 0 0 setup 0 0 0 setup 0 nxsetgray 0 18 243 576 line grestore 0 0 0 setup 0 nxsetgray 27 0 243 585 line 180 243 585 arrow grestore 0 0 0 setup 0 nxsetgray 45 0 297 585 line 0 342 585 arrow grestore 0 0 0 setup gsave 69 execuserobject setfont 0 nxsetgray [1 0 0 -1 0 1167] concat 69 execuserobject setfont 0 nxsetgray 283 589 moveto (s) show grestore grestore grestore 0 0 0 setup gsave 69 execuserobject setfont 0 nxsetgray [1 0 0 -1 0 1024] concat 69 execuserobject setfont 0 nxsetgray 387 508 moveto (1 m) show grestore grestore 0 0 0 setup 0 nxsetgray 0 45 396 459 line 270.002655 396 459 arrow grestore 0 0 0 setup 0 nxsetgray 0 36 396 531 line 90.002655 396 567 arrow grestore 0 0 0 setup 0 nxsetgray 18 0 387 567 line grestore 0 0 0 setup 0 nxsetgray 18 0 387 459 line grestore 0 0 0 setup 0 nxsetgray 66 -53 203 514 line -38.766663 269 461 arrow grestore 0 0 0 setup gsave 69 execuserobject setfont 0 nxsetgray [1 0 0 -1 0 1034] concat 69 execuserobject setfont 0 nxsetgray 135 513 moveto (mass density) show 135 532 moveto (x\(t\) kg/m) show grestore grestore 0 0 0 setup 0 nxsetgray 63 -108 288 567 line grestore 0 0 0 setup 0 nxsetgray 135 -18 167 559 line -7.594868 302 541 arrow grestore 0 0 0 setup gsave 69 execuserobject setfont 0 nxsetgray [1 0 0 -1 0 1170] concat 69 execuserobject setfont 0 nxsetgray 99 581 moveto (Line of force from) show 99 600 moveto (pt \(s, 1\) to pt \(t, 0\)) show grestore grestore 0 0 0 setup 0 nxsetgray 72 18 171 441 line 14.036657 243 459 arrow grestore 0 0 0 setup gsave 69 execuserobject setfont 0 nxsetgray [1 0 0 -1 0 881] concat 69 execuserobject setfont 0 nxsetgray 141 446 moveto (bar) show grestore grestore 0 0 0 setup 0 nxsetgray 0 207 243 459 line 90.002655 243 666 arrow grestore 0 0 0 setup 0 nxsetgray 243 0 243 459 line 0 486 459 arrow grestore 0 0 0 setup gsave 69 execuserobject setfont 0 nxsetgray [1 0 0 -1 0 895] concat 69 execuserobject setfont 0 nxsetgray 492 453 moveto (x) show grestore grestore 0 0 0 setup gsave 69 execuserobject setfont 0 nxsetgray [1 0 0 -1 0 1339] concat 69 execuserobject setfont 0 nxsetgray 234 675 moveto (y) show grestore grestore 0 0 0 setup 0 nxsetgray 0 18 243 405 line grestore 0 0 0 setup 0 nxsetgray 0 18 369 405 line grestore 0 0 0 setup 0 nxsetgray 45 0 243 414 line 180 243 414 arrow grestore 0 0 0 setup 0 nxsetgray 45 0 324 414 line 0 369 414 arrow grestore 0 0 0 setup gsave 69 execuserobject setfont 0 nxsetgray [1 0 0 -1 0 808] concat 69 execuserobject setfont 0 nxsetgray 297 400 moveto (1 m) show grestore grestore grestore 0 setgray 0.333333 setgray gsave 0 0 526 702 rectclip [1 0 0 -1 0 702] concat grestore %%Trailer %%DocumentFonts: Times-Bold %End of Graphics MathPictureEnd :[font = subsection; inactive; Cclosed; preserveAspect; startGroup] Inverse Problem Approach - Problems 2 and 3 :[font = subsubsection; inactive; preserveAspect] Problem 2 :[font = subsubsection; inactive; preserveAspect] The task here is to determine the underlying mass density, x[t] kg/m^3, of the unit long mass at distance t meters from left hand edge - see diagram above. Observations about the vertical component of the gravitational force on a unit point mass of 1 kg which is resting 1 m above the bar where the particle is distance s m from left hand edge are known. Using an understanding of gravitational force, components of the vector of force, and addition of elements of force with the concept of integration determine the parameters of the model x(t) = a + b t which fit the data best. This then would be a reasonable conjecture as to what the underlying mass density really is from the observed surface data. :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] Suppose we have the following experimental data for the magnitude of the vertical component of the force of attraction between the bar and the unit point mass at (s, 1). The first coordinate is s and the second coordinate is the magnitude of the vertical component of the force. :[font = input; preserveAspect; endGroup] data = {{0, 5.8*10^-11}, {0.2, 6.8*10^-11}, {0.4, 7.4*10^-11}, {0.6, 7.6*10^-11}, {0.8, 7.1*10^-11}}; :[font = subsubsection; inactive; preserveAspect] Problem 3 :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] Suppose we have the following experimental data for the magnitude of the vertical component of the force of attraction between the bar and the unit point mass at (s, 1). The first coordinate is s and the second coordinate is the magnitude of the vertical component of the force. :[font = input; preserveAspect; endGroup] data = {{0, 5.1*10^-11}, {0.2, 5.9*10^-11}, {0.4, 6.4*10^-11}, {0.6, 6.3*10^-11}, {0.8, 5.8*10^-11}}; :[font = subsubsection; inactive; preserveAspect; endGroup; endGroup] Using an understanding of gravitational force, components of the vector of force, and addition of elements of force with the concept of integration determine the parameters of the model x(t) = a + b t + c t^2 which fit the data best. This then would be a reasonable conjecture as to what the underlying mass density really is from the observed surface data. :[font = section; inactive; Cclosed; preserveAspect; startGroup] KEYWORDS :[font = subsubsection; inactive; preserveAspect; endGroup] Newton's law of gravitational force of attraction between masses, vectors, vector projection, integration, function of integral, least squares, parameter estimation, inverse problem. :[font = section; inactive; Cclosed; preserveAspect; startGroup] TEACHER NOTES :[font = subsection; inactive; preserveAspect] ISSUES RELATED TO THE PROBLEM :[font = subsection; inactive; Cclosed; preserveAspect; startGroup] Prerequisites :[font = subsubsection; inactive; preserveAspect; endGroup] Concurrent or previous exposure to Newton's law of gravity, notion of force, force vector, projection of force vector, integration as a process of accumulation, least squares optimization for parameter estimation. :[font = subsection; inactive; Cclosed; preserveAspect; startGroup] Time allotment - time management :[font = subsubsection; inactive; preserveAspect; endGroup] With some guidance this problem can be given as a homework problem with about two days to give students chance to clarify their approach. Introducing it in class should take 5 minutes the first day with another 10 minutes the second day for clarifications. :[font = subsection; inactive; Cclosed; preserveAspect; startGroup] Expectations :[font = subsubsection; inactive; preserveAspect; endGroup] We expect students to see the difference between the forward and inverse problems. Expect students to have difficulty with finding and adding up the vertical components of the forces. :[font = subsection; inactive; Cclosed; preserveAspect; startGroup] Future payoffs :[font = subsubsection; inactive; preserveAspect; endGroup] Students will have seen an inverse problem and just how difficult it must be to solve the general inverse problem by doing some suggested modiffied inverse problems. :[font = subsection; inactive; Cclosed; preserveAspect; startGroup] Extensions :[font = subsubsection; inactive; preserveAspect] In the forward problem: One could offer several experimental values for vy[s], say at s = 0, .2, .4, .6, .8, and 1.0. Then after offering up several models of what the underlying mass densities x[t] could be one could ask which model for x[t] is most plausible given the data. One could again pose several experimental values for vy[s], say at s = 0, .2, .4, .6, .8, and 1.0 and offer up one model with an underlying set of parameters, say a quadratic x[t] = a t^2 + b t + c, and ask what are the best fitting parameters? :[font = subsubsection; inactive; preserveAspect; endGroup] In the inverse problem: Certainly the students could by now understand the very difficult nature of the inverse problem, namely, given the observations, what form and exactness could we offer for the underlying mass density distribution? Moreover, this approach does give some insight into the nature of the inverse problem in this case, while actually honing skills from the forward problem. :[font = subsection; inactive; Cclosed; preserveAspect; startGroup] References and Sources :[font = subsubsection; inactive; preserveAspect; endGroup; endGroup] 1. Groetsch, Charles W. 1993. Inverse Problems in the Mathematical Sciences. Weisbaden GERMANY: Vieweg Publishers. pp. 9-10. 2. Halliday, David, Robert Resnick, and Jearl Walker. 1993. Fundamentals of Physics, Fourth Edition. New York: John Wiley & Sons. Chapter 15: Gravitation. :[font = section; inactive; Cclosed; preserveAspect; startGroup] POSSIBLE SOLUTION(S) :[font = subsection; inactive; Cclosed; preserveAspect; startGroup] Problem 1 :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] G is the gravitational constant in N m^2/kg^2 or m^3/(kg s^2). :[font = input; preserveAspect; endGroup] G = 6.67 10^(-11); :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] Let x[t] be the mass density kg/m^3 of the one meter long mass at distance t meters from left hand edge. :[font = input; preserveAspect; endGroup] x[t_] := 100 + 90t^3 :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] Then vy[s] is the vertical component of the gravitational force on a particle of constant mass density, say 1, which is resting 1 m above the bar where s is distance in m from left hand edge. We obtain an integral by first determining the little element of force from (s, 1) to (t, 1) as DF = G 1 x(t) dt/(Sqrt[(s-t)^2 + 1])^2 . We then take the vertical component of this force, obtained by multiplying this magnitude DF by the cosine of the angle between the vertical and the line of direction for this force from (s, 1) to (t, 1), and integrate this over the entire line from (0,0) to (1, 0). :[font = input; preserveAspect; startGroup] vy[s_] = G Integrate[((s-t)^2 + 1)^(-3/2) x[t],{t,0,1}] :[font = output; output; inactive; preserveAspect; endGroup; endGroup] 6.67*10^-11*((-10*(18 - 10*s + 9*s^2 - 9*s^4))/ (1 + s^2)^(1/2) + (10*(37 - 55*s + 9*s^2 + 9*s^3 - 9*s^4))/ (2 - 2*s + s^2)^(1/2) + 270*s*ArcSinh[1 - s] + 270*s*ArcSinh[s]) ;[o] 2 4 -11 -10 (18 - 10 s + 9 s - 9 s ) 6.67 10 (----------------------------- + 2 Sqrt[1 + s ] 2 3 4 10 (37 - 55 s + 9 s + 9 s - 9 s ) ----------------------------------- + 2 Sqrt[2 - 2 s + s ] 270 s ArcSinh[1 - s] + 270 s ArcSinh[s]) :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] We plot a sketch of vy(s) noting its necessary asymmetry due to the non-uniform distribution of mass in our bar. :[font = input; preserveAspect; startGroup] Plot[vy[s],{s,0,1}, AxesLabel->{"Distance-s","MagVertForce"}] :[font = postscript; PostScript; formatAsPostScript; output; inactive; preserveAspect; pictureLeft = 100; pictureWidth = 385; pictureHeight = 238] %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 0.952381 -1.72865 3.20196e+008 [ [(0.2)] .21429 .19252 0 2 Msboxa [(0.4)] .40476 .19252 0 2 Msboxa [(0.6)] .59524 .19252 0 2 Msboxa [(0.8)] .78571 .19252 0 2 Msboxa [(1)] .97619 .19252 0 2 Msboxa [(Distance-s)] 1.025 .19252 -1 0 Msboxa [( -9)(5.5 10)] .01131 .03243 1 0 Msboxa [( -9)(5.75 10)] .01131 .11247 1 0 Msboxa [( -9)(6.25 10)] .01131 .27257 1 0 Msboxa [( -9)(6.5 10)] .01131 .35262 1 0 Msboxa [( -9)(6.75 10)] .01131 .43267 1 0 Msboxa [( -9)(7. 10)] .01131 .51272 1 0 Msboxa [( -9)(7.25 10)] .01131 .59277 1 0 Msboxa [(MagVertForce)] .02381 .61803 0 -4 Msboxa [ -0.001 -0.001 0 0 ] [ 1.001 .61903 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .21429 .19252 m .21429 .19877 L s P [(0.2)] .21429 .19252 0 2 Mshowa p .002 w .40476 .19252 m .40476 .19877 L s P [(0.4)] .40476 .19252 0 2 Mshowa p .002 w .59524 .19252 m .59524 .19877 L s P [(0.6)] .59524 .19252 0 2 Mshowa p .002 w .78571 .19252 m .78571 .19877 L s P [(0.8)] .78571 .19252 0 2 Mshowa p .002 w .97619 .19252 m .97619 .19877 L s P [(1)] .97619 .19252 0 2 Mshowa p .001 w .0619 .19252 m .0619 .19627 L s P p .001 w .1 .19252 m .1 .19627 L s P p .001 w .1381 .19252 m .1381 .19627 L s P p .001 w .17619 .19252 m .17619 .19627 L s P p .001 w .25238 .19252 m .25238 .19627 L s P p .001 w .29048 .19252 m .29048 .19627 L s P p .001 w .32857 .19252 m .32857 .19627 L s P p .001 w .36667 .19252 m .36667 .19627 L s P p .001 w .44286 .19252 m .44286 .19627 L s P p .001 w .48095 .19252 m .48095 .19627 L s P p .001 w .51905 .19252 m .51905 .19627 L s P p .001 w .55714 .19252 m .55714 .19627 L s P p .001 w .63333 .19252 m .63333 .19627 L s P p .001 w .67143 .19252 m .67143 .19627 L s P p .001 w .70952 .19252 m .70952 .19627 L s P p .001 w .74762 .19252 m .74762 .19627 L s P p .001 w .82381 .19252 m .82381 .19627 L s P p .001 w .8619 .19252 m .8619 .19627 L s P p .001 w .9 .19252 m .9 .19627 L s P p .001 w .9381 .19252 m .9381 .19627 L s P [(Distance-s)] 1.025 .19252 -1 0 Mshowa p .002 w 0 .19252 m 1 .19252 L s P p .002 w .02381 .03243 m .03006 .03243 L s P [( -9)(5.5 10)] .01131 .03243 1 0 Mshowa p .002 w .02381 .11247 m .03006 .11247 L s P [( -9)(5.75 10)] .01131 .11247 1 0 Mshowa p .002 w .02381 .27257 m .03006 .27257 L s P [( -9)(6.25 10)] .01131 .27257 1 0 Mshowa p .002 w .02381 .35262 m .03006 .35262 L s P [( -9)(6.5 10)] .01131 .35262 1 0 Mshowa p .002 w .02381 .43267 m .03006 .43267 L s P [( -9)(6.75 10)] .01131 .43267 1 0 Mshowa p .002 w .02381 .51272 m .03006 .51272 L s P [( -9)(7. 10)] .01131 .51272 1 0 Mshowa p .002 w .02381 .59277 m .03006 .59277 L s P [( -9)(7.25 10)] .01131 .59277 1 0 Mshowa p .001 w .02381 .04844 m .02756 .04844 L s P p .001 w .02381 .06445 m .02756 .06445 L s P p .001 w .02381 .08046 m .02756 .08046 L s P p .001 w .02381 .09646 m .02756 .09646 L s P p .001 w .02381 .12848 m .02756 .12848 L s P p .001 w .02381 .14449 m .02756 .14449 L s P p .001 w .02381 .1605 m .02756 .1605 L s P p .001 w .02381 .17651 m .02756 .17651 L s P p .001 w .02381 .20853 m .02756 .20853 L s P p .001 w .02381 .22454 m .02756 .22454 L s P p .001 w .02381 .24055 m .02756 .24055 L s P p .001 w .02381 .25656 m .02756 .25656 L s P p .001 w .02381 .28858 m .02756 .28858 L s P p .001 w .02381 .30459 m .02756 .30459 L s P p .001 w .02381 .3206 m .02756 .3206 L s P p .001 w .02381 .33661 m .02756 .33661 L s P p .001 w .02381 .36863 m .02756 .36863 L s P p .001 w .02381 .38464 m .02756 .38464 L s P p .001 w .02381 .40065 m .02756 .40065 L s P p .001 w .02381 .41666 m .02756 .41666 L s P p .001 w .02381 .44868 m .02756 .44868 L s P p .001 w .02381 .46469 m .02756 .46469 L s P p .001 w .02381 .4807 m .02756 .4807 L s P p .001 w .02381 .49671 m .02756 .49671 L s P p .001 w .02381 .52873 m .02756 .52873 L s P p .001 w .02381 .54474 m .02756 .54474 L s P p .001 w .02381 .56075 m .02756 .56075 L s P p .001 w .02381 .57676 m .02756 .57676 L s P p .001 w .02381 .01642 m .02756 .01642 L s P p .001 w .02381 .00041 m .02756 .00041 L s P p .001 w .02381 .60878 m .02756 .60878 L s P [(MagVertForce)] .02381 .61803 0 -4 Mshowa p .002 w .02381 0 m .02381 .61803 L s P P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath p p .004 w .02381 .01472 m .06349 .08499 L .10317 .15321 L .14286 .21876 L .18254 .28101 L .22222 .33932 L .2619 .39311 L .30159 .44182 L .34127 .48493 L .38095 .52195 L .42063 .55246 L .44048 .56516 L .46032 .57609 L .48016 .58523 L .5 .59253 L .51984 .59798 L .52976 .59999 L .53472 .60082 L .53968 .60153 L .54464 .60213 L .5496 .6026 L .55208 .60279 L .55456 .60296 L .55704 .60309 L .55952 .60319 L .56076 .60323 L .562 .60326 L .56324 .60329 L .56448 .60331 L .56572 .60332 L .56696 .60332 L .5682 .60331 L .56944 .6033 L .57068 .60328 L .57192 .60325 L .57316 .60322 L .5744 .60318 L .57688 .60307 L .57937 .60293 L .58433 .60256 L .58929 .60208 L .59921 .60074 L .60913 .59893 L .61905 .59663 L .63889 .59058 L .65873 .58262 L .69841 .56099 L .7381 .53192 L .77778 .49569 L .81746 .45268 L Mistroke .85714 .40335 L .89683 .34827 L .93651 .28806 L .97619 .22343 L Mfstroke P P % End of Graphics MathPictureEnd :[font = output; output; inactive; preserveAspect; endGroup; endGroup; endGroup] Graphics["<<>>"] ;[o] -Graphics- :[font = subsection; inactive; Cclosed; preserveAspect; startGroup] Problem 2 :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] We have the following experimental data for the magnitude of the vertical component of the force of attraction between the bar and the unit point mass at (s, 1). The first coordinate is s and the second coordinate is the magnitude of the vertical component of the force. :[font = input; preserveAspect; endGroup] data = {{0, 5.8*10^-11}, {0.2, 6.8*10^-11}, {0.4, 7.4*10^-11}, {0.6, 7.6*10^-11}, {0.8, 7.1*10^-11}}; :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] We plot the data: :[font = input; preserveAspect; startGroup] ysp = ListPlot[data,PlotStyle->PointSize[.03]] :[font = postscript; PostScript; formatAsPostScript; output; inactive; preserveAspect; pictureLeft = 100; pictureWidth = 300; pictureHeight = 185] %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 1.19048 -1.8819 3.27002e+010 [ [(0.2)] .2619 .08012 0 2 Msboxa [(0.4)] .5 .08012 0 2 Msboxa [(0.6)] .7381 .08012 0 2 Msboxa [(0.8)] .97619 .08012 0 2 Msboxa [( -11)(6.25 10)] .01131 .16187 1 0 Msboxa [( -11)(6.5 10)] .01131 .24362 1 0 Msboxa [( -11)(6.75 10)] .01131 .32537 1 0 Msboxa [( -11)(7. 10)] .01131 .40712 1 0 Msboxa [( -11)(7.25 10)] .01131 .48887 1 0 Msboxa [( -11)(7.5 10)] .01131 .57062 1 0 Msboxa [ -0.001 -0.001 0 0 ] [ 1.001 .61903 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .2619 .08012 m .2619 .08637 L s P [(0.2)] .2619 .08012 0 2 Mshowa p .002 w .5 .08012 m .5 .08637 L s P [(0.4)] .5 .08012 0 2 Mshowa p .002 w .7381 .08012 m .7381 .08637 L s P [(0.6)] .7381 .08012 0 2 Mshowa p .002 w .97619 .08012 m .97619 .08637 L s P [(0.8)] .97619 .08012 0 2 Mshowa p .001 w .07143 .08012 m .07143 .08387 L s P p .001 w .11905 .08012 m .11905 .08387 L s P p .001 w .16667 .08012 m .16667 .08387 L s P p .001 w .21429 .08012 m .21429 .08387 L s P p .001 w .30952 .08012 m .30952 .08387 L s P p .001 w .35714 .08012 m .35714 .08387 L s P p .001 w .40476 .08012 m .40476 .08387 L s P p .001 w .45238 .08012 m .45238 .08387 L s P p .001 w .54762 .08012 m .54762 .08387 L s P p .001 w .59524 .08012 m .59524 .08387 L s P p .001 w .64286 .08012 m .64286 .08387 L s P p .001 w .69048 .08012 m .69048 .08387 L s P p .001 w .78571 .08012 m .78571 .08387 L s P p .001 w .83333 .08012 m .83333 .08387 L s P p .001 w .88095 .08012 m .88095 .08387 L s P p .001 w .92857 .08012 m .92857 .08387 L s P p .002 w 0 .08012 m 1 .08012 L s P p .002 w .02381 .16187 m .03006 .16187 L s P [( -11)(6.25 10)] .01131 .16187 1 0 Mshowa p .002 w .02381 .24362 m .03006 .24362 L s P [( -11)(6.5 10)] .01131 .24362 1 0 Mshowa p .002 w .02381 .32537 m .03006 .32537 L s P [( -11)(6.75 10)] .01131 .32537 1 0 Mshowa p .002 w .02381 .40712 m .03006 .40712 L s P [( -11)(7. 10)] .01131 .40712 1 0 Mshowa p .002 w .02381 .48887 m .03006 .48887 L s P [( -11)(7.25 10)] .01131 .48887 1 0 Mshowa p .002 w .02381 .57062 m .03006 .57062 L s P [( -11)(7.5 10)] .01131 .57062 1 0 Mshowa p .001 w .02381 .01472 m .02756 .01472 L s P p .001 w .02381 .03107 m .02756 .03107 L s P p .001 w .02381 .04742 m .02756 .04742 L s P p .001 w .02381 .06377 m .02756 .06377 L s P p .001 w .02381 .09647 m .02756 .09647 L s P p .001 w .02381 .11282 m .02756 .11282 L s P p .001 w .02381 .12917 m .02756 .12917 L s P p .001 w .02381 .14552 m .02756 .14552 L s P p .001 w .02381 .17822 m .02756 .17822 L s P p .001 w .02381 .19457 m .02756 .19457 L s P p .001 w .02381 .21092 m .02756 .21092 L s P p .001 w .02381 .22727 m .02756 .22727 L s P p .001 w .02381 .25997 m .02756 .25997 L s P p .001 w .02381 .27632 m .02756 .27632 L s P p .001 w .02381 .29267 m .02756 .29267 L s P p .001 w .02381 .30902 m .02756 .30902 L s P p .001 w .02381 .34172 m .02756 .34172 L s P p .001 w .02381 .35807 m .02756 .35807 L s P p .001 w .02381 .37442 m .02756 .37442 L s P p .001 w .02381 .39077 m .02756 .39077 L s P p .001 w .02381 .42347 m .02756 .42347 L s P p .001 w .02381 .43982 m .02756 .43982 L s P p .001 w .02381 .45617 m .02756 .45617 L s P p .001 w .02381 .47252 m .02756 .47252 L s P p .001 w .02381 .50522 m .02756 .50522 L s P p .001 w .02381 .52157 m .02756 .52157 L s P p .001 w .02381 .53792 m .02756 .53792 L s P p .001 w .02381 .55427 m .02756 .55427 L s P p .001 w .02381 .58697 m .02756 .58697 L s P p .001 w .02381 .60332 m .02756 .60332 L s P p .002 w .02381 0 m .02381 .61803 L s P P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath p p .03 w .02381 .01472 Mdot .2619 .34172 Mdot .5 .53792 Mdot .7381 .60332 Mdot .97619 .43982 Mdot P P % End of Graphics MathPictureEnd :[font = output; output; inactive; preserveAspect; endGroup; endGroup] Graphics["<<>>"] ;[o] -Graphics- :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] We now attempt a model function x(t) to see if we can "determine" the parameters which fit the z(s) function best. :[font = input; preserveAspect; endGroup] x[t_] = a + b t; :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] We compute the vertical component of the force for this function x(t). :[font = input; preserveAspect; startGroup] yy[s_] = G Integrate[((s-t)^2 + 1)^(-3/2) x[t],{t,0,1}] :[font = output; output; inactive; preserveAspect; endGroup; endGroup] 6.67*10^-11*((a - b - a*s + b*s - b*s^2)/ (2 - 2*s + s^2)^(1/2) + (b + a*s + b*s^2)/(1 + s^2)^(1/2)) ;[o] 2 2 -11 a - b - a s + b s - b s b + a s + b s 6.67 10 (------------------------ + --------------) 2 2 Sqrt[2 - 2 s + s ] Sqrt[1 + s ] :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] We now determine the least square function of the parameters a and b in our function x(t) = a + b t. We seek to minimize ss(a, b) . :[font = input; preserveAspect; endGroup] ss[a_,b_] = Sum[(yy[data[[i,1]]] - data[[i,2]])^2, {i,1,Length[data]}]//N; :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] We use Mathematica's FindMinimum command to determine the values of a and b which fit the z(s) data best. :[font = input; preserveAspect; startGroup] sol = FindMinimum[ss[a,b],{a,1},{b,1}] :[font = output; output; inactive; preserveAspect; endGroup; endGroup] {3.665469379236343*10^-25, {a -> 1.018313115915117, b -> 0.5031310546632801}} ;[o] -25 {3.66547 10 , {a -> 1.01831, b -> 0.503131}} :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] We create our best fitting model function. :[font = input; preserveAspect; startGroup] xs[t_] = x[t]/.sol[[2]] :[font = output; output; inactive; preserveAspect; endGroup; endGroup] 1.018313115915117 + 0.5031310546632801*t ;[o] 1.01831 + 0.503131 t :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] So how good is our model for x(t), i.e., xs(t)? We put it into our z(s) function and plot the theoretical z's against the observed ones. :[font = input; preserveAspect; startGroup] z[s_] = G Integrate[((s-t)^2 + 1)^(-3/2) xs[t],{t,0,1}] :[font = output; output; inactive; preserveAspect; endGroup] 6.67*10^-11*(-((-0.5031310546632801 - 1.018313115915117*s - 0.5031310546632801*s^2) /(1 + s^2)^(1/2)) + (0.5151820612518372 - 0.5151820612518372*s - 0.5031310546632801*s^2)/(2 - 2*s + s^2)^(1/2)) ;[o] 2 -11 -0.503131 - 1.01831 s - 0.503131 s 6.67 10 (-(-----------------------------------) + 2 Sqrt[1 + s ] 2 0.515182 - 0.515182 s - 0.503131 s -----------------------------------) 2 Sqrt[2 - 2 s + s ] :[font = input; preserveAspect; startGroup] zplot = Plot[z[s],{s,0,1}] :[font = postscript; PostScript; formatAsPostScript; output; inactive; preserveAspect; pictureLeft = 100; pictureWidth = 300; pictureHeight = 185] %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 0.952381 -1.87457 3.26545e+010 [ [(0.2)] .21429 .0847 0 2 Msboxa [(0.4)] .40476 .0847 0 2 Msboxa [(0.6)] .59524 .0847 0 2 Msboxa [(0.8)] .78571 .0847 0 2 Msboxa [(1)] .97619 .0847 0 2 Msboxa [( -11)(5.75 10)] .01131 .00306 1 0 Msboxa [( -11)(6.25 10)] .01131 .16633 1 0 Msboxa [( -11)(6.5 10)] .01131 .24797 1 0 Msboxa [( -11)(6.75 10)] .01131 .32961 1 0 Msboxa [( -11)(7. 10)] .01131 .41124 1 0 Msboxa [( -11)(7.25 10)] .01131 .49288 1 0 Msboxa [( -11)(7.5 10)] .01131 .57451 1 0 Msboxa [ -0.001 -0.001 0 0 ] [ 1.001 .61903 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .21429 .0847 m .21429 .09095 L s P [(0.2)] .21429 .0847 0 2 Mshowa p .002 w .40476 .0847 m .40476 .09095 L s P [(0.4)] .40476 .0847 0 2 Mshowa p .002 w .59524 .0847 m .59524 .09095 L s P [(0.6)] .59524 .0847 0 2 Mshowa p .002 w .78571 .0847 m .78571 .09095 L s P [(0.8)] .78571 .0847 0 2 Mshowa p .002 w .97619 .0847 m .97619 .09095 L s P [(1)] .97619 .0847 0 2 Mshowa p .001 w .0619 .0847 m .0619 .08845 L s P p .001 w .1 .0847 m .1 .08845 L s P p .001 w .1381 .0847 m .1381 .08845 L s P p .001 w .17619 .0847 m .17619 .08845 L s P p .001 w .25238 .0847 m .25238 .08845 L s P p .001 w .29048 .0847 m .29048 .08845 L s P p .001 w .32857 .0847 m .32857 .08845 L s P p .001 w .36667 .0847 m .36667 .08845 L s P p .001 w .44286 .0847 m .44286 .08845 L s P p .001 w .48095 .0847 m .48095 .08845 L s P p .001 w .51905 .0847 m .51905 .08845 L s P p .001 w .55714 .0847 m .55714 .08845 L s P p .001 w .63333 .0847 m .63333 .08845 L s P p .001 w .67143 .0847 m .67143 .08845 L s P p .001 w .70952 .0847 m .70952 .08845 L s P p .001 w .74762 .0847 m .74762 .08845 L s P p .001 w .82381 .0847 m .82381 .08845 L s P p .001 w .8619 .0847 m .8619 .08845 L s P p .001 w .9 .0847 m .9 .08845 L s P p .001 w .9381 .0847 m .9381 .08845 L s P p .002 w 0 .0847 m 1 .0847 L s P p .002 w .02381 .00306 m .03006 .00306 L s P [( -11)(5.75 10)] .01131 .00306 1 0 Mshowa p .002 w .02381 .16633 m .03006 .16633 L s P [( -11)(6.25 10)] .01131 .16633 1 0 Mshowa p .002 w .02381 .24797 m .03006 .24797 L s P [( -11)(6.5 10)] .01131 .24797 1 0 Mshowa p .002 w .02381 .32961 m .03006 .32961 L s P [( -11)(6.75 10)] .01131 .32961 1 0 Mshowa p .002 w .02381 .41124 m .03006 .41124 L s P [( -11)(7. 10)] .01131 .41124 1 0 Mshowa p .002 w .02381 .49288 m .03006 .49288 L s P [( -11)(7.25 10)] .01131 .49288 1 0 Mshowa p .002 w .02381 .57451 m .03006 .57451 L s P [( -11)(7.5 10)] .01131 .57451 1 0 Mshowa p .001 w .02381 .01939 m .02756 .01939 L s P p .001 w .02381 .03572 m .02756 .03572 L s P p .001 w .02381 .05204 m .02756 .05204 L s P p .001 w .02381 .06837 m .02756 .06837 L s P p .001 w .02381 .10102 m .02756 .10102 L s P p .001 w .02381 .11735 m .02756 .11735 L s P p .001 w .02381 .13368 m .02756 .13368 L s P p .001 w .02381 .15001 m .02756 .15001 L s P p .001 w .02381 .18266 m .02756 .18266 L s P p .001 w .02381 .19899 m .02756 .19899 L s P p .001 w .02381 .21532 m .02756 .21532 L s P p .001 w .02381 .23164 m .02756 .23164 L s P p .001 w .02381 .2643 m .02756 .2643 L s P p .001 w .02381 .28062 m .02756 .28062 L s P p .001 w .02381 .29695 m .02756 .29695 L s P p .001 w .02381 .31328 m .02756 .31328 L s P p .001 w .02381 .34593 m .02756 .34593 L s P p .001 w .02381 .36226 m .02756 .36226 L s P p .001 w .02381 .37859 m .02756 .37859 L s P p .001 w .02381 .39491 m .02756 .39491 L s P p .001 w .02381 .42757 m .02756 .42757 L s P p .001 w .02381 .4439 m .02756 .4439 L s P p .001 w .02381 .46022 m .02756 .46022 L s P p .001 w .02381 .47655 m .02756 .47655 L s P p .001 w .02381 .50921 m .02756 .50921 L s P p .001 w .02381 .52553 m .02756 .52553 L s P p .001 w .02381 .54186 m .02756 .54186 L s P p .001 w .02381 .55819 m .02756 .55819 L s P p .001 w .02381 .59084 m .02756 .59084 L s P p .001 w .02381 .60717 m .02756 .60717 L s P p .002 w .02381 0 m .02381 .61803 L s P P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath p p .004 w .02381 .01472 m .06349 .08957 L .10317 .16188 L .14286 .23092 L .18254 .29598 L .22222 .35636 L .2619 .4114 L .30159 .46048 L .34127 .50304 L .38095 .53859 L .42063 .56671 L .44048 .57786 L .46032 .58703 L .48016 .59418 L .49008 .59699 L .5 .59929 L .50992 .60108 L .51488 .60177 L .51984 .60234 L .52232 .60258 L .5248 .60278 L .52728 .60295 L .52976 .60309 L .531 .60315 L .53224 .6032 L .53348 .60324 L .53472 .60327 L .53596 .60329 L .5372 .60331 L .53844 .60332 L .53968 .60332 L .54092 .60331 L .54216 .60329 L .5434 .60327 L .54464 .60324 L .54712 .60315 L .5496 .60303 L .55208 .60287 L .55456 .60268 L .55952 .60221 L .56944 .60087 L .57937 .59902 L .59921 .59374 L .61905 .58639 L .65873 .56555 L .69841 .5367 L .7381 .50014 L .77778 .45628 L .81746 .40561 L .85714 .34873 L Mistroke .89683 .28629 L .93651 .21901 L .97619 .14766 L Mfstroke P P % End of Graphics MathPictureEnd :[font = output; output; inactive; preserveAspect; endGroup; endGroup] Graphics["<<>>"] ;[o] -Graphics- :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] Not all that bad!!!!! :[font = input; preserveAspect; startGroup] Show[ysp,zplot] :[font = postscript; PostScript; formatAsPostScript; output; inactive; preserveAspect; pictureLeft = 100; pictureWidth = 300; pictureHeight = 185] %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 0.952381 -1.8623 3.24423e+010 [ [(0.2)] .21429 .08424 0 2 Msboxa [(0.4)] .40476 .08424 0 2 Msboxa [(0.6)] .59524 .08424 0 2 Msboxa [(0.8)] .78571 .08424 0 2 Msboxa [(1)] .97619 .08424 0 2 Msboxa [( -11)(5.75 10)] .01131 .00314 1 0 Msboxa [( -11)(6.25 10)] .01131 .16535 1 0 Msboxa [( -11)(6.5 10)] .01131 .24645 1 0 Msboxa [( -11)(6.75 10)] .01131 .32756 1 0 Msboxa [( -11)(7. 10)] .01131 .40867 1 0 Msboxa [( -11)(7.25 10)] .01131 .48977 1 0 Msboxa [( -11)(7.5 10)] .01131 .57088 1 0 Msboxa [ -0.001 -0.001 0 0 ] [ 1.001 .61903 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .21429 .08424 m .21429 .09049 L s P [(0.2)] .21429 .08424 0 2 Mshowa p .002 w .40476 .08424 m .40476 .09049 L s P [(0.4)] .40476 .08424 0 2 Mshowa p .002 w .59524 .08424 m .59524 .09049 L s P [(0.6)] .59524 .08424 0 2 Mshowa p .002 w .78571 .08424 m .78571 .09049 L s P [(0.8)] .78571 .08424 0 2 Mshowa p .002 w .97619 .08424 m .97619 .09049 L s P [(1)] .97619 .08424 0 2 Mshowa p .001 w .0619 .08424 m .0619 .08799 L s P p .001 w .1 .08424 m .1 .08799 L s P p .001 w .1381 .08424 m .1381 .08799 L s P p .001 w .17619 .08424 m .17619 .08799 L s P p .001 w .25238 .08424 m .25238 .08799 L s P p .001 w .29048 .08424 m .29048 .08799 L s P p .001 w .32857 .08424 m .32857 .08799 L s P p .001 w .36667 .08424 m .36667 .08799 L s P p .001 w .44286 .08424 m .44286 .08799 L s P p .001 w .48095 .08424 m .48095 .08799 L s P p .001 w .51905 .08424 m .51905 .08799 L s P p .001 w .55714 .08424 m .55714 .08799 L s P p .001 w .63333 .08424 m .63333 .08799 L s P p .001 w .67143 .08424 m .67143 .08799 L s P p .001 w .70952 .08424 m .70952 .08799 L s P p .001 w .74762 .08424 m .74762 .08799 L s P p .001 w .82381 .08424 m .82381 .08799 L s P p .001 w .8619 .08424 m .8619 .08799 L s P p .001 w .9 .08424 m .9 .08799 L s P p .001 w .9381 .08424 m .9381 .08799 L s P p .002 w 0 .08424 m 1 .08424 L s P p .002 w .02381 .00314 m .03006 .00314 L s P [( -11)(5.75 10)] .01131 .00314 1 0 Mshowa p .002 w .02381 .16535 m .03006 .16535 L s P [( -11)(6.25 10)] .01131 .16535 1 0 Mshowa p .002 w .02381 .24645 m .03006 .24645 L s P [( -11)(6.5 10)] .01131 .24645 1 0 Mshowa p .002 w .02381 .32756 m .03006 .32756 L s P [( -11)(6.75 10)] .01131 .32756 1 0 Mshowa p .002 w .02381 .40867 m .03006 .40867 L s P [( -11)(7. 10)] .01131 .40867 1 0 Mshowa p .002 w .02381 .48977 m .03006 .48977 L s P [( -11)(7.25 10)] .01131 .48977 1 0 Mshowa p .002 w .02381 .57088 m .03006 .57088 L s P [( -11)(7.5 10)] .01131 .57088 1 0 Mshowa p .001 w .02381 .01936 m .02756 .01936 L s P p .001 w .02381 .03558 m .02756 .03558 L s P p .001 w .02381 .0518 m .02756 .0518 L s P p .001 w .02381 .06802 m .02756 .06802 L s P p .001 w .02381 .10046 m .02756 .10046 L s P p .001 w .02381 .11668 m .02756 .11668 L s P p .001 w .02381 .13291 m .02756 .13291 L s P p .001 w .02381 .14913 m .02756 .14913 L s P p .001 w .02381 .18157 m .02756 .18157 L s P p .001 w .02381 .19779 m .02756 .19779 L s P p .001 w .02381 .21401 m .02756 .21401 L s P p .001 w .02381 .23023 m .02756 .23023 L s P p .001 w .02381 .26267 m .02756 .26267 L s P p .001 w .02381 .2789 m .02756 .2789 L s P p .001 w .02381 .29512 m .02756 .29512 L s P p .001 w .02381 .31134 m .02756 .31134 L s P p .001 w .02381 .34378 m .02756 .34378 L s P p .001 w .02381 .36 m .02756 .36 L s P p .001 w .02381 .37622 m .02756 .37622 L s P p .001 w .02381 .39244 m .02756 .39244 L s P p .001 w .02381 .42489 m .02756 .42489 L s P p .001 w .02381 .44111 m .02756 .44111 L s P p .001 w .02381 .45733 m .02756 .45733 L s P p .001 w .02381 .47355 m .02756 .47355 L s P p .001 w .02381 .50599 m .02756 .50599 L s P p .001 w .02381 .52221 m .02756 .52221 L s P p .001 w .02381 .53843 m .02756 .53843 L s P p .001 w .02381 .55466 m .02756 .55466 L s P p .001 w .02381 .5871 m .02756 .5871 L s P p .001 w .02381 .60332 m .02756 .60332 L s P p .002 w .02381 0 m .02381 .61803 L s P P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath p p p .03 w .02381 .01936 Mdot .21429 .34378 Mdot .40476 .53843 Mdot .59524 .60332 Mdot .78571 .44111 Mdot P P p p .004 w .02381 .01472 m .06349 .08908 L .10317 .16092 L .14286 .22952 L .18254 .29415 L .22222 .35414 L .2619 .40882 L .30159 .45758 L .34127 .49987 L .38095 .53519 L .42063 .56312 L .44048 .5742 L .46032 .58331 L .48016 .59042 L .49008 .59321 L .5 .59549 L .50992 .59727 L .51488 .59796 L .51984 .59852 L .52232 .59876 L .5248 .59896 L .52728 .59913 L .52976 .59927 L .531 .59932 L .53224 .59937 L .53348 .59941 L .53472 .59945 L .53596 .59947 L .5372 .59949 L .53844 .59949 L .53968 .59949 L .54092 .59949 L .54216 .59947 L .5434 .59944 L .54464 .59941 L .54712 .59932 L .5496 .5992 L .55208 .59905 L .55456 .59886 L .55952 .59839 L .56944 .59706 L .57937 .59522 L .59921 .58998 L .61905 .58268 L .65873 .56197 L .69841 .53331 L .7381 .49699 L .77778 .45341 L .81746 .40307 L .85714 .34656 L Mistroke .89683 .28452 L .93651 .21768 L .97619 .1468 L Mfstroke P P P % End of Graphics MathPictureEnd :[font = output; output; inactive; preserveAspect; endGroup; endGroup; endGroup] Graphics["<<>>"] ;[o] -Graphics- :[font = subsection; inactive; Cclosed; preserveAspect; startGroup] Problem 3 :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] Suppose we have the following experimental data for the magnitude of the vertical component of the force of attraction between the bar and the unit point mass at (s, 1). The first coordinate is s and the second cooredinate is the magnitude of the vertical component of the force. :[font = input; preserveAspect; endGroup] data = {{0, 5.1*10^-11}, {0.2, 5.9*10^-11}, {0.4, 6.4*10^-11}, {0.6, 6.3*10^-11}, {0.8, 5.8*10^-11}}; :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] Then ys[s] is the vertical component of the gravitational force on a particle resting 1 m above the bar where s is distance in m from left hand edge. :[font = input; preserveAspect; endGroup] ys[s_] := G Integrate[((s-t)^2 + 1)^(-3/2) x[t],{t,0,1}] :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] We plot the data: :[font = input; preserveAspect; startGroup] ysp = ListPlot[data,PlotStyle->PointSize[.03]] :[font = postscript; PostScript; formatAsPostScript; output; inactive; preserveAspect; pictureLeft = 100; pictureWidth = 300; pictureHeight = 185] %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 1.19048 -2.29442 4.52772e+010 [ [(0.2)] .2619 .42221 0 2 Msboxa [(0.4)] .5 .42221 0 2 Msboxa [(0.6)] .7381 .42221 0 2 Msboxa [(0.8)] .97619 .42221 0 2 Msboxa [( -11)(5.2 10)] .01131 .05999 1 0 Msboxa [( -11)(5.4 10)] .01131 .15055 1 0 Msboxa [( -11)(5.6 10)] .01131 .2411 1 0 Msboxa [( -11)(5.8 10)] .01131 .33166 1 0 Msboxa [( -11)(6.2 10)] .01131 .51276 1 0 Msboxa [( -11)(6.4 10)] .01131 .60332 1 0 Msboxa [ -0.001 -0.001 0 0 ] [ 1.001 .61903 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .2619 .42221 m .2619 .42846 L s P [(0.2)] .2619 .42221 0 2 Mshowa p .002 w .5 .42221 m .5 .42846 L s P [(0.4)] .5 .42221 0 2 Mshowa p .002 w .7381 .42221 m .7381 .42846 L s P [(0.6)] .7381 .42221 0 2 Mshowa p .002 w .97619 .42221 m .97619 .42846 L s P [(0.8)] .97619 .42221 0 2 Mshowa p .001 w .07143 .42221 m .07143 .42596 L s P p .001 w .11905 .42221 m .11905 .42596 L s P p .001 w .16667 .42221 m .16667 .42596 L s P p .001 w .21429 .42221 m .21429 .42596 L s P p .001 w .30952 .42221 m .30952 .42596 L s P p .001 w .35714 .42221 m .35714 .42596 L s P p .001 w .40476 .42221 m .40476 .42596 L s P p .001 w .45238 .42221 m .45238 .42596 L s P p .001 w .54762 .42221 m .54762 .42596 L s P p .001 w .59524 .42221 m .59524 .42596 L s P p .001 w .64286 .42221 m .64286 .42596 L s P p .001 w .69048 .42221 m .69048 .42596 L s P p .001 w .78571 .42221 m .78571 .42596 L s P p .001 w .83333 .42221 m .83333 .42596 L s P p .001 w .88095 .42221 m .88095 .42596 L s P p .001 w .92857 .42221 m .92857 .42596 L s P p .002 w 0 .42221 m 1 .42221 L s P p .002 w .02381 .05999 m .03006 .05999 L s P [( -11)(5.2 10)] .01131 .05999 1 0 Mshowa p .002 w .02381 .15055 m .03006 .15055 L s P [( -11)(5.4 10)] .01131 .15055 1 0 Mshowa p .002 w .02381 .2411 m .03006 .2411 L s P [( -11)(5.6 10)] .01131 .2411 1 0 Mshowa p .002 w .02381 .33166 m .03006 .33166 L s P [( -11)(5.8 10)] .01131 .33166 1 0 Mshowa p .002 w .02381 .51276 m .03006 .51276 L s P [( -11)(6.2 10)] .01131 .51276 1 0 Mshowa p .002 w .02381 .60332 m .03006 .60332 L s P [( -11)(6.4 10)] .01131 .60332 1 0 Mshowa p .001 w .02381 .0781 m .02756 .0781 L s P p .001 w .02381 .09621 m .02756 .09621 L s P p .001 w .02381 .11432 m .02756 .11432 L s P p .001 w .02381 .13244 m .02756 .13244 L s P p .001 w .02381 .16866 m .02756 .16866 L s P p .001 w .02381 .18677 m .02756 .18677 L s P p .001 w .02381 .20488 m .02756 .20488 L s P p .001 w .02381 .22299 m .02756 .22299 L s P p .001 w .02381 .25921 m .02756 .25921 L s P p .001 w .02381 .27732 m .02756 .27732 L s P p .001 w .02381 .29543 m .02756 .29543 L s P p .001 w .02381 .31354 m .02756 .31354 L s P p .001 w .02381 .34977 m .02756 .34977 L s P p .001 w .02381 .36788 m .02756 .36788 L s P p .001 w .02381 .38599 m .02756 .38599 L s P p .001 w .02381 .4041 m .02756 .4041 L s P p .001 w .02381 .44032 m .02756 .44032 L s P p .001 w .02381 .45843 m .02756 .45843 L s P p .001 w .02381 .47654 m .02756 .47654 L s P p .001 w .02381 .49465 m .02756 .49465 L s P p .001 w .02381 .53088 m .02756 .53088 L s P p .001 w .02381 .54899 m .02756 .54899 L s P p .001 w .02381 .5671 m .02756 .5671 L s P p .001 w .02381 .58521 m .02756 .58521 L s P p .001 w .02381 .04188 m .02756 .04188 L s P p .001 w .02381 .02377 m .02756 .02377 L s P p .001 w .02381 .00566 m .02756 .00566 L s P p .002 w .02381 0 m .02381 .61803 L s P P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath p p .03 w .02381 .01472 Mdot .2619 .37693 Mdot .5 .60332 Mdot .7381 .55804 Mdot .97619 .33166 Mdot P P % End of Graphics MathPictureEnd :[font = output; output; inactive; preserveAspect; endGroup; endGroup] Graphics["<<>>"] ;[o] -Graphics- :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] We now attempt a model function x(t) to see if we can "determine" the parameters which fit the ys(s) function best. :[font = input; preserveAspect; startGroup] x[t_] = a + b t + c t^2 :[font = output; output; inactive; preserveAspect; endGroup; endGroup] a + b*t + c*t^2 ;[o] 2 a + b t + c t :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] We compute the vertical component for this function x(t). :[font = input; preserveAspect; startGroup] zz[s_] = G Integrate[((s-t)^2 + 1)^(-3/2) x[t],{t,0,1}] :[font = output; output; inactive; preserveAspect; endGroup; endGroup] 6.67*10^-11*((a - b - c - a*s + b*s - c*s - b*s^2 + c*s^2 - c*s^3)/(2 - 2*s + s^2)^(1/2) + (b + a*s + c*s + b*s^2 + c*s^3)/(1 + s^2)^(1/2) + c*ArcSinh[1 - s] + c*ArcSinh[s]) ;[o] -11 6.67 10 ( 2 2 3 a - b - c - a s + b s - c s - b s + c s - c s ------------------------------------------------\ 2 Sqrt[2 - 2 s + s ] 2 3 b + a s + c s + b s + c s + --------------------------- + 2 Sqrt[1 + s ] c ArcSinh[1 - s] + c ArcSinh[s]) :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] We now determine the least square function of the parameters a and b in our function x(t) = a + b t + c t^2. We seek to minimize ss(a, b,c) . :[font = input; preserveAspect; endGroup] ss[a_,b_,c_] = Sum[(zz[data[[i,1]]] - data[[i,2]])^2, {i,1,Length[data]}]//N; :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] We use Mathematica's FindMinimum command to determine the values of a and b which fit the z(s) data best. :[font = input; preserveAspect; startGroup] sol = FindMinimum[ss[a,b,c],{a,3},{b,2},{c,4}] :[font = output; output; inactive; preserveAspect; endGroup; endGroup] {1.824807083322384*10^-25, {a -> 0.839821527694304, b -> 1.732541874516921, c -> -1.936505133064068}} ;[o] -25 {1.82481 10 , {a -> 0.839822, b -> 1.73254, c -> -1.93651}} :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] We create our best fitting model function. :[font = input; preserveAspect; startGroup] xb[t_] = x[t]/.sol[[2]] :[font = output; output; inactive; preserveAspect; endGroup; endGroup] 0.839821527694304 + 1.732541874516921*t - 1.936505133064068*t^2 ;[o] 2 0.839822 + 1.73254 t - 1.93651 t :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] So how good is our model for x(t), i.e. xb(t)? We put it into our z(s) function and plot the theoretical z's against the observed ones. :[font = input; preserveAspect; startGroup] z[s_] = G Integrate[((s-t)^2 + 1)^(-3/2) xb[t],{t,0,1}] :[font = output; output; inactive; preserveAspect; endGroup] 6.67*10^-11*((1.043784786241451 + 2.829225479886685*s - 3.669047007580988*s^2 + 1.936505133064068*s^3)/(2 - 2*s + s^2)^(1/2) - (-1.732541874516921 + 1.096683605369764*s - 1.732541874516921*s^2 + 1.936505133064068*s^3)/ (1 + s^2)^(1/2) - 1.936505133064068*ArcSinh[1 - s] - 1.936505133064068*ArcSinh[s]) ;[o] -11 6.67 10 ( 2 3 1.04378 + 2.82923 s - 3.66905 s + 1.93651 s --------------------------------------------- - 2 Sqrt[2 - 2 s + s ] 2 3 -1.73254 + 1.09668 s - 1.73254 s + 1.93651 s ---------------------------------------------- - 2 Sqrt[1 + s ] 1.93651 ArcSinh[1 - s] - 1.93651 ArcSinh[s]) :[font = input; preserveAspect; startGroup] zplot = Plot[z[s],{s,0,1}] :[font = postscript; PostScript; formatAsPostScript; output; inactive; preserveAspect; pictureLeft = 100; pictureWidth = 300; pictureHeight = 185] %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 0.952381 -1.94007 3.96535e+010 [ [(0.2)] .21429 .04261 0 2 Msboxa [(0.4)] .40476 .04261 0 2 Msboxa [(0.6)] .59524 .04261 0 2 Msboxa [(0.8)] .78571 .04261 0 2 Msboxa [(1)] .97619 .04261 0 2 Msboxa [( -11)(5.2 10)] .01131 .12191 1 0 Msboxa [( -11)(5.4 10)] .01131 .20122 1 0 Msboxa [( -11)(5.6 10)] .01131 .28053 1 0 Msboxa [( -11)(5.8 10)] .01131 .35983 1 0 Msboxa [( -11)(6. 10)] .01131 .43914 1 0 Msboxa [( -11)(6.2 10)] .01131 .51845 1 0 Msboxa [( -11)(6.4 10)] .01131 .59776 1 0 Msboxa [ -0.001 -0.001 0 0 ] [ 1.001 .61903 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .21429 .04261 m .21429 .04886 L s P [(0.2)] .21429 .04261 0 2 Mshowa p .002 w .40476 .04261 m .40476 .04886 L s P [(0.4)] .40476 .04261 0 2 Mshowa p .002 w .59524 .04261 m .59524 .04886 L s P [(0.6)] .59524 .04261 0 2 Mshowa p .002 w .78571 .04261 m .78571 .04886 L s P [(0.8)] .78571 .04261 0 2 Mshowa p .002 w .97619 .04261 m .97619 .04886 L s P [(1)] .97619 .04261 0 2 Mshowa p .001 w .0619 .04261 m .0619 .04636 L s P p .001 w .1 .04261 m .1 .04636 L s P p .001 w .1381 .04261 m .1381 .04636 L s P p .001 w .17619 .04261 m .17619 .04636 L s P p .001 w .25238 .04261 m .25238 .04636 L s P p .001 w .29048 .04261 m .29048 .04636 L s P p .001 w .32857 .04261 m .32857 .04636 L s P p .001 w .36667 .04261 m .36667 .04636 L s P p .001 w .44286 .04261 m .44286 .04636 L s P p .001 w .48095 .04261 m .48095 .04636 L s P p .001 w .51905 .04261 m .51905 .04636 L s P p .001 w .55714 .04261 m .55714 .04636 L s P p .001 w .63333 .04261 m .63333 .04636 L s P p .001 w .67143 .04261 m .67143 .04636 L s P p .001 w .70952 .04261 m .70952 .04636 L s P p .001 w .74762 .04261 m .74762 .04636 L s P p .001 w .82381 .04261 m .82381 .04636 L s P p .001 w .8619 .04261 m .8619 .04636 L s P p .001 w .9 .04261 m .9 .04636 L s P p .001 w .9381 .04261 m .9381 .04636 L s P p .002 w 0 .04261 m 1 .04261 L s P p .002 w .02381 .12191 m .03006 .12191 L s P [( -11)(5.2 10)] .01131 .12191 1 0 Mshowa p .002 w .02381 .20122 m .03006 .20122 L s P [( -11)(5.4 10)] .01131 .20122 1 0 Mshowa p .002 w .02381 .28053 m .03006 .28053 L s P [( -11)(5.6 10)] .01131 .28053 1 0 Mshowa p .002 w .02381 .35983 m .03006 .35983 L s P [( -11)(5.8 10)] .01131 .35983 1 0 Mshowa p .002 w .02381 .43914 m .03006 .43914 L s P [( -11)(6. 10)] .01131 .43914 1 0 Mshowa p .002 w .02381 .51845 m .03006 .51845 L s P [( -11)(6.2 10)] .01131 .51845 1 0 Mshowa p .002 w .02381 .59776 m .03006 .59776 L s P [( -11)(6.4 10)] .01131 .59776 1 0 Mshowa p .001 w .02381 .05847 m .02756 .05847 L s P p .001 w .02381 .07433 m .02756 .07433 L s P p .001 w .02381 .09019 m .02756 .09019 L s P p .001 w .02381 .10605 m .02756 .10605 L s P p .001 w .02381 .13777 m .02756 .13777 L s P p .001 w .02381 .15364 m .02756 .15364 L s P p .001 w .02381 .1695 m .02756 .1695 L s P p .001 w .02381 .18536 m .02756 .18536 L s P p .001 w .02381 .21708 m .02756 .21708 L s P p .001 w .02381 .23294 m .02756 .23294 L s P p .001 w .02381 .2488 m .02756 .2488 L s P p .001 w .02381 .26467 m .02756 .26467 L s P p .001 w .02381 .29639 m .02756 .29639 L s P p .001 w .02381 .31225 m .02756 .31225 L s P p .001 w .02381 .32811 m .02756 .32811 L s P p .001 w .02381 .34397 m .02756 .34397 L s P p .001 w .02381 .3757 m .02756 .3757 L s P p .001 w .02381 .39156 m .02756 .39156 L s P p .001 w .02381 .40742 m .02756 .40742 L s P p .001 w .02381 .42328 m .02756 .42328 L s P p .001 w .02381 .455 m .02756 .455 L s P p .001 w .02381 .47086 m .02756 .47086 L s P p .001 w .02381 .48673 m .02756 .48673 L s P p .001 w .02381 .50259 m .02756 .50259 L s P p .001 w .02381 .53431 m .02756 .53431 L s P p .001 w .02381 .55017 m .02756 .55017 L s P p .001 w .02381 .56603 m .02756 .56603 L s P p .001 w .02381 .58189 m .02756 .58189 L s P p .001 w .02381 .02675 m .02756 .02675 L s P p .001 w .02381 .01088 m .02756 .01088 L s P p .001 w .02381 .61362 m .02756 .61362 L s P p .002 w .02381 0 m .02381 .61803 L s P P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath p p .004 w .02381 .08016 m .06349 .15645 L .10317 .2292 L .14286 .29757 L .18254 .36075 L .22222 .41796 L .2619 .46849 L .30159 .51171 L .34127 .54705 L .36111 .56163 L .38095 .57408 L .40079 .58437 L .42063 .59245 L .44048 .59832 L .4504 .60041 L .45536 .60124 L .46032 .60194 L .46528 .60249 L .47024 .6029 L .47272 .60306 L .47396 .60312 L .4752 .60318 L .47644 .60322 L .47768 .60326 L .47892 .60329 L .48016 .60331 L .4814 .60332 L .48264 .60332 L .48388 .60331 L .48512 .6033 L .48636 .60327 L .4876 .60324 L .49008 .60315 L .49256 .60302 L .49504 .60285 L .5 .60242 L .50496 .60185 L .50992 .60114 L .51984 .59929 L .52976 .59688 L .53968 .59392 L .55952 .58633 L .57937 .57655 L .61905 .55055 L .65873 .51628 L .69841 .47419 L .7381 .42481 L .77778 .36877 L .81746 .30678 L .85714 .23958 L Mistroke .89683 .16796 L .93651 .09273 L .97619 .01472 L Mfstroke P P % End of Graphics MathPictureEnd :[font = output; output; inactive; preserveAspect; endGroup; endGroup] Graphics["<<>>"] ;[o] -Graphics- :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] Not all that bad!!!!! :[font = input; preserveAspect; startGroup] Show[ysp,zplot] :[font = postscript; PostScript; formatAsPostScript; output; inactive; preserveAspect; pictureLeft = 100; pictureWidth = 300; pictureHeight = 185] %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 0.952381 -1.94007 3.96535e+010 [ [(0.2)] .21429 .04261 0 2 Msboxa [(0.4)] .40476 .04261 0 2 Msboxa [(0.6)] .59524 .04261 0 2 Msboxa [(0.8)] .78571 .04261 0 2 Msboxa [(1)] .97619 .04261 0 2 Msboxa [( -11)(5.2 10)] .01131 .12191 1 0 Msboxa [( -11)(5.4 10)] .01131 .20122 1 0 Msboxa [( -11)(5.6 10)] .01131 .28053 1 0 Msboxa [( -11)(5.8 10)] .01131 .35983 1 0 Msboxa [( -11)(6. 10)] .01131 .43914 1 0 Msboxa [( -11)(6.2 10)] .01131 .51845 1 0 Msboxa [( -11)(6.4 10)] .01131 .59776 1 0 Msboxa [ -0.001 -0.001 0 0 ] [ 1.001 .61903 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .21429 .04261 m .21429 .04886 L s P [(0.2)] .21429 .04261 0 2 Mshowa p .002 w .40476 .04261 m .40476 .04886 L s P [(0.4)] .40476 .04261 0 2 Mshowa p .002 w .59524 .04261 m .59524 .04886 L s P [(0.6)] .59524 .04261 0 2 Mshowa p .002 w .78571 .04261 m .78571 .04886 L s P [(0.8)] .78571 .04261 0 2 Mshowa p .002 w .97619 .04261 m .97619 .04886 L s P [(1)] .97619 .04261 0 2 Mshowa p .001 w .0619 .04261 m .0619 .04636 L s P p .001 w .1 .04261 m .1 .04636 L s P p .001 w .1381 .04261 m .1381 .04636 L s P p .001 w .17619 .04261 m .17619 .04636 L s P p .001 w .25238 .04261 m .25238 .04636 L s P p .001 w .29048 .04261 m .29048 .04636 L s P p .001 w .32857 .04261 m .32857 .04636 L s P p .001 w .36667 .04261 m .36667 .04636 L s P p .001 w .44286 .04261 m .44286 .04636 L s P p .001 w .48095 .04261 m .48095 .04636 L s P p .001 w .51905 .04261 m .51905 .04636 L s P p .001 w .55714 .04261 m .55714 .04636 L s P p .001 w .63333 .04261 m .63333 .04636 L s P p .001 w .67143 .04261 m .67143 .04636 L s P p .001 w .70952 .04261 m .70952 .04636 L s P p .001 w .74762 .04261 m .74762 .04636 L s P p .001 w .82381 .04261 m .82381 .04636 L s P p .001 w .8619 .04261 m .8619 .04636 L s P p .001 w .9 .04261 m .9 .04636 L s P p .001 w .9381 .04261 m .9381 .04636 L s P p .002 w 0 .04261 m 1 .04261 L s P p .002 w .02381 .12191 m .03006 .12191 L s P [( -11)(5.2 10)] .01131 .12191 1 0 Mshowa p .002 w .02381 .20122 m .03006 .20122 L s P [( -11)(5.4 10)] .01131 .20122 1 0 Mshowa p .002 w .02381 .28053 m .03006 .28053 L s P [( -11)(5.6 10)] .01131 .28053 1 0 Mshowa p .002 w .02381 .35983 m .03006 .35983 L s P [( -11)(5.8 10)] .01131 .35983 1 0 Mshowa p .002 w .02381 .43914 m .03006 .43914 L s P [( -11)(6. 10)] .01131 .43914 1 0 Mshowa p .002 w .02381 .51845 m .03006 .51845 L s P [( -11)(6.2 10)] .01131 .51845 1 0 Mshowa p .002 w .02381 .59776 m .03006 .59776 L s P [( -11)(6.4 10)] .01131 .59776 1 0 Mshowa p .001 w .02381 .05847 m .02756 .05847 L s P p .001 w .02381 .07433 m .02756 .07433 L s P p .001 w .02381 .09019 m .02756 .09019 L s P p .001 w .02381 .10605 m .02756 .10605 L s P p .001 w .02381 .13777 m .02756 .13777 L s P p .001 w .02381 .15364 m .02756 .15364 L s P p .001 w .02381 .1695 m .02756 .1695 L s P p .001 w .02381 .18536 m .02756 .18536 L s P p .001 w .02381 .21708 m .02756 .21708 L s P p .001 w .02381 .23294 m .02756 .23294 L s P p .001 w .02381 .2488 m .02756 .2488 L s P p .001 w .02381 .26467 m .02756 .26467 L s P p .001 w .02381 .29639 m .02756 .29639 L s P p .001 w .02381 .31225 m .02756 .31225 L s P p .001 w .02381 .32811 m .02756 .32811 L s P p .001 w .02381 .34397 m .02756 .34397 L s P p .001 w .02381 .3757 m .02756 .3757 L s P p .001 w .02381 .39156 m .02756 .39156 L s P p .001 w .02381 .40742 m .02756 .40742 L s P p .001 w .02381 .42328 m .02756 .42328 L s P p .001 w .02381 .455 m .02756 .455 L s P p .001 w .02381 .47086 m .02756 .47086 L s P p .001 w .02381 .48673 m .02756 .48673 L s P p .001 w .02381 .50259 m .02756 .50259 L s P p .001 w .02381 .53431 m .02756 .53431 L s P p .001 w .02381 .55017 m .02756 .55017 L s P p .001 w .02381 .56603 m .02756 .56603 L s P p .001 w .02381 .58189 m .02756 .58189 L s P p .001 w .02381 .02675 m .02756 .02675 L s P p .001 w .02381 .01088 m .02756 .01088 L s P p .001 w .02381 .61362 m .02756 .61362 L s P p .002 w .02381 0 m .02381 .61803 L s P P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath p p p .03 w .02381 .08226 Mdot .21429 .39949 Mdot .40476 .59776 Mdot .59524 .5581 Mdot .78571 .35983 Mdot P P p p .004 w .02381 .08016 m .06349 .15645 L .10317 .2292 L .14286 .29757 L .18254 .36075 L .22222 .41796 L .2619 .46849 L .30159 .51171 L .34127 .54705 L .36111 .56163 L .38095 .57408 L .40079 .58437 L .42063 .59245 L .44048 .59832 L .4504 .60041 L .45536 .60124 L .46032 .60194 L .46528 .60249 L .47024 .6029 L .47272 .60306 L .47396 .60312 L .4752 .60318 L .47644 .60322 L .47768 .60326 L .47892 .60329 L .48016 .60331 L .4814 .60332 L .48264 .60332 L .48388 .60331 L .48512 .6033 L .48636 .60327 L .4876 .60324 L .49008 .60315 L .49256 .60302 L .49504 .60285 L .5 .60242 L .50496 .60185 L .50992 .60114 L .51984 .59929 L .52976 .59688 L .53968 .59392 L .55952 .58633 L .57937 .57655 L .61905 .55055 L .65873 .51628 L .69841 .47419 L .7381 .42481 L .77778 .36877 L .81746 .30678 L .85714 .23958 L Mistroke .89683 .16796 L .93651 .09273 L .97619 .01472 L Mfstroke P P P % End of Graphics MathPictureEnd :[font = output; output; inactive; preserveAspect; endGroup; endGroup; endGroup; endGroup] Graphics["<<>>"] ;[o] -Graphics- :[font = section; inactive; Cclosed; preserveAspect; startGroup] ISSUES IN SOLUTION :[font = subsection; inactive; preserveAspect; endGroup; endGroup] Students do have difficulty with finding and adding up the vertical components of the forces. There are a number of different concepts running through this problem and students need to keep straight where they are in the analysis. Summarizing what they have accomplished at each stage is important before they proceed. ^*)