Flux Invariants for Shape
|
|
- Sarah Todd
- 5 years ago
- Views:
Transcription
1 Flux Invariants for Shape avel Dimitrov James N. Damon Kaleem Siddiqi School of Computer Science Department of Mathematics McGill University University of North Carolina, Chapel Hill Montreal, Canada Chapel Hill, North Carolina, USA Abstract We consider the average outward flux through a Jordan curve of the gradient vector field of the Euclidean distance function to the boundary of a 2D shape. Using an alternate form of the divergence theorem, we show that in the it as the area of the region enclosed by such a curve shrinks to zero, this measure has very different behaviours at medial points than at non-medial ones, providing a theoretical justification for its use in the Hamilton-Jacobi skeletonization algorithm of [7]. We then specialize to the case of shrinking circular neighborhoods and show that the average outward flux measure also reveals the object angle at skeletal points. Hence, formulae for obtaining the boundary curves, their curvatures, and other geometric quantities of interest, can be written in terms of the average outward flux it values at skeletal points. Thus this measure can be viewed as a Euclidean invariant for shape description: it can be used to both detect the skeleton from the Euclidean distance function, as well as to explicitly reconstruct the boundary from it. We illustrate our results with several numerical simulations. 1. Introduction In recent years there has been a resurgence of interest in representations based on the Blum skeleton, and extensions to 3D, for shape representation and analysis have been developed [9, 5, 8]. A number of abstractions related to the skeleton have been proposed, and algorithms for a diverse range of computer vision and medical imaging applications have been developed. These include shape segmentation, shape matching, view-based object recognition and indexing. Although the possibility of using the Blum skeleton in such domains has long been considered, only recently have robust and efficient algorithms for computing such representations been developed. This paper focuses on one class of such algorithms, those which utilize the average outward flux of the gradient vector field of the signed Euclidean distance function [7, 3, 6]. This paper makes several contributions. First, an alternate form of the divergence theorem is introduced in order to investigate the behaviour of the gradient vector field of the Euclidean distance function at medial points, where it is discontinuous. Second, we show that the iting behaviour of the average outward flux of this vector field through a Jordan curve as the area it encloses shrinks to zero, is very different at medial points than at non-medial ones, providing a strong theoretical justification for its use in the Hamilton- Jacobi skeletonization algorithm [7]. Third, we specialize to the case of shrinking circular neighborhoods and show that the average outward flux measure also reveals the object angle at each of the three types of skeletal points: 1) regular points, 2) branch points and 3) end- points. Hence, formulae for obtaining the boundary curves, their curvatures, and other geometric quantities of interest, can be written in terms of the average outward flux it values at skeletal points. We illustrate our results with several numerical simulations. 2. Flux and Divergence 2.1. reinaries We shall consider a shape to be any path connected space X in the plane, for which the boundary X is given by curves. For each point 2 in the plane we define a distance function D : 2 D( ) = inf d(, Q) Q X where d(, Q) is the Euclidean distance between and Q. Given our assumptions on X, D is a continuously differentiable function everywhere except on the skeleton. Now, for each point 2, define the set C = {(x, y) X D(x, y) = D( )}. Thus, C contains all points on the boundary that are closest to. The skeleton of a shape X, denoted S(X), is the set of points for which 1
2 there are two or more closest boundary points S(X) = { X C 2}. Technically, this excludes end-points of the skeleton since an end-point coincides with the center of curvature of a boundary point Q X that is a local maximum of positive curvature. The distance function gradient vector field q = D, is continuous and has magnitude one, except at skeletal points where it is not continuous. Hence, with the exception of skeletal points, D can be parametrized by its orientation at each point in the plane. With θ(x, y) : 2 [, 2π] being the orientation of the gradient vector at (x, y) and (θ) being the rotation matrix with argument θ, we have D(x, y) = (θ(x, y))e 1, where e 1 = (1, ), and D(x, y) v = cos( θ(x, y)), where v is any unit vector with orientation. We now define the outward flux and the average outward flux of the gradient vector field q through a Jordan curve. Let N be the outward normal at each point on the boundary of such a curve, and ds be an arc-length element. Definition 1 The outward flux of q through is given by q, N ds. Definition 2 The average outward flux of q through is given by q, N ds. ds An Extension of the Divergence Theorem If q is differentiable over a region with boundary, the standard form of the divergence theorem states that div( q) dv q, N ds. (1) Here dv is an area element. In other words, the integral of the divergence of the vector field over a region is given by the outward flux of the vector field through that region s bounding contour. We shall now develop an extension of the divergence theorem which can be applied to investigate properties of the vector field q = D at skeletal points, where it is discontinuous. Figure 1 illustrates the set up for the calculation which follows. Let S be a branch of the skeleton and let = 1 2 be a path connected region which intersects it. Let = C 1 C 2 and = S. Let t, C 3t be parallel curves to which approach as t. Let 1t and 2t be the regions obtained from 1 and 2 by removing the region between the curves C 3t and C 3t Finally, let q + denote q above S and q denote q below S. 2 1 C 1 C 2 C 3t C 3t Figure 1. A region which intersects a branch of the skeleton S. The outward flux of q through is given by q, N ds = q, N ds + q, N ds. C 1 C 2 Applying the divergence theorem to 1t and 2t div( q) dv = q, N ds + q, N ds, 1t C 1t C 3t div( q) dv = q, N ds + q, N ds. 2t C 2t S C 3t Adding the above two equations we have div( q) dv + div( q) dv = 1t 2t q, N ds + q, N ds + C 1t C 2t q, N ds + q, N ds. C 3t C 3t It is a standard property that the tangent to the skeleton bisects the the angle between q + and q at a skeletal point (see Figure 2). Thus, on we have q +, N + = q, N, (2) where N +, N denote the normals to from above and from below, respectively. Thus, one can take the it as t of both sides of the above equation to obtain the following extension of the divergence theorem Theorem 1 For a path connected region which contains part of a skeletal curve S, the divergence of the vector field q is related to its flux through by the following equation div( q) dv = = 1 2 q, N ds + 2 q, N C3 ds. 2
3 Although the divergence theorem fails for such regions because q is discontinuous on, the last integral is well defined due to Eq. 2; we either take q = q + and N C3 = N + or q = q and N C3 = N. Because q, N C3 is in the interval (, 1] on the skeleton, it also follows that q, N ds div( q) dv 2 length( ) Intrinsic Meaning of the Outward Flux It is a standard fact that the quantity div( q) dv, and hence the outward flux of q through, measures the degree to which the flow produced by the vector field q is area preserving in 2D (or volume preserving in 3D). This can be seen by applying the transport theorem [1] (p. 1) to the area enclosed by. This quantity is highly nonintrinsic in that it is strongly dependent on, and is zero only for very special cases, those where the vector field q is locally parallel. Thus, calculations such as those reported in [1], where the outward flux measure introduced in [6] is related to a notion of Blum s ligature, do not hold in the general case. This also raises the question of how such a measure could be used to distinguish skeletal points from non-skeletal ones, as suggested in [6, 3]. The answer, as we will now see, is to consider the iting behavior of the average outward flux of q = D, as the region enclosed shrinks to zero. This measure is at the heart of the Hamilton-Jacobi skeletonization algorithm in [7] Limiting Behaviour We now consider the it values of the outward flux and the average outward flux, of the vector field q = D through a convex curve that is the boundary of a region, as the region shrinks to a point. The results reported here are actually a special case of a more general result which applies to the case of any shrinking convex domain, and which extends to higher dimensions [2]. We begin by considering the case where the it point does not lie on the skeleton and hence Eq. (1) applies. We can write L div( q) dv q + δ(s), N ds where L is the Euclidean length of, q is the value of the vector field at the it point and δ(s) is the quantity added to get the value of q at neighboring points on ( q is continuous at ). In the it as L we have ) L ( L q, N ds + L δ(s), N ds It is easy to see that the second integral goes to zero due to the fact that δ(s) goes to zero. The first integral also goes to. zero due to the fact that is a closed curve. Thus, the it value obtained for the outward flux value at a non-skeletal point is zero. By essentially the same argument, the it value obtained for the average outward flux value at such a point is also zero. We now consider the second case where the it point is a skeletal point and hence Theorem 1 applies, which we rewrite as q, N ds = div( q) dv 2 q, N C3 ds. By the same argument presented for the first case above, the first integral on the right hand side goes to zero in the it as. The second integral on the right hand side can be written as 2 q + δ(s), N C3 ds where is the it point and δ(s) is the quantity added to get the value of q at neighboring points on. Owing to the fact that the integrand is in (, 1] (a special property of skeletal points), the value of this second integral is bounded above by 2 (inf C3 q C3, N C3 ) length( ) and below by 2 length( ). Thus we deduce that the it value of the outward flux at a skeletal point is also zero, due to the multiplicative length term. On the other hand, when the average outward flux is considered (see Definition 2), the it value reached as the region shrinks to a skeletal point is 2 q, N. Summarizing the above results, we have the property that whereas the it value of the outward flux is zero for both skeletal and non-skeletal points, the average outward flux has a different iting behaviours at skeletal points than at non-skeletal ones, providing a theoretical justification for its use in the Hamilton-Jacobi skeletonization algorithm [7]. 3. Circular Neighborhoods We now specialize the average outward flux calculation to the case of circular neighborhoods shrinking to a skeletal point. We shall treat the three cases of regular points, branch points and end-points of the skeleton separately egular Skeletal oints A regular skeletal point is one for which C = {Q 1, Q 2 } for Q 1 Q 2. Let n 1 and n 2 be the unit inward normals to the boundary at Q 1 and Q 2 respectively. Let t be the unit tangent vector to the skeleton at and define the object angle at to be ( ) [, π/2], such that n 1 n 2 = cos 2( ). It follows that n i t = cos ( ) for i = 1, 2 (see Figure 2). 3
4 S(t) Q 1 t Q 2 Figure 2. The object angle = ( ) at a simple skeletal point. Here S(t) is a parameterization of the skeleton curve. Hence, t = S (t ) is the tangent at t, i.e. where = S(t ). 1 C ε S(t) Figure 3. The distance function gradient vector field in the ε-neighborhood of is given by a step function one value for the top semi-circle and another for the bottom one. Both these vectors form an angle of = ( ) with t, since the skeleton is assumed to cut Cε in half at and 1. Now, let Cε be the circle with radius ε centered at. Let Cε : [, 2π] 2 be given by C ε (s) = ε(cos(s + θ(t ), sin(s + θ(t )) +, (3) where Cε () = + εt and Cε (π) = εt. Now consider Figure 3. Here, it is assumed that the gradient field has one value along Cε (s) for s (, π) and another for s (π, 2π). Also, both Cε () = and Cε (π) = 1 are on the skeleton. Let the outward normal of this circle at s be N (s). Hence, the outward flux of D though Cε (s) is F ε ( ) = 2πε = ε π = 4ε sin() D(C ε (s)), N (s) ds cos( s) ds ε 2π π cos( s) ds Notice that this calculation holds regardless of the orientation of t. However, it makes very strict assumptions that do not hold in most situations. Fortunately, the general case is similar to this one. There are only two differences: (1) C ε () and C ε (π) may not be on the skeleton, and (2) the distance function gradient field may take on more than two values along C ε (s) for s [, 2π]. For small enough ε, the circle will intersect the skeleton at precisely two points, which we label = C ε (δ ) and 1 = C ε (π + δ 1 ). Thus, the distance function gradient field is continuous on C ε (s) for s I = (δ, π+δ 1 ) and also for s I 1 = (π+δ 1, 2π δ ) 1. However, it may take on more than one value in the intervals I and I 1. Define β (s) and β 1 (s) on I and I 1 respectively, to account for such eventualities: t θ(c ε (s)) = cos (( ) + β (s)), s I t θ(c ε (s)) = cos ( ( ) + β 1 (s)), s I 1. Therefore, the outward flux calculation for regular skeletal points becomes F ε ( ) = 2πε π+δ1 = ε ε D(C ε (s)), N (s) ds δ 2π δ π+δ 1 cos( + β (s) s) ds cos( + β 1 (s) s) ds. The continuity of the distance function gradient field along the circle implies that both β (s) and β 1 (s) are continuous functions. Further, as ε, necessarily sup ε s [δ,π+δ 1] sup ε s [π+δ 1,2π δ ] β (s) = β 1 (s) =. Also, since the skeletal curve has continuous tangents, we must have that δ i = for i =, 1. Therefore the average outward flux through a shrinking circular region is ε given by F ε ( ) = 2 sin. ε 2πε π 3.2. Skeletal End-oints Let be a skeletal end-point. Let the point Q ε be on the branch which is at distance ε from. Choose ε small enough so that Q ε is a regular skeletal point. Thus, the object angle is well defined for Q ε. Now, let = ε (Q ε ). This it makes sense, because the circle 2 Cε intersects the skeleton at a single point. Also, the object angle varies continuously along a skeletal branch. 1 However, it is not necessarily continuous on the closure of I I 1. 2 Here C ε is as defined in Eq. (3) but t = ε t Qε. 4
5 S(t) C ε ε 2 Figure 4. A circular neighborhood of radius ε around the end-point. Along the arc of angle 2 the gradient vectors agree (in orientation) with the normals to C ε. Along the arc above S(t) the gradient vectors all form an angle of with S () = t. Similarly, for the arc below, this angle is. Now consider Figure 4. Along the arc arc opposite to the skeleton curve, the distance function gradient field must coincide with the inner normals of the circle. This is because the end-point results from the collapse of a circular arc (possibly a point if = ) on the boundary. On the rest of the circle, the distance function gradient field behaves as if were a regular skeletal point. Thus, F ε ( ) = ε ε π+δ ε ds cos( + β (s) s) ds 2π ε cos( + β 1 (s) s) ds π+δ where δ and β i (s) account for the circle not intersecting the skeleton midway and the distance function gradient field not being strictly a step function on Cε arc. Therefore, F ε ( ) = 1 ε 2πε π ( + sin ) since, as ε, δ, β (s) and β 1 (s) vanish. Notice, however, that = if the end-point is generated from a contour segment where the curvature is continuous Skeletal Junction oints Let be a skeletal junction point; that is where n skeletal curves meet. Let these curves be given by parameterizations S i (t) so that S i () =. Consider a circle of radius ε centered at. Denote it Cε. For small enough ε, Cε intersects the skeleton at precisely n regular points. efer to them as Q i ε = S i (ε). Hence, to each there is a corresponding object angle. Define i as i = ε Q i ε. S 1 (t) C ε S 2 (t) S3(t) Figure 5. A circular neighborhood of radius ε around the junction point. There are three skeletal curves denoted by S 1(t), S 2(t) and S 3(t) respectively. The dashed lines link and its closest points on the boundary (i.e. points in C). Note that = π. Now consider Figure 5. It suggests that i 2 i = 2π. Indeed, i is the angle between S i () 3 and the line joining to some point in C. To compute the outward flux through C ε, we can divide the circle into n arcs, each corresponding to a skeletal curve. In particular, for S i (t) this would be the arc of angle 2 i. For example, in Figure 5, the arc corresponding to S 1 (t) is the union of the two arcs of angle 1. Notice that the distance function gradient field along that arc behaves like that of a regular skeletal point with object angle i. Hence, the outward flux through it is i F ε (arc i ) = ε cos( i + β,i (s) s) ds δ i δi ε cos( i + β 1,i (s) s) ds i where δ i, β,i (s) and β,i (s) all vanish as ε. Thus, the total outward flux is F ε ( ) = n i=1 F ε(arc i ) and the average outward flux becomes F ε ( ) ε 2πε 3.4. Non-Skeletal oints = 1 π n sin i. i=1 Now, let be a non-skeletal point. In particular, there exists an ε small enough, so that Cε contains no skeletal points. Hence, the distance function gradient field along the circle is continuous. Thus, we can write F ε ( ) = ε 2π 3 Here S i () = t + S i (t). cos( + β(s) s) ds, 5
6 OINT TYE egular oints (Figure 2) F ε ( ) ε 2πε 2 π sin End-oints (Figure 4) 1 π (sin ) Junction oints (Figure 5) 1 π n i=1 sin i Non-Skeletal oints Table 1. A summary of results relating the it values of the average outward flux to the object angle for shrinking circular neighborhoods. Note that for contours of type, will be zero for the case of end-points. where ε sup s [,2π] β(s) =. Hence, F ε ( ) = 1 ε 2πε 2π 2π 3.5. Boundary econstruction cos( s) ds =. These results also suggest a method for parameterizing the boundary curves given a parameterization of the skeletal curves. Let S(t) be an arc-length parameterization of a skeletal curve and denote by (t) and r(t) the object angle at S(t) and distance to boundary, respectively (see Figure 2; if S(t) = then r(t) = d(, Q 1 ) = d(, Q 2 )). Note that the object angle can in turn be written in terms of the average outward flux it values, as summarized in Table 1. Assume that r(t) is monotonically decreasing, i.e., segment the skeletal curve if necessary, The two boundary curves that created this skeleton curve are given by C 1,2 (t) = S(t) + r(t)(±(t))s (t), (4) where S (t) is the tangent to the skeleton at S(t) and () denotes the rotation of S (t) by. 4. Experimental esults We have performed a several numerical simulations to corroborate findings as well as to demonstrate the applicability of our results. Figure 6 presents an overview of the steps involved in the experiments. We begin with a binary image and trace its discrete boundary to get an ordered list of pixels connected with line segments. Thus, a closed curve C(t) : 2 is obtained 4 (step 1). It is now possible to define a Euclidean distance function D : 2 on the whole plane. In this particular setting, D(x, y) can be computed exactly by considering all line segments and returning the distance to the segment closest to a point (x, y). Thus, it is also possible to approximate the average outward flux of D through a circle of radius ε > centered anywhere in 2. Furthermore, the precision of the computation depends only on the hardware. Hence, the results summarized in Table 1 can be used to obtain the skeleton of the curve. The exact method is described in [3]. This algorithm returns a discrete sampling of the skeletal curves as sequences of points (step 2). Thus, tangents at skeletal points may be approximated. The radius function is available through D(x, y), and the object angle can be obtained from the flux at regular points. This information can then be used to reconstruct boundary sections corresponding to regular skeletal points (step 3), see Section 3.5. Figure 7 provides additional examples of the above computation. Figure 8 compares the accuracy of the method to that of an exact calculation. The profile on the left uses straight lines to show the association of regular skeletal points with their bi-tangent points on the contour (the black circles). Here the association has been determined by using the average outward flux it values to obtain the object angle. The profile on the right demonstrates an exact computation, where the bi-tangent points are obtained by connecting each regular skeletal point to its two closest contour points. Notice how similar the two computations are. In fact, the average outward flux calculation can be performed with arbitrary precision 5. It should be pointed out that certain portions of the contour have not been reconstructed because the end-points of the skeleton have been chosen to satisfy an object angle threshold (above 3 ). Consequently, the end-points shown, although very close to the real ones 6, may not be actual end-points of the skeleton. Thus, to approximate the missing portions of the contour, it would be necessary to draw the circular arcs corresponding to the approximate end-points. 5. Conclusion We have derived an alternate form of the divergence theorem that accounts for singularities in the gradient vector field of the Euclidean distance function to the boundary of a shape and have used it to establish a criterion for detect- 4 The pixel locations are smoothed to account for jaggedness inherent in all discrete images. 5 Given a sub-sampling of a skeletal curve, it can be improved by numerical methods. One can use the fact that a skeletal point is a maximum along the gradient line leading to it. Thus, the Nelder-Mead Simplex method in 1D (see [4]) may be applied to precisely interpolate new points. 6 In fact it can be shown that over 9% of the area due to protrusions is kept, although a detailed analysis is beyond the scope of this paper. 6
7 1 2 3 Figure 6. (1) From a binary image, the boundary is extracted and represented as a continuous curve. (2) The skeleton is computed. (3) Using the average outward flux and radius values along the skeleton, the boundary is reconstructed. Figure 8. Bi-tangent points associated with regular skeletal points: (LEFT) computed using average outward flux information and (IGHT) computed explicitly. Figure 7. For each shape the original is shown in gray, the skeleton obtained using the average outward flux is shown with curves and the boundary points estimated from the skeleton using the relationship between the average outward flux and the object angle are shown with black circles. ing skeletal points. In particular, we have shown that the average outward flux through a shrinking region containing a skeletal point tends to a nonzero value, whereas it vanishes for non-skeletal points. An analysis assuming shrinking circular regions shows that the average outward flux also reveals the object angle for the three types of skeletal points, allowing for the boundary to be explicitly reconstructed given a parametrization of the skeleton. Thus the average outward flux measure can be viewed as a Euclidean invariant for shape description: it can be used to both detect the skeleton as well as to explicitly reconstruct the boundary from it. We have presented several numerical experiments to corroborate the theory. eferences [1] A. J. Chorin and J. E. Marden. A Mathematical Introduction to Fluid Mechanics. Springer-Verlag, [2] J. N. Damon. Global geometry of regions and boundaries via skeletal and medial integrals. manuscript in preparation, 23. [3]. Dimitrov, C. hillips, and K. Siddiqi. obust and efficient skeletal graphs. In CV, volume 1, pages , 2. [4] J. C. Lagarias, J. A. eeds, M. H. Wright, and. E. Wright. Convergence properties of the nelder-mead simplex algorithm in low dimensions. SIAM Journal on Optimization, 9: , [5] T. Sebastian,. Kline, and B. B. Kimia. ecognition of shapes by editing shock graphs. In International Conference on Computer Vision, pages , Vancouver, Canada, July 21. [6] K. Siddiqi, S. Bouix, A. Tannenbaum, and S. W. Zucker. The hamilton-jacobi skeleton. In International Conference on Computer Vision, pages , Kerkyra, Greece, September [7] K. Siddiqi, S. Bouix, A. Tannenbaum, and S. W. Zucker. Hamilton-jacobi skeletons. International Journal of Computer Vision, 48(3): , 22. [8] K. Siddiqi, A. Shokoufandeh, S. J. Dickinson, and S. W. Zucker. Shock graphs and shape matching. International Journal of Computer Vision, 35(1):13 32, [9] G. D. Stetten and S. M. izer. Automated identification and measurement of objects via populations of medial primitives, with application to real time 3d echocardiography. In International Conference on Information rocessing in Medical Imaging, pages 84 97, [1] A. Torsello and E.. Hancock. A skeletal measure of 2d shape similarity. In roceedings of the International Workshop On Visual Form, volume LNCS 259, pages , 21. 7
On a Differential Equation Arising in Plant Vascular Biology
TR roof of Result by. Dimitrov and S. W. Zucker draft February 10, 2006 p. 1 On a Differential Equation Arising in lant Vascular Biology avel Dimitrov and Steven W. Zucker February 10, 2006 Abstract The
More informationThe Hamilton-Jacobi Skeleton
The Hamilton-Jacobi Skeleton Kaleem Siddiqi Λ Sylvain Bouix y Allen Tannenbaum z Steven W. Zucker x Abstract The eikonal equation and variants of it are of significant interest for problems in computer
More informationMath 263 Final. (b) The cross product is. i j k c. =< c 1, 1, 1 >
Math 63 Final Problem 1: [ points, 5 points to each part] Given the points P : (1, 1, 1), Q : (1,, ), R : (,, c 1), where c is a parameter, find (a) the vector equation of the line through P and Q. (b)
More informationENGI 4430 Line Integrals; Green s Theorem Page 8.01
ENGI 443 Line Integrals; Green s Theorem Page 8. 8. Line Integrals Two applications of line integrals are treated here: the evaluation of work done on a particle as it travels along a curve in the presence
More informationContents. MATH 32B-2 (18W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables. 1 Multiple Integrals 3. 2 Vector Fields 9
MATH 32B-2 (8W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables Contents Multiple Integrals 3 2 Vector Fields 9 3 Line and Surface Integrals 5 4 The Classical Integral Theorems 9 MATH 32B-2 (8W)
More informationThe Minimum Speed for a Blocking Problem on the Half Plane
The Minimum Speed for a Blocking Problem on the Half Plane Alberto Bressan and Tao Wang Department of Mathematics, Penn State University University Park, Pa 16802, USA e-mails: bressan@mathpsuedu, wang
More informationMATH 332: Vector Analysis Summer 2005 Homework
MATH 332, (Vector Analysis), Summer 2005: Homework 1 Instructor: Ivan Avramidi MATH 332: Vector Analysis Summer 2005 Homework Set 1. (Scalar Product, Equation of a Plane, Vector Product) Sections: 1.9,
More informationMATHEMATICS 317 December 2010 Final Exam Solutions
MATHEMATI 317 December 1 Final Eam olutions 1. Let r(t) = ( 3 cos t, 3 sin t, 4t ) be the position vector of a particle as a function of time t. (a) Find the velocity of the particle as a function of time
More informationGLOBAL GEOMETRY OF REGIONS AND BOUNDARIES VIA SKELETAL AND MEDIAL INTEGRALS
GLOBAL GEOMETRY OF REGIONS AND BOUNDARIES VIA SKELETAL AND MEDIAL INTEGRALS JAMES DAMON Department of Mathematics University of North Carolina at Chapel Hill, Chapel Hill, NC 27599-3250, USA PRELIMINARY
More informationHOMEWORK 2 SOLUTIONS
HOMEWORK SOLUTIONS MA11: ADVANCED CALCULUS, HILARY 17 (1) Find parametric equations for the tangent line of the graph of r(t) = (t, t + 1, /t) when t = 1. Solution: A point on this line is r(1) = (1,,
More informationVector Calculus, Maths II
Section A Vector Calculus, Maths II REVISION (VECTORS) 1. Position vector of a point P(x, y, z) is given as + y and its magnitude by 2. The scalar components of a vector are its direction ratios, and represent
More informationLECTURE-15 : LOGARITHMS AND COMPLEX POWERS
LECTURE-5 : LOGARITHMS AND COMPLEX POWERS VED V. DATAR The purpose of this lecture is twofold - first, to characterize domains on which a holomorphic logarithm can be defined, and second, to show that
More informationGauss s Law & Potential
Gauss s Law & Potential Lecture 7: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay Flux of an Electric Field : In this lecture we introduce Gauss s law which happens to
More informationPlanar interpolation with a pair of rational spirals T. N. T. Goodman 1 and D. S. Meek 2
Planar interpolation with a pair of rational spirals T N T Goodman and D S Meek Abstract Spirals are curves of one-signed monotone increasing or decreasing curvature Spiral segments are fair curves with
More informationMATH 2433 Homework 1
MATH 433 Homework 1 1. The sequence (a i ) is defined recursively by a 1 = 4 a i+1 = 3a i find a closed formula for a i in terms of i.. In class we showed that the Fibonacci sequence (a i ) defined by
More informationThe Divergence Theorem Stokes Theorem Applications of Vector Calculus. Calculus. Vector Calculus (III)
Calculus Vector Calculus (III) Outline 1 The Divergence Theorem 2 Stokes Theorem 3 Applications of Vector Calculus The Divergence Theorem (I) Recall that at the end of section 12.5, we had rewritten Green
More informationVectors. January 13, 2013
Vectors January 13, 2013 The simplest tensors are scalars, which are the measurable quantities of a theory, left invariant by symmetry transformations. By far the most common non-scalars are the vectors,
More informationDefinite integrals. We shall study line integrals of f (z). In order to do this we shall need some preliminary definitions.
5. OMPLEX INTEGRATION (A) Definite integrals Integrals are extremely important in the study of functions of a complex variable. The theory is elegant, and the proofs generally simple. The theory is put
More informationChapter 31. The Laplace Transform The Laplace Transform. The Laplace transform of the function f(t) is defined. e st f(t) dt, L[f(t)] =
Chapter 3 The Laplace Transform 3. The Laplace Transform The Laplace transform of the function f(t) is defined L[f(t)] = e st f(t) dt, for all values of s for which the integral exists. The Laplace transform
More informationMATH 0350 PRACTICE FINAL FALL 2017 SAMUEL S. WATSON. a c. b c.
MATH 35 PRACTICE FINAL FALL 17 SAMUEL S. WATSON Problem 1 Verify that if a and b are nonzero vectors, the vector c = a b + b a bisects the angle between a and b. The cosine of the angle between a and c
More informationII. DIFFERENTIABLE MANIFOLDS. Washington Mio CENTER FOR APPLIED VISION AND IMAGING SCIENCES
II. DIFFERENTIABLE MANIFOLDS Washington Mio Anuj Srivastava and Xiuwen Liu (Illustrations by D. Badlyans) CENTER FOR APPLIED VISION AND IMAGING SCIENCES Florida State University WHY MANIFOLDS? Non-linearity
More informationSOLUTIONS TO THE FINAL EXAM. December 14, 2010, 9:00am-12:00 (3 hours)
SOLUTIONS TO THE 18.02 FINAL EXAM BJORN POONEN December 14, 2010, 9:00am-12:00 (3 hours) 1) For each of (a)-(e) below: If the statement is true, write TRUE. If the statement is false, write FALSE. (Please
More informationPractice Problems for Exam 3 (Solutions) 1. Let F(x, y) = xyi+(y 3x)j, and let C be the curve r(t) = ti+(3t t 2 )j for 0 t 2. Compute F dr.
1. Let F(x, y) xyi+(y 3x)j, and let be the curve r(t) ti+(3t t 2 )j for t 2. ompute F dr. Solution. F dr b a 2 2 F(r(t)) r (t) dt t(3t t 2 ), 3t t 2 3t 1, 3 2t dt t 3 dt 1 2 4 t4 4. 2. Evaluate the line
More informationFraunhofer Institute for Computer Graphics Research Interactive Graphics Systems Group, TU Darmstadt Fraunhoferstrasse 5, Darmstadt, Germany
Scale Space and PDE methods in image analysis and processing Arjan Kuijper Fraunhofer Institute for Computer Graphics Research Interactive Graphics Systems Group, TU Darmstadt Fraunhoferstrasse 5, 64283
More informationMA 323 Geometric Modelling Course Notes: Day 07 Parabolic Arcs
MA 323 Geometric Modelling Course Notes: Day 07 Parabolic Arcs David L. Finn December 9th, 2004 We now start considering the basic curve elements to be used throughout this course; polynomial curves and
More informationCalculus III. Math 233 Spring Final exam May 3rd. Suggested solutions
alculus III Math 33 pring 7 Final exam May 3rd. uggested solutions This exam contains twenty problems numbered 1 through. All problems are multiple choice problems, and each counts 5% of your total score.
More information8. THE FARY-MILNOR THEOREM
Math 501 - Differential Geometry Herman Gluck Tuesday April 17, 2012 8. THE FARY-MILNOR THEOREM The curvature of a smooth curve in 3-space is 0 by definition, and its integral w.r.t. arc length, (s) ds,
More informationAppendix E : Note on regular curves in Euclidean spaces
Appendix E : Note on regular curves in Euclidean spaces In Section III.5 of the course notes we posed the following question: Suppose that U is a connected open subset of R n and x, y U. Is there a continuous
More informationSolutions for the Practice Final - Math 23B, 2016
olutions for the Practice Final - Math B, 6 a. True. The area of a surface is given by the expression d, and since we have a parametrization φ x, y x, y, f x, y with φ, this expands as d T x T y da xy
More informationPRACTICE PROBLEMS. Please let me know if you find any mistakes in the text so that i can fix them. 1. Mixed partial derivatives.
PRACTICE PROBLEMS Please let me know if you find any mistakes in the text so that i can fix them. 1.1. Let Show that f is C 1 and yet How is that possible? 1. Mixed partial derivatives f(x, y) = {xy x
More information18.02 Multivariable Calculus Fall 2007
MIT OpenourseWare http://ocw.mit.edu 18.02 Multivariable alculus Fall 2007 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 18.02 Lecture 21. Test for
More informationSet, functions and Euclidean space. Seungjin Han
Set, functions and Euclidean space Seungjin Han September, 2018 1 Some Basics LOGIC A is necessary for B : If B holds, then A holds. B A A B is the contraposition of B A. A is sufficient for B: If A holds,
More informationPhysics 2400 Midterm I Sample March 2017
Physics 4 Midterm I Sample March 17 Question: 1 3 4 5 Total Points: 1 1 1 1 6 Gamma function. Leibniz s rule. 1. (1 points) Find positive x that minimizes the value of the following integral I(x) = x+1
More informationSummary of various integrals
ummary of various integrals Here s an arbitrary compilation of information about integrals Moisés made on a cold ecember night. 1 General things o not mix scalars and vectors! In particular ome integrals
More informationLECTURE 15: COMPLETENESS AND CONVEXITY
LECTURE 15: COMPLETENESS AND CONVEXITY 1. The Hopf-Rinow Theorem Recall that a Riemannian manifold (M, g) is called geodesically complete if the maximal defining interval of any geodesic is R. On the other
More informationName: Instructor: Lecture time: TA: Section time:
Math 222 Final May 11, 29 Name: Instructor: Lecture time: TA: Section time: INSTRUCTIONS READ THIS NOW This test has 1 problems on 16 pages worth a total of 2 points. Look over your test package right
More informationMath 233. Practice Problems Chapter 15. i j k
Math 233. Practice Problems hapter 15 1. ompute the curl and divergence of the vector field F given by F (4 cos(x 2 ) 2y)i + (4 sin(y 2 ) + 6x)j + (6x 2 y 6x + 4e 3z )k olution: The curl of F is computed
More informationGeometric approximation of curves and singularities of secant maps Ghosh, Sunayana
University of Groningen Geometric approximation of curves and singularities of secant maps Ghosh, Sunayana IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish
More informationTHEODORE VORONOV DIFFERENTIABLE MANIFOLDS. Fall Last updated: November 26, (Under construction.)
4 Vector fields Last updated: November 26, 2009. (Under construction.) 4.1 Tangent vectors as derivations After we have introduced topological notions, we can come back to analysis on manifolds. Let M
More informationMAC2313 Final A. (5 pts) 1. How many of the following are necessarily true? i. The vector field F = 2x + 3y, 3x 5y is conservative.
MAC2313 Final A (5 pts) 1. How many of the following are necessarily true? i. The vector field F = 2x + 3y, 3x 5y is conservative. ii. The vector field F = 5(x 2 + y 2 ) 3/2 x, y is radial. iii. All constant
More informationGLOBAL GEOMETRY OF REGIONS AND BOUNDARIES VIA SKELETAL AND MEDIAL INTEGRALS
GLOBAL GEOMETRY OF REGIONS AND BOUNDARIES VIA SKELETAL AND MEDIAL INTEGRALS JAMES DAMON Abstract. For a compact region Ω in R n+ with smooth generic boundary B, the Blum medial axis M is the locus of centers
More informationMathematics (Course B) Lent Term 2005 Examples Sheet 2
N12d Natural Sciences, Part IA Dr M. G. Worster Mathematics (Course B) Lent Term 2005 Examples Sheet 2 Please communicate any errors in this sheet to Dr Worster at M.G.Worster@damtp.cam.ac.uk. Note that
More informationCorrections to the First Printing: Chapter 6. This list includes corrections and clarifications through November 6, 1999.
June 2,2 1 Corrections to the First Printing: Chapter 6 This list includes corrections and clarifications through November 6, 1999. We should have mentioned that our treatment of differential forms, especially
More information1 Differentiable manifolds and smooth maps
1 Differentiable manifolds and smooth maps Last updated: April 14, 2011. 1.1 Examples and definitions Roughly, manifolds are sets where one can introduce coordinates. An n-dimensional manifold is a set
More informationThe Extended Evolute as the Transition Set of Distance Functions
The Extended Evolute as the Transition Set of Distance Functions Xianming Chen Richard F. Riesenfeld Elaine Cohen Abstract This paper addresses three practical issues arising from applying singularity
More informationImage Analysis. Feature extraction: corners and blobs
Image Analysis Feature extraction: corners and blobs Christophoros Nikou cnikou@cs.uoi.gr Images taken from: Computer Vision course by Svetlana Lazebnik, University of North Carolina at Chapel Hill (http://www.cs.unc.edu/~lazebnik/spring10/).
More informationPrint Your Name: Your Section:
Print Your Name: Your Section: Mathematics 1c. Practice Final Solutions This exam has ten questions. J. Marsden You may take four hours; there is no credit for overtime work No aids (including notes, books,
More information53. Flux Integrals. Here, R is the region over which the double integral is evaluated.
53. Flux Integrals Let be an orientable surface within 3. An orientable surface, roughly speaking, is one with two distinct sides. At any point on an orientable surface, there exists two normal vectors,
More information4B. Line Integrals in the Plane
4. Line Integrals in the Plane 4A. Plane Vector Fields 4A-1 Describe geometrically how the vector fields determined by each of the following vector functions looks. Tell for each what the largest region
More informationMultiple Integrals and Vector Calculus (Oxford Physics) Synopsis and Problem Sets; Hilary 2015
Multiple Integrals and Vector Calculus (Oxford Physics) Ramin Golestanian Synopsis and Problem Sets; Hilary 215 The outline of the material, which will be covered in 14 lectures, is as follows: 1. Introduction
More informationCURRENT MATERIAL: Vector Calculus.
Math 275, section 002 (Ultman) Spring 2012 FINAL EXAM REVIEW The final exam will be held on Wednesday 9 May from 8:00 10:00am in our regular classroom. You will be allowed both sides of two 8.5 11 sheets
More informationG G. G. x = u cos v, y = f(u), z = u sin v. H. x = u + v, y = v, z = u v. 1 + g 2 x + g 2 y du dv
1. Matching. Fill in the appropriate letter. 1. ds for a surface z = g(x, y) A. r u r v du dv 2. ds for a surface r(u, v) B. r u r v du dv 3. ds for any surface C. G x G z, G y G z, 1 4. Unit normal N
More informationISOPERIMETRIC INEQUALITY FOR FLAT SURFACES
Proceedings of The Thirteenth International Workshop on Diff. Geom. 3(9) 3-9 ISOPERIMETRIC INEQUALITY FOR FLAT SURFACES JAIGYOUNG CHOE Korea Institute for Advanced Study, Seoul, 3-7, Korea e-mail : choe@kias.re.kr
More informationTHE FOUR VERTEX THEOREM AND ITS CONVERSE in honor of Björn Dahlberg
Writhe - Four Vertex Theorem.14 June, 2006 THE FOUR VERTEX THEOREM AND ITS CONVERSE in honor of Björn Dahlberg The Four Vertex Theorem says that a simple closed curve in the plane, other than a circle,
More informationChaotic Billiards. Part I Introduction to Dynamical Billiards. 1 Review of the Dierent Billiard Systems. Ben Parker and Alex Riina.
Chaotic Billiards Ben Parker and Alex Riina December 3, 2009 Part I Introduction to Dynamical Billiards 1 Review of the Dierent Billiard Systems In investigating dynamical billiard theory, we focus on
More informationMath 212-Lecture 20. P dx + Qdy = (Q x P y )da. C
15. Green s theorem Math 212-Lecture 2 A simple closed curve in plane is one curve, r(t) : t [a, b] such that r(a) = r(b), and there are no other intersections. The positive orientation is counterclockwise.
More informationENGI 4430 Line Integrals; Green s Theorem Page 8.01
ENGI 4430 Line Integrals; Green s Theorem Page 8.01 8. Line Integrals Two applications of line integrals are treated here: the evaluation of work done on a particle as it travels along a curve in the presence
More informationexample consider flow of water in a pipe. At each point in the pipe, the water molecule has a velocity
Module 1: A Crash Course in Vectors Lecture 1: Scalar and Vector Fields Objectives In this lecture you will learn the following Learn about the concept of field Know the difference between a scalar field
More informationMinimal Surface equations non-solvability strongly convex functional further regularity Consider minimal surface equation.
Lecture 7 Minimal Surface equations non-solvability strongly convex functional further regularity Consider minimal surface equation div + u = ϕ on ) = 0 in The solution is a critical point or the minimizer
More informationArc Length. Philippe B. Laval. Today KSU. Philippe B. Laval (KSU) Arc Length Today 1 / 12
Philippe B. Laval KSU Today Philippe B. Laval (KSU) Arc Length Today 1 / 12 Introduction In this section, we discuss the notion of curve in greater detail and introduce the very important notion of arc
More informationMaxwell s equations for electrostatics
Maxwell s equations for electrostatics October 6, 5 The differential form of Gauss s law Starting from the integral form of Gauss s law, we treat the charge as a continuous distribution, ρ x. Then, letting
More informationMath 31CH - Spring Final Exam
Math 3H - Spring 24 - Final Exam Problem. The parabolic cylinder y = x 2 (aligned along the z-axis) is cut by the planes y =, z = and z = y. Find the volume of the solid thus obtained. Solution:We calculate
More information49. Green s Theorem. The following table will help you plan your calculation accordingly. C is a simple closed loop 0 Use Green s Theorem
49. Green s Theorem Let F(x, y) = M(x, y), N(x, y) be a vector field in, and suppose is a path that starts and ends at the same point such that it does not cross itself. Such a path is called a simple
More informationCreated by T. Madas LINE INTEGRALS. Created by T. Madas
LINE INTEGRALS LINE INTEGRALS IN 2 DIMENSIONAL CARTESIAN COORDINATES Question 1 Evaluate the integral ( x + 2y) dx, C where C is the path along the curve with equation y 2 = x + 1, from ( ) 0,1 to ( )
More informationVector Fields and Line Integrals The Fundamental Theorem for Line Integrals
Math 280 Calculus III Chapter 16 Sections: 16.1, 16.2 16.3 16.4 16.5 Topics: Vector Fields and Line Integrals The Fundamental Theorem for Line Integrals Green s Theorem Curl and Divergence Section 16.1
More informationarxiv: v3 [math.dg] 19 Jun 2017
THE GEOMETRY OF THE WIGNER CAUSTIC AND AFFINE EQUIDISTANTS OF PLANAR CURVES arxiv:160505361v3 [mathdg] 19 Jun 017 WOJCIECH DOMITRZ, MICHA L ZWIERZYŃSKI Abstract In this paper we study global properties
More informationAn Investigation of the Four Vertex Theorem and its Converse
Union College Union Digital Works Honors Theses Student Work 6-2017 An Investigation of the Four Vertex Theorem and its Converse Rebeka Kelmar Union College - Schenectady, NY Follow this and additional
More informationModule 16 : Line Integrals, Conservative fields Green's Theorem and applications. Lecture 48 : Green's Theorem [Section 48.
Module 16 : Line Integrals, Conservative fields Green's Theorem and applications Lecture 48 : Green's Theorem [Section 48.1] Objectives In this section you will learn the following : Green's theorem which
More informationComplex Analysis MATH 6300 Fall 2013 Homework 4
Complex Analysis MATH 6300 Fall 2013 Homework 4 Due Wednesday, December 11 at 5 PM Note that to get full credit on any problem in this class, you must solve the problems in an efficient and elegant manner,
More informationMorphological image processing
INF 4300 Digital Image Analysis Morphological image processing Fritz Albregtsen 09.11.2017 1 Today Gonzalez and Woods, Chapter 9 Except sections 9.5.7 (skeletons), 9.5.8 (pruning), 9.5.9 (reconstruction)
More informationConductors and Insulators
Conductors and Insulators Lecture 11: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay Self Energy of a Charge Distribution : In Lecture 1 we briefly discussed what we called
More informationApproximation of Circular Arcs by Parametric Polynomials
Approximation of Circular Arcs by Parametric Polynomials Emil Žagar Lecture on Geometric Modelling at Charles University in Prague December 6th 2017 1 / 44 Outline Introduction Standard Reprezentations
More informationMath 20C Homework 2 Partial Solutions
Math 2C Homework 2 Partial Solutions Problem 1 (12.4.14). Calculate (j k) (j + k). Solution. The basic properties of the cross product are found in Theorem 2 of Section 12.4. From these properties, we
More informationLecture notes: Introduction to Partial Differential Equations
Lecture notes: Introduction to Partial Differential Equations Sergei V. Shabanov Department of Mathematics, University of Florida, Gainesville, FL 32611 USA CHAPTER 1 Classification of Partial Differential
More informationCorner. Corners are the intersections of two edges of sufficiently different orientations.
2D Image Features Two dimensional image features are interesting local structures. They include junctions of different types like Y, T, X, and L. Much of the work on 2D features focuses on junction L,
More informationMATH H53 : Final exam
MATH H53 : Final exam 11 May, 18 Name: You have 18 minutes to answer the questions. Use of calculators or any electronic items is not permitted. Answer the questions in the space provided. If you run out
More informationIn English, this means that if we travel on a straight line between any two points in C, then we never leave C.
Convex sets In this section, we will be introduced to some of the mathematical fundamentals of convex sets. In order to motivate some of the definitions, we will look at the closest point problem from
More information1 The Differential Geometry of Surfaces
1 The Differential Geometry of Surfaces Three-dimensional objects are bounded by surfaces. This section reviews some of the basic definitions and concepts relating to the geometry of smooth surfaces. 1.1
More informationẋ = f(x, y), ẏ = g(x, y), (x, y) D, can only have periodic solutions if (f,g) changes sign in D or if (f,g)=0in D.
4 Periodic Solutions We have shown that in the case of an autonomous equation the periodic solutions correspond with closed orbits in phase-space. Autonomous two-dimensional systems with phase-space R
More informationPractice Final Solutions
Practice Final Solutions Math 1, Fall 17 Problem 1. Find a parameterization for the given curve, including bounds on the parameter t. Part a) The ellipse in R whose major axis has endpoints, ) and 6, )
More informationLength of parallel curves and rotation index
Length of parallel curves and rotation index E. Macías-Virgós 1 Institute of Mathematics. Department of Geometry and Topology. University of Santiago de Compostela. 15782- SPAIN Abstract We prove that
More informationClassical transcendental curves
Classical transcendental curves Reinhard Schultz May, 2008 In his writings on coordinate geometry, Descartes emphasized that he was only willing to work with curves that could be defined by algebraic equations.
More informationSolutions to old Exam 3 problems
Solutions to old Exam 3 problems Hi students! I am putting this version of my review for the Final exam review here on the web site, place and time to be announced. Enjoy!! Best, Bill Meeks PS. There are
More informationDIFFERENTIAL GEOMETRY, LECTURE 16-17, JULY 14-17
DIFFERENTIAL GEOMETRY, LECTURE 16-17, JULY 14-17 6. Geodesics A parametrized line γ : [a, b] R n in R n is straight (and the parametrization is uniform) if the vector γ (t) does not depend on t. Thus,
More information3 Chapter. Gauss s Law
3 Chapter Gauss s Law 3.1 Electric Flux... 3-2 3.2 Gauss s Law (see also Gauss s Law Simulation in Section 3.10)... 3-4 Example 3.1: Infinitely Long Rod of Uniform Charge Density... 3-9 Example 3.2: Infinite
More informationManifold Regularization
9.520: Statistical Learning Theory and Applications arch 3rd, 200 anifold Regularization Lecturer: Lorenzo Rosasco Scribe: Hooyoung Chung Introduction In this lecture we introduce a class of learning algorithms,
More informationfor all subintervals I J. If the same is true for the dyadic subintervals I D J only, we will write ϕ BMO d (J). In fact, the following is true
3 ohn Nirenberg inequality, Part I A function ϕ L () belongs to the space BMO() if sup ϕ(s) ϕ I I I < for all subintervals I If the same is true for the dyadic subintervals I D only, we will write ϕ BMO
More informationA DIFFERENTIAL GEOMETRIC APPROACH TO FLUID MECHANICS
International Journal of Scientific and Research Publications Volume 5 Issue 9 September 2015 1 A DIFFERENTIAL GEOMETRIC APPROACH TO FLUID MECHANICS Mansour Hassan Mansour * M A Bashir ** * Department
More information(x 1, y 1 ) = (x 2, y 2 ) if and only if x 1 = x 2 and y 1 = y 2.
1. Complex numbers A complex number z is defined as an ordered pair z = (x, y), where x and y are a pair of real numbers. In usual notation, we write z = x + iy, where i is a symbol. The operations of
More informationHIGHER ORDER RIGIDITY - WHAT IS THE PROPER DEFINITION?
HIGHER ORDER RIGIDITY - WHAT IS THE PROPER DEFINITION? ROBERT CONNELLY AND HERMAN SERVATIUS Abstract. We show that there is a bar and joint framework G(p) which has a configuration p in the plane such
More informationarxiv: v1 [math.dg] 11 Nov 2007
Length of parallel curves arxiv:711.167v1 [math.dg] 11 Nov 27 E. Macías-Virgós Abstract We prove that the length difference between a closed periodic curve and its parallel curve at a sufficiently small
More informationLECTURE-13 : GENERALIZED CAUCHY S THEOREM
LECTURE-3 : GENERALIZED CAUCHY S THEOREM VED V. DATAR The aim of this lecture to prove a general form of Cauchy s theorem applicable to multiply connected domains. We end with computations of some real
More informationThe fundamental theorem of calculus for definite integration helped us to compute If has an anti-derivative,
Module 16 : Line Integrals, Conservative fields Green's Theorem and applications Lecture 47 : Fundamental Theorems of Calculus for Line integrals [Section 47.1] Objectives In this section you will learn
More informationMath 5378, Differential Geometry Solutions to practice questions for Test 2
Math 5378, Differential Geometry Solutions to practice questions for Test 2. Find all possible trajectories of the vector field w(x, y) = ( y, x) on 2. Solution: A trajectory would be a curve (x(t), y(t))
More informationMcGill University April 16, Advanced Calculus for Engineers
McGill University April 16, 2014 Faculty of cience Final examination Advanced Calculus for Engineers Math 264 April 16, 2014 Time: 6PM-9PM Examiner: Prof. R. Choksi Associate Examiner: Prof. A. Hundemer
More informationMORE CONSEQUENCES OF CAUCHY S THEOREM
MOE CONSEQUENCES OF CAUCHY S THEOEM Contents. The Mean Value Property and the Maximum-Modulus Principle 2. Morera s Theorem and some applications 3 3. The Schwarz eflection Principle 6 We have stated Cauchy
More informationSTOKES THEOREM ON MANIFOLDS
STOKES THEOREM ON MANIFOLDS GIDEON DRESDNER Abstract. The generalization of the Fundamental Theorem of Calculus to higher dimensions requires fairly sophisticated geometric and algebraic machinery. In
More informationWORKSHEET #13 MATH 1260 FALL 2014
WORKSHEET #3 MATH 26 FALL 24 NOT DUE. Short answer: (a) Find the equation of the tangent plane to z = x 2 + y 2 at the point,, 2. z x (, ) = 2x = 2, z y (, ) = 2y = 2. So then the tangent plane equation
More informationPhys 122 Lecture 3 G. Rybka
Phys 122 Lecture 3 G. Rybka A few more Demos Electric Field Lines Example Calculations: Discrete: Electric Dipole Overview Continuous: Infinite Line of Charge Next week Labs and Tutorials begin Electric
More informationName: SOLUTIONS Date: 11/9/2017. M20550 Calculus III Tutorial Worksheet 8
Name: SOLUTIONS Date: /9/7 M55 alculus III Tutorial Worksheet 8. ompute R da where R is the region bounded by x + xy + y 8 using the change of variables given by x u + v and y v. Solution: We know R is
More information