396(1).wxm 1 / 14. (%i3) dr_nabla1(psi) := dx*diff(psi,x) + dy*diff(psi,y) + dz*diff(psi,z); d. d Y +dz d
|
|
- Sophie Terry
- 6 years ago
- Views:
Transcription
1 396(1).wxm 1 / 14 (%i1) kill(all); (%o0) done 1 Definitions (%i1) cross(a,b) := [a[]*b[3] - a[3]*b[], a[3]*b[1] - a[1]*b[3], a[1]*b[] - a[]*b[1]]; (%o1) cross( a,b ):=[a b 3 a 3 b,a 3 b 1 a 1 b 3,a 1 b a b 1 ] (%i) dr: [,, d]; (%o) [,,d] 1.1 Components of Eq.(14) (%i3) dr_nabla1(psi) := *diff(psi,x) + *diff(psi,y) + d*diff(psi,); d (%o3) dr_nabla1( ):= d X + d d Y +d d d (%i4) dr_nabla(psi) := (block [fh], fh: dr_nabla1(psi), *diff(fh,x) + *diff(fh,y) + d*diff(fh,) ); (%o4) dr_nabla( ):= d block fh,fh:dr_nabla1( ), d X fh + d d Y fh +d d d fh (%i5) dr_nabla3(psi) := (block [fh], fh: dr_nabla(psi), *diff(fh,x) + *diff(fh,y) + d*diff(fh,) ); (%o5) dr_nabla3( ):= d block fh,fh:dr_nabla( ), d X fh + d d Y fh +d d d fh (%i6) dr_nabla4(psi) := (block [fh], fh: dr_nabla3(psi), *diff(fh,x) + *diff(fh,y) + d*diff(fh,) ); (%o6) dr_nabla4( ):= d block fh,fh:dr_nabla3( ), d X fh + d d Y fh +d d d fh (%i7) dr_nabla5(psi) := (block [fh], fh: dr_nabla4(psi), *diff(fh,x) + *diff(fh,y) + d*diff(fh,) ); (%o7) dr_nabla5( ):= d block fh,fh:dr_nabla4( ), d X fh + d d Y fh +d d d fh
2 396(1).wxm / 14 (%i8) dr_nabla6(psi) := (block [fh], fh: dr_nabla5(psi), *diff(fh,x) + *diff(fh,y) + d*diff(fh,) ); (%o8) dr_nabla6( ):= d block fh,fh:dr_nabla5( ), d X fh + d d Y fh +d d d fh Averaging functions (%i9) /* Function for averaging delta r terms in an expression */ averagem(f, dr) := block([arglist, s0, den, i, j, res, l, el], /* make list of all summands in the formula */ arglist: [], s0: string(0), arglist: args(expand(f)), /* check if dr[i] appears in any denominator */ for j:1 thru length(arglist) do ( den: denom(arglist[j]), /*print ("DEN:", den),*/ for i:1 thru length(dr) do ( if not freeof(dr[i],den) then print(i, "Denominator error!") )), /* remove all list elements containing occurences of dr[i], dr[i]^3, for j:1 thru length(arglist) do ( for i:1 thru length(dr) do ( if not equal(string(coeff(arglist[j], dr[i], 1)), s0) or not equal(string(coeff(arglist[j], dr[i], 3)), s0) or not equal(string(coeff(arglist[j], dr[i], 5)), s0) then arglist[j]: 0 /*print("l(", dr[i], ") ", arglist[j])*/ )), res: 0, /* construct result expression */ l: length(arglist), for j:1 thru l do ( el: pop(arglist), res: res + el /*print(j,el,res)*/ ), res /*print(res)*/ )$
3 396(1).wxm 3 / 14 (%i10) /* Function for replacing deltax^ etc. to delta r^ terms in cartesian coordinates */ replacer(f,dri) := block([f1, i, dr, len], f1: f, len: length(dri), for i:1 thru len do ( f1: ratsubst(1/len*dr[av]^, dri[i]^, f1) ), factorsum(f1) )$ 3 General formulas for Taylor expansion with aver 3.1 1st order (%i11) depends(f, [X,Y,]); (%o11) [f( X,Y, )] (%i1) dr_nabla1(f); d (%o1) d d f + d d Y f + d d X f (%i14) f1: ratsimp(averagem(dr_nabla1(f), dr)); (%o13) 0 (%o14) 0 3. nd order (%i15) depends(f, [X,Y,]); (%o15) [f( X,Y, )] (%i16) Df: 1/!*dr_nabla(f); (%o16) (d d d d f + d d f +d d Y d d X f +d d d Y d f + d Y d f + d d X d Y f d d X d f + d d X d f + d d X d Y f )/ +
4 396(1).wxm 4 / 14 (%i18) f1: ratsimp(averagem(df, dr)); d d d f (%o17) + d d Y f + d d X f (%o18) dr av d d f rd order d d Y f + (%i19) depends(f, [X,Y,]); (%o19) [f( X,Y, )] 6 d d X f (%i0) Df3: 1/3!*dr_nabla3(f); (%o0) (d (d d d 3 f + d Y d f + f + d X d d Y d f +d d Y d f + d X d Y d f + d X d f +d d X d f + d X d Y d f )+ ( d f +d Y 3 d Y d f + d X d Y f +d d Y d f +d d Y d f + d X d Y d f + d X d f + Y d X d f +d Y d X d Y d f )+ ( d f +d X 3 d X d f + d X d f +d Y d X d f +d d X d f + d X d Y d f + d X d f + Y d X d Y f +d d X d Y d f ))/ 6 (%i) f1: ratsimp(averagem(df3, dr)); (%o1) 0 (%o) th order
5 396(1).wxm 5 / 14 (%i3) depends(f, [X,Y,]); (%o3) [f( X,Y, )]
6 396(1).wxm 6 / 14 (%i4) Df4: 1/4!*dr_nabla4(f); (%o4) (d (d (d d d 4 f + d Y d 3 f + d X d 3 f + d Y d f +d d Y d f + 3 d X d Y d f + d X d f +d d X d f + 3 d X d Y d f )+ ( d Y 3 d f +d d Y d f + d X d Y d f +d d Y d f +d d Y d f + 3 d X d Y d f + d X d Y d f + d X d Y d f +d d X d Y d f )+ ( f +d f + f +d d X 3 d d X d d X d Y d d X d f +d d X d f + 3 d X d Y d f + f + f +d f ))+ ( ( d X d Y d d X d Y d d X d Y d d Y 4 f +d d Y 3 d f + d X d Y 3 f +d d Y 3 d f +d d Y d f + d X d Y d f + d X d Y f + d X d Y 3 f +d d X d Y d f )+d ( d Y 3 d f +d d Y d f + d X d Y d f +d d Y d f +d d Y d f + 3 d X d Y d f + d X d Y d f + d X d Y d f +d d X d Y d f )+ ( d X 3 d f + Y d X d Y f +d d X d Y d f + d X d f + Y d X d f +d Y 3 d X d Y d f +d d X f + d Y d d X d Y f +d f ))+ ( ( d d X d Y d d X 4 f +d d X 3 d f + d X 3 d f +d Y f +d f + f + d X 3 d d X d d X d Y d
7 396(1).wxm 7 / 14 (%i6) f1: ratsimp(averagem(df4, dr)); (%o5) (d 4 d f 4 +4 d f +6 Y 4 d d f +6 X 4 d d X d f +6 (%o6) dr4 av d 4 f th order d4 d Y f +6 4 (%i7) depends(f, [X,Y,]); (%o7) [f( X,Y, )] (%i8) Df5: 1/5!*dr_nabla5(f)$ d Y d f +4 f )/ 4 d X d Y d Y d f + d X d f + 16 (%i30) f1: ratsimp(averagem(df5, dr)); (%o9) 0 (%o30) th order (%i31) depends(f, [X,Y,]); (%o31) [f( X,Y, )] (%i3) Df6: 1/6!*dr_nabla6(f)$ d X d Y f + d4 d X 4 f
8 396(1).wxm 8 / 14 (%i34) f1: ratsimp(averagem(df6, dr)); (%o33) (d 6 d 6 d f 6 +6 d 6 d f +15 Y 6 4 d d 6 f +15 d Y 4 d d 4 d 6 d Y d f 4 +6 d 6 d f +15 X 6 4 d d 6 d X 4 d f d 6 d X 4 d f +15 Y d 4 d 6 d X d f d 6 f +90 d X d Y4 d d 6 f )/ 70 d X d Y d (%o34) (dr 6 av ( d6 d 6 f + d 6 d Y 4 d f + d 6 f +90 d X 6 d6 d Y d 4 f + d 6 d6 d Y 6 f +15 d6 d X 4 d f + f ))/ d X d Y d d6 d X 4 d Y f + d6 d X d 4 f + 4 Vector potential of magnetic dipole (%i35) m: [m_x, m_y, m_]; (%o35) [m X,m Y,m ] d6 d X d Y 4 f + (%i36) r: [X, Y, ]; (%o36) [X,Y,] (%i37) dr: [,, d]; (%o37) [,,d] (%i38) r: r.r; (%o38) +Y +X (%i39) A: cross(m, r)/r^(3/); (%o39) [ m Y Y m, 3/ X m m X, 3/ (%i40) [m_x, m_y, m_]: [0,0,1]; (%o40) [0,0,1] (%i41) A: ev(a); (%o41) [ Y, 3/ X,0] 3/ Y m X X m Y ] 3/
9 396(1).wxm 9 / DA, nd order (%i43) DAX: 1/!*dr_nabla(A[1]); DAY: 1/!*dr_nabla(A[]); 3 d 15 Y d (%o4) ( ( + 5/ 9 Y 15 Y 3 A 3 X 15 X Y + )+d 5/ 7/ 5/ 7/ 3 Y d 15 Y d ( 5/ 3 15 Y A 7/ 5/ 7/ 15 X Y 15 X Y d / )+ ( 3 X A / / Y 5/ 7 15 X Y ))/ 7/ 5 15 X Y 7/ (%o43) ( ( 15 X Y 7 3 X d ( ) / + / + 3 d ( ) / + 15 X d ( ) / 3 Y ( ) +Y +X 5 +Y +X 7 +Y +X 15 X 3 9 X A / / )+d ( 7 15 X d ( ) / X Y ( ) / + +Y +X 5 +Y +X 7 +Y +X 15 X 3 A 15 X Y d )+ ( 7/ 5/ 7/ + 15 X Y A))/ / / / 15 X Y Y 5 / + 3 X 5/
10 396(1).wxm 10 / 14 (%i47) f1: ratsimp(averagem(dax, dr)); f1: ratsimp(averagem(day, dr)); p (%o44) ( +Y +X (( 1 Y 3 Y 3 3 X Y) d +( 9 Y +6 Y 3 9 X Y) +( 3 Y 3 Y 3 +1 X Y) ))/( 8 +( 8 Y +8 X ) 6 + ( 1 Y 4 +4 X Y +1 X 4 ) 4 +( ) +1 X 4 Y 4 +8 X 6 Y + X 8 ) (%o45) 0 p 8 Y 6 +4 X Y 4 +4 X 4 Y +8 X 6 + Y 8 +8 X Y 6 (%o46) ( +Y +X (( 1 X 3 X Y 3 X 3 ) d +( 3 X +1 X Y 3 X 3 ) +( 9 X 9 X Y +6 X 3 ) ))/( 8 +( 8 Y +8 X ) 6 + ( 1 Y 4 +4 X Y +1 X 4 ) 4 +( ) +1 X 4 Y 4 +8 X 6 Y + X 8 ) (%o47) 0 4. DA, 4th order (%i49) DA4X: 1/4!*dr_nabla4(A[1])$ DA4Y: 1/4!*dr_nabla4(A[])$ (%i53) f1: ratsimp(averagem(da4x, dr))$ f1: ratsimp(averagem(da4y, dr))$ 8 Y 6 +4 X Y 4 +4 X 4 Y +8 X 6 + Y 8 +8 X Y 6 (%o51) 35 Y ( Y 3 X +Y 4 5 X Y +3 X 4) dr4 av 18 11/ (%o53) 35 X ( Y 5 X +3 Y 4 5 X Y +X 4) dr4 av 18 11/ 4.3 DA, 6th order (%i55) DA6X: 1/6!*dr_nabla6(A[1])$ DA6Y: 1/6!*dr_nabla6(A[])$ (%i59) f1: ratsimp(averagem(da6x, dr))$ f1: ratsimp(averagem(da6y, dr))$ (%o57) (7 Y ( Y 4 75 X 4 1 Y X Y 75 X 4 + Y 6 1 X Y 4 15 X 4 Y +8 X 6 ) dr av +Y +X 15/ ) 6 )/(9 ( ) (%o59) (7 X ( Y 4 15 X 4 75 Y X Y 1 X 4 +8 Y 6 15 X Y 4 1 X 4 Y + X 6 ) dr 6 av )/(9 15/ )
11 396(1).wxm 11 / 14 5 Plots magnetic vector potential (%i60) A_XY: [ev(a[1], [=0]), ev(a[], [=0])]; (%o60) [ Y, ( Y +X ) 3/ X ] ( Y +X ) 3/ (%i61) plotdf(a_xy, [X,Y], [X,-3,3], [Y,-3,3])$ (%i6) A_X: [ev(a[1], [Y=1]), ev(a[3], [Y=1])]; (%o6) [ 1,0] ( +X +1) 3/ (%i63) plotdf(a_x, [X,], [X,-3,3], [,-3,3])$ 6 Coulomb potential (scalar) (%i64) r: [X, Y, ]; (%o64) [X,Y,] (%i65) dr: [,, d]; (%o65) [,,d] (%i66) r: r.r; (%o66) +Y +X (%i67) U: -1/sqrt(r); 1 (%o67) p +Y +X 6.1 DU, nd order (%i68) DU: 1/!*dr_nabla(U); 0 (%o68) (d d 3 d 3 Y 3 X A 3/ 5/ 5/ + 3 Y d 3 Y 3 X Y + A 5/ 3/ 5/ + A)/ / 3 X d 5/ 3 X Y + 5/ 3/ 3 X 5
12 396(1).wxm 1 / 14 (%i70) f1: ratsimp(averagem(du, dr)); (%o69) p +Y +X Y + X ( Y X ) d +( + Y X ) +( ) 6 +( 6 Y +6 X ) 4 +( ) (%o70) 0 6. DU, 4th order (%i71) DU4: 1/4!*dr_nabla4(U)$ 6 Y 4 +1 X Y +6 X 4 + Y 6 +6 X Y 4 +6 X 4 Y + X 6 (%i73) f1: ratsimp(averagem(du4, dr)); f: DU4Rav: replacer(f1,dr); (%o7) ( 8 4 +( 4 Y 4 X ) +3 Y 4 +6 X Y +3 X 4 d 4 +( 4 4 +( 16 Y 18 X ) 4 Y 4 18 X Y +6 X ( 16 X 18 Y ) +6 Y 4 18 X Y 4 X 4 ) d ( 6 X 4 Y ) +8 Y 4 4 X Y +3 X ( 18 Y 18 X ) 4 Y X Y 4 X 4 + p 3 4 +( 6 Y 4 X ) +3 Y 4 4 X Y +8 X 4 4 )/( +Y +X (8 8 + ( 3 Y +3 X ) 6 +( 48 Y X Y +48 X 4 ) 4 + ( 3 Y X Y X 4 Y +3 X 6 ) +8 Y 8 +3 X Y X 4 Y 4 +3 X 6 Y +8 X 8 )) (%o73) 7 ( 4 3 Y 3 X +Y 4 3 X Y +X 4) dr4 av 18 9/ 6.3 DU, 6th order (%i74) DU6: 1/6!*dr_nabla6(U)$ (%i76) f1: ratsimp(averagem(du6, dr))$ f: DU6Rav: replacer(f1,dr); (%o76) (( 6 15 Y 4 15 X 4 15 Y X Y 15 X 4 + Y 6 15 X Y 4 15 X 4 Y + X 6 ) dr 6 av )/(9 13/ ) 7 Plots Coulomb potential (%i77) assume(x>0,y>0,>0); (%o77) [X>0,Y>0,>0]
13 396(1).wxm 13 / 14 (%i81) UX: ev(u, [Y=0, =0]); DU4X: ev(du4rav/dr[av]^4*dra^4, [dra=., Y=0, =0]); DU6X: ev(du6rav/dr[av]^6*dra^6, [dra=., Y=0, =0]); Utot: UX+DU4X+DU6X; (%o78) 1 X (%o79) X 5 (%o80) X 7 (%o81) 1 X X X 7 (%i8) wxplotd([ux,du4x,du6x], [X,0,1], [y,-30,30], [legend, "UX", "D4UX", "D6UX"]); plotd: expression evaluates to non-numeric value somewhere in plotting range plotd: some values were clipped. plotd: expression evaluates to non-numeric value somewhere in plotting range plotd: some values were clipped. plotd: expression evaluates to non-numeric value somewhere in plotting range plotd: some values were clipped. (%t8) (%o8)
14 396(1).wxm 14 / 14 (%i83) wxplotd([ux,utot], [X,0,1], [y,-30,0], [legend, "UX", "U_{tot}"]); plotd: expression evaluates to non-numeric value somewhere in plotting range plotd: some values were clipped. plotd: expression evaluates to non-numeric value somewhere in plotting range plotd: some values were clipped. (%t83) (%o83)
Chapter 1: Precalculus Review
: Precalculus Review Math 115 17 January 2018 Overview 1 Important Notation 2 Exponents 3 Polynomials 4 Rational Functions 5 Cartesian Coordinates 6 Lines Notation Intervals: Interval Notation (a, b) (a,
More informationChapter 6: Vector Analysis
Chapter 6: Vector Analysis We use derivatives and various products of vectors in all areas of physics. For example, Newton s 2nd law is F = m d2 r. In electricity dt 2 and magnetism, we need surface and
More informationRJP Problems. Mathematica problems MM1, MM2, & MM3 are found under the Mathematical Summary document.
RJP Problems Mathematica problems MM1, MM2, & MM3 are found under the Mathematical Summary document. RJP-10: Write a Fortran program that computes the x and y coordinates of the points on a circle of arbitrary
More informationMultipole moments. November 9, 2015
Multipole moments November 9, 5 The far field expansion Suppose we have a localized charge distribution, confined to a region near the origin with r < R. Then for values of r > R, the electric field must
More informationMATRICES EXAM QUESTIONS
MATRICES EXAM QUESTIONS (Part One) Question 1 (**) The matrices A, B and C are given below in terms of the scalar constants a, b, c and d, by 2 3 A =, 1 a b 1 B =, 2 4 1 c C =. d 4 Given that A + B = C,
More informationThe Particle-Field Hamiltonian
The Particle-Field Hamiltonian For a fundamental understanding of the interaction of a particle with the electromagnetic field we need to know the total energy of the system consisting of particle and
More informationThe Helmholtz theorem at last!
Problem. The Helmholtz theorem at last! Recall in class the Helmholtz theorem that says that if if E =0 then E can be written as E = φ () if B =0 then B can be written as B = A (2) (a) Let n be a unit
More informationLinear Algebra. 1.1 Introduction to vectors 1.2 Lengths and dot products. January 28th, 2013 Math 301. Monday, January 28, 13
Linear Algebra 1.1 Introduction to vectors 1.2 Lengths and dot products January 28th, 2013 Math 301 Notation for linear systems 12w +4x + 23y +9z =0 2u + v +5w 2x +2y +8z =1 5u + v 6w +2x +4y z =6 8u 4v
More informationAnnouncements Wednesday, October 10
Announcements Wednesday, October 10 The second midterm is on Friday, October 19 That is one week from this Friday The exam covers 35, 36, 37, 39, 41, 42, 43, 44 (through today s material) WeBWorK 42, 43
More informationAnnouncements Monday, October 02
Announcements Monday, October 02 Please fill out the mid-semester survey under Quizzes on Canvas WeBWorK 18, 19 are due Wednesday at 11:59pm The quiz on Friday covers 17, 18, and 19 My office is Skiles
More informationQuartic Equation. By CH vd Westhuizen A unique Solution assuming Complex roots. Ax^4 + Bx^3 + Cx^2 + Dx + E = 0
Quartic Equation By CH vd Westhuizen A unique Solution assuming Complex roots The general Quartic is given by Ax^4 + Bx^3 + Cx^ + Dx + E = 0 As in the third order polynomial we are first going to reduce
More informationMATHEMATICS (2) Before you begin read these instructions carefully:
NATURAL SCIENCES TRIPOS Part IA Wednesday, 13 June, 2018 9:00 am to 12:00 pm MATHEMATICS (2) Before you begin read these instructions carefully: The paper has two sections, A and B. Section A contains
More information1 Equation set of 7 equations
380(4).wxm 1 / 8 (%i1) kill(all); (%o0) one (%i1) cross(a,b) := [a[2]*b[3] - a[3]*b[2], a[3]*b[1] - a[1]*b[3], a[1]*b[2] - a[2]*b[1]]; (%o1) cross( a,b ):=[a 2 b 3 a 3 b 2,a 3 b 1 a 1 b 3,a 1 b 2 a 2 b
More informationVECTOR SPACES & SUBSPACES
VECTOR SPACES & SUBSPACES Definition: The real Vector Space " V " is the set of all entities called "vectors" with real entries that satisfy two closure properties and obey a set of eight rules. If "x"
More information11 a 12 a 21 a 11 a 22 a 12 a 21. (C.11) A = The determinant of a product of two matrices is given by AB = A B 1 1 = (C.13) and similarly.
C PROPERTIES OF MATRICES 697 to whether the permutation i 1 i 2 i N is even or odd, respectively Note that I =1 Thus, for a 2 2 matrix, the determinant takes the form A = a 11 a 12 = a a 21 a 11 a 22 a
More informationCh 4 Differentiation
Ch 1 Partial fractions Ch 6 Integration Ch 2 Coordinate geometry C4 Ch 5 Vectors Ch 3 The binomial expansion Ch 4 Differentiation Chapter 1 Partial fractions We can add (or take away) two fractions only
More informationDIPOLES III. q const. The voltage produced by such a charge distribution is given by. r r'
DIPOLES III We now consider a particularly important charge configuration a dipole. This consists of two equal but opposite charges separated by a small distance. We define the dipole moment as p lim q
More informationMATH 105: PRACTICE PROBLEMS FOR CHAPTER 3: SPRING 2010
MATH 105: PRACTICE PROBLEMS FOR CHAPTER 3: SPRING 010 INSTRUCTOR: STEVEN MILLER (SJM1@WILLIAMS.EDU Question 1 : Compute the partial derivatives of order 1 order for: (1 f(x, y, z = e x+y cos(x sin(y. Solution:
More information( ) is called the dependent variable because its
page 1 of 16 CLASS NOTES: 3 8 thru 4 3 and 11 7 Functions, Exponents and Polynomials 3 8: Function Notation A function is a correspondence between two sets, the domain (x) and the range (y). An example
More informationMath M111: Lecture Notes For Chapter 3
Section 3.1: Math M111: Lecture Notes For Chapter 3 Note: Make sure you already printed the graphing papers Plotting Points, Quadrant s signs, x-intercepts and y-intercepts Example 1: Plot the following
More informationMathsGeeks. Everything You Need to Know A Level Edexcel C4. March 2014 MathsGeeks Copyright 2014 Elite Learning Limited
Everything You Need to Know A Level Edexcel C4 March 4 Copyright 4 Elite Learning Limited Page of 4 Further Binomial Expansion: Make sure it starts with a e.g. for ( x) ( x ) then use ( + x) n + nx + n(n
More informationMATHEMATICS 132 Applied Mathematics 1A (Engineering) EXAMINATION
MATHEMATICS 132 Applied Mathematics 1A (Engineering) EXAMINATION DURATION: 3 HOURS 26TH MAY 2011 MAXIMUM MARKS: 100 LECTURERS: PROF J. VAN DEN BERG AND DR J. M. T. NGNOTCHOUYE EXTERNAL EXAMINER: DR K.
More informationarxiv: v4 [math.nt] 31 Jul 2017
SERIES REPRESENTATION OF POWER FUNCTION arxiv:1603.02468v4 [math.nt] 31 Jul 2017 KOLOSOV PETRO Abstract. In this paper described numerical expansion of natural-valued power function x n, in point x = x
More informationChapter. Algebra techniques. Syllabus Content A Basic Mathematics 10% Basic algebraic techniques and the solution of equations.
Chapter 2 Algebra techniques Syllabus Content A Basic Mathematics 10% Basic algebraic techniques and the solution of equations. Page 1 2.1 What is algebra? In order to extend the usefulness of mathematical
More informationCM2104: Computational Mathematics General Maths: 2. Algebra - Factorisation
CM204: Computational Mathematics General Maths: 2. Algebra - Factorisation Prof. David Marshall School of Computer Science & Informatics Factorisation Factorisation is a way of simplifying algebraic expressions.
More informationWelcome back to PHY 3305
Welcome back to PHY 3305 Today s Lecture: Hydrogen Atom Part I John von Neumann 1903-1957 One-Dimensional Atom To analyze the hydrogen atom, we must solve the Schrodinger equation for the Coulomb potential
More informationENGI 4430 PDEs - d Alembert Solutions Page 11.01
ENGI 4430 PDEs - d Alembert Solutions Page 11.01 11. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More informationTwo dimensional oscillator and central forces
Two dimensional oscillator and central forces September 4, 04 Hooke s law in two dimensions Consider a radial Hooke s law force in -dimensions, F = kr where the force is along the radial unit vector and
More information7.5 Rationalizing Denominators and Numerators of Radical Expressions
440 CHAPTER Rational Exponents, Radicals, and Complex Numbers 86. Find the area and perimeter of the trapezoid. (Hint: The area of a trapezoid is the product of half the height 6 meters and the sum of
More informationAE/ME 339. Computational Fluid Dynamics (CFD) K. M. Isaac. Momentum equation. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.
AE/ME 339 Computational Fluid Dynamics (CFD) 9//005 Topic7_NS_ F0 1 Momentum equation 9//005 Topic7_NS_ F0 1 Consider the moving fluid element model shown in Figure.b Basis is Newton s nd Law which says
More informationRadiation Integrals and Auxiliary Potential Functions
Radiation Integrals and Auxiliary Potential Functions Ranga Rodrigo June 23, 2010 Lecture notes are fully based on Balanis [?]. Some diagrams and text are directly from the books. Contents 1 The Vector
More informationSOLUTIONS FOR PROBLEMS 1-30
. Answer: 5 Evaluate x x + 9 for x SOLUTIONS FOR PROBLEMS - 0 When substituting x in x be sure to do the exponent before the multiplication by to get (). + 9 5 + When multiplying ( ) so that ( 7) ( ).
More informationAN INTRODUCTION TO CURVILINEAR ORTHOGONAL COORDINATES
AN INTRODUCTION TO CURVILINEAR ORTHOGONAL COORDINATES Overview Throughout the first few weeks of the semester, we have studied vector calculus using almost exclusively the familiar Cartesian x,y,z coordinate
More informationChapter 3 - Vector Calculus
Chapter 3 - Vector Calculus Gradient in Cartesian coordinate system f ( x, y, z,...) dr ( dx, dy, dz,...) Then, f f f f,,,... x y z f f f df dx dy dz... f dr x y z df 0 (constant f contour) f dr 0 or f
More informationmultiply both sides of eq. by a and projection overlap
Fourier Series n x n x f xa ancos bncos n n periodic with period x consider n, sin x x x March. 3, 7 Any function with period can be represented with a Fourier series Examples (sawtooth) (square wave)
More informationMath 101 Study Session Spring 2016 Test 4 Chapter 10, Chapter 11 Chapter 12 Section 1, and Chapter 12 Section 2
Math 101 Study Session Spring 2016 Test 4 Chapter 10, Chapter 11 Chapter 12 Section 1, and Chapter 12 Section 2 April 11, 2016 Chapter 10 Section 1: Addition and Subtraction of Polynomials A monomial is
More informationE = 1 c. where ρ is the charge density. The last equality means we can solve for φ in the usual way using Coulomb s Law:
Photons In this section, we want to quantize the vector potential. We will work in the Coulomb gauge, where A = 0. In gaussian units, the fields are defined as follows, E = 1 c A t φ, B = A, where A and
More informationGame Engineering: 2D
Game Engineering: 2D CS420-2010F-06 2D Math David Galles Department of Computer Science University of San Francisco 06-0: Back to Basics A Vector is a displacement Vector has both direction and length
More informationChapter-2 Relations and Functions. Miscellaneous
1 Chapter-2 Relations and Functions Miscellaneous Question 1: The relation f is defined by The relation g is defined by Show that f is a function and g is not a function. The relation f is defined as It
More informationChapter 6: Rational Expr., Eq., and Functions Lecture notes Math 1010
Section 6.1: Rational Expressions and Functions Definition of a rational expression Let u and v be polynomials. The algebraic expression u v is a rational expression. The domain of this rational expression
More informationSolutions: Homework 5
Ex. 5.1: Capacitor Solutions: Homework 5 (a) Consider a parallel plate capacitor with large circular plates, radius a, a distance d apart, with a d. Choose cylindrical coordinates (r,φ,z) and let the z
More informationSolutions to the Calculus and Linear Algebra problems on the Comprehensive Examination of January 28, 2011
Solutions to the Calculus and Linear Algebra problems on the Comprehensive Examination of January 8, Solutions to Problems 5 are omitted since they involve topics no longer covered on the Comprehensive
More informationl=0 The expansion coefficients can be determined, for example, by finding the potential on the z-axis and expanding that result in z.
Electrodynamics I Exam - Part A - Closed Book KSU 15/11/6 Name Electrodynamic Score = 14 / 14 points Instructions: Use SI units. Where appropriate, define all variables or symbols you use, in words. Try
More informationThe Magnetic Field
4-50 4. The Magnetic Field 4.9 Mathematical Field Theory A field is defined as a special domain in which a physical quantity is depicted as a function of space and time. Mathematics supports the analytical
More informationMath 4377/6308 Advanced Linear Algebra
2.3 Composition Math 4377/6308 Advanced Linear Algebra 2.3 Composition of Linear Transformations Jiwen He Department of Mathematics, University of Houston jiwenhe@math.uh.edu math.uh.edu/ jiwenhe/math4377
More informationLinear Equations and Matrix
1/60 Chia-Ping Chen Professor Department of Computer Science and Engineering National Sun Yat-sen University Linear Algebra Gaussian Elimination 2/60 Alpha Go Linear algebra begins with a system of linear
More informationg(2, 1) = cos(2π) + 1 = = 9
1. Let gx, y 2x 2 cos2πy 2 + y 2. You can use the fact that Dg2, 1 [8 2]. a Find an equation for the tangent plane to the graph z gx, y at the point 2, 1. There are two key parts to this problem. The first,
More informationVector Potential for the Magnetic Field
Vector Potential for the Magnetic Field Let me start with two two theorems of Vector Calculus: Theorem 1: If a vector field has zero curl everywhere in space, then that field is a gradient of some scalar
More informationSummary: Curvilinear Coordinates
Physics 2460 Electricity and Magnetism I, Fall 2007, Lecture 10 1 Summary: Curvilinear Coordinates 1. Summary of Integral Theorems 2. Generalized Coordinates 3. Cartesian Coordinates: Surfaces of Constant
More informationCreated by T. Madas VECTOR OPERATORS. Created by T. Madas
VECTOR OPERATORS GRADIENT gradϕ ϕ Question 1 A surface S is given by the Cartesian equation x 2 2 + y = 25. a) Draw a sketch of S, and describe it geometrically. b) Determine an equation of the tangent
More informationChapter 4: Partial differentiation
Chapter 4: Partial differentiation It is generally the case that derivatives are introduced in terms of functions of a single variable. For example, y = f (x), then dy dx = df dx = f. However, most of
More informationChapter 2: Unconstrained Extrema
Chapter 2: Unconstrained Extrema Math 368 c Copyright 2012, 2013 R Clark Robinson May 22, 2013 Chapter 2: Unconstrained Extrema 1 Types of Sets Definition For p R n and r > 0, the open ball about p of
More information(x 3)(x + 5) = (x 3)(x 1) = x + 5. sin 2 x e ax bx 1 = 1 2. lim
SMT Calculus Test Solutions February, x + x 5 Compute x x x + Answer: Solution: Note that x + x 5 x x + x )x + 5) = x )x ) = x + 5 x x + 5 Then x x = + 5 = Compute all real values of b such that, for fx)
More informationThe Derivative. Appendix B. B.1 The Derivative of f. Mappings from IR to IR
Appendix B The Derivative B.1 The Derivative of f In this chapter, we give a short summary of the derivative. Specifically, we want to compare/contrast how the derivative appears for functions whose domain
More informationBrief Review of Vector Algebra
APPENDIX Brief Review of Vector Algebra A.0 Introduction Vector algebra is used extensively in computational mechanics. The student must thus understand the concepts associated with this subject. The current
More informationGrade 9 Linear Equations in Two Variables
ID : cn-9-linear-equations-in-two-variables [1] Grade 9 Linear Equations in Two Variables For more such worksheets visit www.edugain.com Answer the questions (1) In the graph of the linear equation 4x
More informationJUST THE MATHS UNIT NUMBER 1.5. ALGEBRA 5 (Manipulation of algebraic expressions) A.J.Hobson
JUST THE MATHS UNIT NUMBER 1.5 ALGEBRA 5 (Manipulation of algebraic expressions) by A.J.Hobson 1.5.1 Simplification of expressions 1.5.2 Factorisation 1.5.3 Completing the square in a quadratic expression
More informationSTRAND J: TRANSFORMATIONS, VECTORS and MATRICES
Mathematics SKE, Strand J STRAND J: TRANSFORMATIONS, VECTORS and MATRICES J4 Matrices Text Contents * * * * Section J4. Matrices: Addition and Subtraction J4.2 Matrices: Multiplication J4.3 Inverse Matrices:
More informationb n x n + b n 1 x n b 1 x + b 0
Math Partial Fractions Stewart 7.4 Integrating basic rational functions. For a function f(x), we have examined several algebraic methods for finding its indefinite integral (antiderivative) F (x) = f(x)
More informationSection September 6, If n = 3, 4, 5,..., the polynomial is called a cubic, quartic, quintic, etc.
Section 2.1-2.2 September 6, 2017 1 Polynomials Definition. A polynomial is an expression of the form a n x n + a n 1 x n 1 + + a 1 x + a 0 where each a 0, a 1,, a n are real numbers, a n 0, and n is a
More informationKevin James. MTHSC 206 Section 12.5 Equations of Lines and Planes
MTHSC 206 Section 12.5 Equations of Lines and Planes Definition A line in R 3 can be described by a point and a direction vector. Given the point r 0 and the direction vector v. Any point r on the line
More informationarxiv:math.ca/ v2 17 Jul 2000
NECESSARY AND SUFFICIENT CONDITIONS FOR DIFFERENTIABILITY OF A FUNCTION OF SEVERAL VARIABLES. R.P. VENKATARAMAN, #1371, 13'TH MAIN ROAD, II STAGE, FIRST PHASE, B.T.M. LAYOUT, BANGALORE 560 076.
More informationLearn how to use Desmos
Learn how to use Desmos Maclaurin and Taylor series 1 Go to www.desmos.com. Create an account (click on bottom near top right of screen) Change the grid settings (click on the spanner) to 1 x 3, 1 y 12
More informationELECTROSTATIC POTENTIAL FROM DISTRIBUTED CHARGES
MISN--147 ELECTROSTATIC POTENTIAL FROM DISTRIBUTED CHARGES z ELECTROSTATIC POTENTIAL FROM DISTRIBUTED CHARGES by J. Kovacs and P. Signell Michigan State University 1. Introduction..............................................
More informationMath Euler Cauchy Equations
1 Math 371 - Euler Cauchy Equations Erik Kjær Pedersen November 29, 2005 We shall now consider the socalled Euler Cauchy equation x 2 y + axy + by = 0 where a and b are constants. To solve this we put
More informationCHAPTER 11 RADIATION 4/13/2017. Outlines. 1. Electric Dipole radiation. 2. Magnetic Dipole Radiation. 3. Point Charge. 4. Synchrotron Radiation
CHAPTER 11 RADIATION Outlines 1. Electric Dipole radiation 2. Magnetic Dipole Radiation 3. Point Charge Lee Chow Department of Physics University of Central Florida Orlando, FL 32816 4. Synchrotron Radiation
More informationδ Substituting into the differential equation gives: x i+1 x i δ f(t i,x i ) (3)
Solving Differential Equations Numerically Differential equations are ubiquitous in Physics since the laws of nature often take on a simple form when expressed in terms of infinitesimal changes of the
More information4 ENERGY Last Updated: July 18, 2012
Problem List 4.1 Drawing equipotentials 4.2 Particle in a 2D well 4.3 Computing line integrals of the magnetic field 4.4 Shrinking an orbit 4.5 You re a collision investigator! 4.6 Choose an initial velocity
More information3/31/ Product of Inertia. Sample Problem Sample Problem 10.6 (continue)
/1/01 10.6 Product of Inertia Product of Inertia: I xy = xy da When the x axis, the y axis, or both are an axis of symmetry, the product of inertia is zero. Parallel axis theorem for products of inertia:
More informationEcon 101A Problem Set 1 Solution
Econ 101A Problem Set 1 Solution Problem 1. Univariate unconstrained maximization. (10 points) Consider the following maximization problem: max x f(x; x 0)=exp( (x x 0 ) 2 ) 1. Write down the first order
More informationVectors. September 2, 2015
Vectors September 2, 2015 Our basic notion of a vector is as a displacement, directed from one point of Euclidean space to another, and therefore having direction and magnitude. We will write vectors in
More informationTutorial Divergence. (ii) Explain why four of these integrals are zero, and calculate the other two.
(1) Below is a graphical representation of a vector field v with a z-component equal to zero. (a) Draw a box somewhere inside this vector field. The box is 3-dimensional. To make things easy, it is a good
More informationA Brief Revision of Vector Calculus and Maxwell s Equations
A Brief Revision of Vector Calculus and Maxwell s Equations Debapratim Ghosh Electronic Systems Group Department of Electrical Engineering Indian Institute of Technology Bombay e-mail: dghosh@ee.iitb.ac.in
More information1 Expectation of a continuously distributed random variable
OCTOBER 3, 204 LECTURE 9 EXPECTATION OF A CONTINUOUSLY DISTRIBUTED RANDOM VARIABLE, DISTRIBUTION FUNCTION AND CHANGE-OF-VARIABLE TECHNIQUES Expectation of a continuously distributed random variable Recall
More informationB = 0. E = 1 c. E = 4πρ
Photons In this section, we will treat the electromagnetic field quantum mechanically. We start by recording the Maxwell equations. As usual, we expect these equations to hold both classically and quantum
More informationMTH MTH Lecture 6. Yevgeniy Kovchegov Oregon State University
MTH 306 0 MTH 306 - Lecture 6 Yevgeniy Kovchegov Oregon State University MTH 306 1 Topics Lines and planes. Systems of linear equations. Systematic elimination of unknowns. Coe cient matrix. Augmented
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationFermat's Last Theorem, Revisited. Ernest R. Lucier
Fermat's Last Theorem, Revisited Ernest R. Lucier Abstract A non-zero integer solution set is derived for C 2 = A 2 + B 2 where n = 2. The solution set for n = 2 is applied to n > 2. The result is a proof
More information26.2. Cauchy-Riemann Equations and Conformal Mapping. Introduction. Prerequisites. Learning Outcomes
Cauchy-Riemann Equations and Conformal Mapping 26.2 Introduction In this Section we consider two important features of complex functions. The Cauchy-Riemann equations provide a necessary and sufficient
More informationVorticity in natural coordinates
Vorticity in natural coordinates (see Holton pg 95, section 4.2.) Let s consider the vertical vorticity component only, i.e. ζ kˆ ω, we have ω u dl kˆ ω lim --- lim ----------------- curve is in xy plane
More information1. FUNDAMENTAL CONCEPTS AND MATH REVIEW
1. FUNDAMENTAL CONCEPTS AND MATH REVIEW 1.1. Introduction Here we provide for your reading pleasure a review of some of the math concepts used in part of this course. Most of this falls under the heading
More informationLinear Equations. Find the domain and the range of the following set. {(4,5), (7,8), (-1,3), (3,3), (2,-3)}
Linear Equations Domain and Range Domain refers to the set of possible values of the x-component of a point in the form (x,y). Range refers to the set of possible values of the y-component of a point in
More informationLinear algebra I Homework #1 due Thursday, Oct. 5
Homework #1 due Thursday, Oct. 5 1. Show that A(5,3,4), B(1,0,2) and C(3, 4,4) are the vertices of a right triangle. 2. Find the equation of the plane that passes through the points A(2,4,3), B(2,3,5),
More information3 Applications of partial differentiation
Advanced Calculus Chapter 3 Applications of partial differentiation 37 3 Applications of partial differentiation 3.1 Stationary points Higher derivatives Let U R 2 and f : U R. The partial derivatives
More informationf n+2 f n = 2 = = 40 H z
Coordinator: Dr.I.M.Nasser Monday, May 18, 2015 Page: 1 Q1. Q2. One end of a stretched string vibrates with a period of 1.5 s. This results in a wave propagating at a speed of 8.0 m/s along the string.
More informationChapter 2. Matrix Arithmetic. Chapter 2
Matrix Arithmetic Matrix Addition and Subtraction Addition and subtraction act element-wise on matrices. In order for the addition/subtraction (A B) to be possible, the two matrices A and B must have the
More informationRydberg States of ThIII and UV
Rydberg States of ThIII and UV W. R. Johnson Department of Physics, Notre Dame University June 25, 2003 Abstract We examine the Rydberg spectrum of the two actinide ions Th 2+ and U 5+. The states of interest
More informationPHYS 281: Midterm Exam
PHYS 28: Midterm Exam October 28, 200, 8:00-9:20 Last name (print): Initials: No calculator or other aids allowed PHYS 28: Midterm Exam Instructor: B. R. Sutherland Date: October 28, 200 Time: 8:00-9:20am
More information1.2 Graphs and Lines. Cartesian Coordinate System
1.2 Graphs and Lines Cartesian Coordinate System Note that there is a one-to-one correspondence between the points in a plane and the elements in the set of all ordered pairs (a, b) of real numbers. Graphs
More informationHomework 1/Solutions. Graded Exercises
MTH 310-3 Abstract Algebra I and Number Theory S18 Homework 1/Solutions Graded Exercises Exercise 1. Below are parts of the addition table and parts of the multiplication table of a ring. Complete both
More informationPoisson s equation. Lab 1. Poisson s equation in two dimensions
Lab 1 Poisson s equation Suppose that we want to describe the distribution of heat throughout a region Ω Let h(x) represent the temperature on the boundary of Ω ( Ω), and let g(x) represent the initial
More informationPaul Schrimpf. October 18, UBC Economics 526. Unconstrained optimization. Paul Schrimpf. Notation and definitions. First order conditions
Unconstrained UBC Economics 526 October 18, 2013 .1.2.3.4.5 Section 1 Unconstrained problem x U R n F : U R. max F (x) x U Definition F = max x U F (x) is the maximum of F on U if F (x) F for all x U and
More informationOR MSc Maths Revision Course
OR MSc Maths Revision Course Tom Byrne School of Mathematics University of Edinburgh t.m.byrne@sms.ed.ac.uk 15 September 2017 General Information Today JCMB Lecture Theatre A, 09:30-12:30 Mathematics revision
More informationQuaternion Dynamics, Part 1 Functions, Derivatives, and Integrals. Gary D. Simpson. rev 00 Dec 27, 2014.
Quaternion Dynamics, Part 1 Functions, Derivatives, and Integrals Gary D. Simpson gsim100887@aol.com rev 00 Dec 27, 2014 Summary Definitions are presented for "quaternion functions" of a quaternion. Polynomial
More informationClass 9 Linear Equations in Two Variables
ID : in-9-linear-equations-in-two-variables [1] Class 9 Linear Equations in Two Variables For more such worksheets visit www.edugain.com Answer the questions (1) A telecom operator charges Rs. 1.1 for
More informationMath Precalculus I University of Hawai i at Mānoa Spring
Math 135 - Precalculus I University of Hawai i at Mānoa Spring - 2013 Created for Math 135, Spring 2008 by Lukasz Grabarek and Michael Joyce Send comments and corrections to lukasz@math.hawaii.edu Contents
More informationSOIL MECHANICS AND PLASTIC ANALYSIS OR LIMIT DESIGN*
157 SOIL MECHANICS AND PLASTIC ANALYSIS OR LIMIT DESIGN* BY D. C. DRUCKER and W. PRAGER Brown University 1. Introduction. Problems of soil mechanics involving stability of slopes, bearing capacity of foundation
More information7.6 The Inverse of a Square Matrix
7.6 The Inverse of a Square Matrix Copyright Cengage Learning. All rights reserved. What You Should Learn Verify that two matrices are inverses of each other. Use Gauss-Jordan elimination to find inverses
More informationENERGY IN ELECTROSTATICS
ENERGY IN ELECTROSTATICS We now turn to the question of energy in electrostatics. The first question to consider is whether or not the force is conservative. You will recall from last semester that a conservative
More informationComparative Statics. Autumn 2018
Comparative Statics Autumn 2018 What is comparative statics? Contents 1 What is comparative statics? 2 One variable functions Multiple variable functions Vector valued functions Differential and total
More information