4.2 Final exam

  4.2.1 Problems listing
  4.2.2 Problem 1
  4.2.3 Problem 2
  4.2.4 Problem 3
  4.2.5 Problem 4
  4.2.6 Problem 5
  4.2.7 Problem 6
  4.2.8 Problem 7
  4.2.9 Problem 8
  4.2.10 Problem 9

4.2.1 Problems listing

PDF

PDF (letter size)

PDF (legal size)

4.2.2 Problem 1

Solve the following system of equations and write the solution as a parametric vector formx+2y3z=52x+y3z=13x+y=8

Solution

In matrix form Ax=b, the above system is(1)[123213110][xyz]=[5138] The augmented matrix is[1235213131108] R2=R22R1[123503331108] R3=R3+R1[123503330333] R3=R3+R2[123503330000] Hence original system (1) becomes(2)[123033000][xyz]=[530] The above shows that z=t is a free variable and x,y are basic variables. Second row gives 3y+3t=3 or y+t=1 or y=t1. First row gives x+2y3t=5 or x=5+3t2(t1) or x=7+t. Hence the solution is[xyz]=[7+tt1t]=[710]+[ttt]=[710]+t[111]

The above is the solution in parametric vector form. For any value of the parameter t, a solution exist.

4.2.3 Problem 2

Compute the determinant using a cofactor expansionA=[150241020] solution

Expanding along the last row since it has most number of zeros gives (only the element A(3,2)0)det(A)=(1)3+2(2)|1021|=2|1021|=2(1)=2

4.2.4 Problem 3

   4.2.4.1 Part a
   4.2.4.2 Part b
   4.2.4.3 Part c

Let A=[134462376]u=[334]

a

is u in NullSpace of A ? Justify your answer.

b

Is u in columnspace of A ? Justify your answer.

c

Determine the rank A and the Nullity of A. Show your work

solution

4.2.4.1 Part a

For an m×n matrix, the solution set corresponding to Ax=0 is called the NullSpace(A). Therefore we need to first find this solution set by solving (1)[134462376][xyz]=[000] The augmented matrix is[134046203760] R2=R2+4R1[1340061803760] R3=R3+3R1[1340061800260] R3=R313R2[1340061800000] Hence (1) becomes(2)[1340618000][xyz]=[000] The above shows that z=t is free variable and x,y are basic variables. Second row gives 6y18t=0 or y=3t. First row gives x3y4t=0 or x=3(3t)+4t or x=5t. Hence the solution is[xyz]=[5t3tt]=t[531]

Now we are ready to answer the question if u=[334] is in the NullSpace(A). In other words, does there exist t which make t[531]=[334]. It is clear there is no such t. To show this, looking at the second row, it says 3t=3 or t=1. But third row says t=4. Therefore there is no t which makes u in NullSpace(A). Hence u is not in NullSpace(A)

4.2.4.2 Part b

The columnspace of A is the set of all linear combinations of the columns of A. The basis for the columnspace are columns of A that correspond to the pivot columns are doing the above REF.  From part A we found that column 1,2 are the pivot columns. Hence the basis of columnspace of A are{[143],[367]} Hence the columnspace of A is two dimensional subspace of R3. To find if u is in columnspace of A, we need to find if there exists a linear combination of these basis which gives u. Therefore we need to solvec1[143]+c2[367]=[334] For c1,c2 to see if a solution exist. In matrix form the above becomes(1)[134637][c1c2]=[334] The augmented matrix is[133463374] R2=R2+4R1[1330615374] R3=R3+3R1[1330615025] R3=R313R2[1330615000] R2=R26[1330152000] R1=R1+3R2[10920152000] Hence (1) becomes[100100][c1c2]=[92520] Therefore the solution is[c1c2]=[9252] Therefore a solution exists. This means 92[143]52[367]=[334] This means u is in the columnspace of A.

4.2.4.3 Part c

The rank of A is the dimension of the columnspace of A. Which is the same as the number of pivot columns found. In this case, it is 2 as found in part b above. Hence rank(A)=2.

Nullity of A is the dimension of the nullspace of A. From part (a) we found that the nullspace of A is given by the one parameter vector t[531]. Hence the dimension is 1. It is the number of the free variables. Therefore Nullity of A=1.  To verify this, we can use the rank-nullity theorem, which says for a matrix m×n, rank(A)+nullity(A)=n Since n=3 and since rank(A)=2 then nullity(A)=1.

4.2.5 Problem 4

   4.2.5.1 Part a
   4.2.5.2 Part b

a

Using the definition, verify that the given transformation is linear transformation T:C2(I)C0(I) defined by T(y)=y+y

b

Find the kernel of T

solution

4.2.5.1 Part a

The transformation T is linear if

  1. T(u+v)=T(u)+T(v) for all u,v in C2(I)

  2. T(cu)=cT(u) for all scalars c and u in C2(I)

To show property 1:T(u+v)=(u+v)+(u+v) By linearity of second derivatives (and since u,v are in C2(I)) the above becomesT(u+v)=(u+v)+(u+v)=u+u+v+v=(u+u)+(v+v)

But (u+u)=T(u) and (v+v)=T(v), Hence the above becomesT(u+v)=T(u)+T(v) To show property 2: T(cu)=(cu)+(cu) But since c is a scalar, we can move it outside the derivative and the above becomesT(cu)=c(u)+cu=c(u+u)

But u+u=T(u). Hence the above becomesT(cu)=cT(u) Both properties are satisfied. Hence T is linear transformation.

4.2.5.2 Part b

The kernel of T:VW is ker(T)={uV:Tu=0}. In this case V=C2(I) and W=C0(I). Hence we need to find all u, such that T(u)=0. Which is the same as saying all u which satisfies u+u=0. Hence ker(T) is the solution of this ode.

This is linear constant coefficient ode. The characteristic equation is λ2+1=0. The solutions are λ=±i. Hence the basis functions are {eix,eix} (assuming the independent variable is x), or using Euler relation {cosx,sinx}. Therefore the solution is linear combination of these basis given by u=c1cosx+c2sinx where c1,c2 are arbitrary constants. Henceker(T)={u:u=c1cosx+c2sinx}

4.2.6 Problem 5

Solve (y+3x2)dx+xdy=0

solution

Writing the ODE as Mdx+Ndy=0 Where M=y+3x2,N=x. Checking if the ODE is exactMy=1Nx=1

Since My=Nx, then it is exact in some region R. Let there exists constant function Φ(x,y)=c which satisfies(1)Φx=M=y+3x2(2)Φy=N=x

For all (x,y) in R.  Integrating (1) w.r.t. x givesΦxdx=y+3x2dx(3)Φ=yx+x3+g(y)

Taking derivative of the above w.r.t. y gives(4)Φy=x+g(y) Comparing (4) and (2) givesx+g(y)=xg(y)=0

Hence g(y)=c1 a constant. Substituting this in (3) givesΦ=yx+x3+c1 But Φ=c. Combining the constants c,c1 into one constant, say C, the above becomesC=yx+x3 Solving for y givesyx=Cx3 For x0y=Cx3x

4.2.7 Problem 6

Using the method of undermined coefficients, find the general solution of the given differential equation(1)yy2y=ex+2cosx solution

The solution is y=yh+yp Where yh is the solution to the homogenous ode yy2y=0 and yp is a particular solution to the ode. We start by finding yh. Since this is linear second order with constant coefficient, then the characteristic equation method is used. The characteristic equation for yy2y=0 is λ2λ2=0(λ2)(λ+1)=0

Hence the roots are λ1=2,λ2=1. Therefore the basis set of solutions for yh is the set (2){e2x,ex} And yh is linear combination of these basis. Therefore(3)yh=c1e2x+c2ex Looking at RHS of (1) shows it is linear combination of basis [ex,cosx]. For each basis in this list, we generate all possible derivatives. Which gives (ignoring sign changes and any leading constants as they will be parts of the unknowns to be found later on). This results in the following list(4)[{ex},{cosx,sinx}] Now we compare each basis in (2) with each basis in (4) to see if there is any duplication. We see that ex is in (4) as well in (2). We now multiply ex in (4) by an extra x and obtain new list(4A)[{xex},{cosx,sinx}] We repeat this process again, checking if (2) still has any duplication in (4A). There are no duplication now. Hence the trial solution is linear combination of the basis in (4A). Which gives(5)yp=Axex+Bcosx+Csinx To determine A,B,C, we substitute yp back in the ODE (1) and solve for these unknowns by comparing terms. (6)yp=AexAxexBsinx+Ccosxyp=AexAex+AxexBcosxCsinx(7)=2Aex+AxexBcosxCsinx

Substituting (5,6,7) into the ODE (1) gives

(2Aex+AxexBcosxCsinx)(AexAxexBsinx+Ccosx)2(Axex+Bcosx+Csinx)=ex+2cosx2Aex+AxexBcosxCsinxAex+Axex+BsinxCcosx2Axex2Bcosx2Csinx=ex+2cosx

Which simplifies to3Aex3Bcosx3Csinx+BsinxCcosx=ex+2cosx3Aex+cosx(3BC)+sinx(3C+B)=ex+2cosx

Comparing terms on each side gives 3 equations to solve for A,B,C3A=13BC=23C+B=0

First equation gives A=13. Multiplying second equation by 3 and adding the result to third equation gives 9B+3C=63C+B=0

Adding gives9B+B=610B=6B=610=35

From 3BC=2 we now find 3(610)C=2, or C=15. Hence the particular solution (5) becomes(8)yp=13xex35cosx15sinx Substituting (8) and (3) in y=yh+yp gives the final solution asy=c1e2x+c2ex13xex35cosx15sinx

4.2.8 Problem 7

Use the Laplace transform to solve the given initial-value problems. You can use the table of transformationy+y=e2ty(0)=0y(0)=1

solution

Taking the Laplace transform of both sides of y+y=e2t gives (using linearity)(1)L(y)+L(y)=L(e2t) Assuming L(y)=Y(s), and using the property thatL(y)=s2L(y)sy(0)y(0) And from table 10.2.1 L(e2t)=1s2,s>2, then the ode becomes(s2L(y)sy(0)y(0))+L(y)=1s2(s2Ys(0)1)+Y=1s2s2Y1+Y=1s2Y(s2+1)=1s2+1(1A)Y=1(s2)(s2+1)+1s2+1

Using partial fractions on the first term in the RHS above gives1(s2)(s2+1)=As2+Bs+Cs2+1=A(s2+1)+(Bs+C)(s2)(s2)(s2+1)=As2+A+(Bs22Bs+Cs2C)(s2)(s2+1)=As2+A+Bs22Bs+Cs2C(s2)(s2+1)=s2(A+B)+s(C2B)+(A2C)(s2)(s2+1)

Therefore1=s2(A+B)+s(C2B)+(A2C) Comparing terms givesA+B=0C2B=0A2C=1

In matrix form the above is(2)[110021102][ABC]=[001] The augmented matrix is[110002101021] R3=R3R1[110002100121] R3=2R3[110002100242] R3=R3R2[110002100052] R3=15R3,R2=12R2 gives[11000112000125] R2=R2+12R3[11000101500125] R1=R1R2[100150101500125] Hence (2) becomes(3)[100010001][ABC]=[151525] Therefore, since now in RREF form, the solution is [ABC]=[151525] Hence1(s2)(s2+1)=As2+Bs+Cs2+1=151s2+15s25s2+1=151s215s+2s2+1(4)=151s215ss2+1251s2+1

Substituting (4) back in (1A) givesY=151s215ss2+1251s2+1+1s2+1(5)=151s215ss2+1+351s2+1

Now we will use tables to do the inverse Laplace transform. From table 10.2.1L1(1s2)=e2ts>2L1(ss2+1)=costs>0L1(1s2+1)=sints>0

Using this result in (5) and since L1(Y(s))=y(t) then (5) becomesy(t)=15e2t15cost+35sint

4.2.9 Problem 8

Find a series solution in powers of x of the differential equation y+x2y+y=0

solution

Let the solution be (1)y(x)=n=0anxn Theny(x)=n=0nanxn1(2)=n=1nanxn1

Andy(x)=n=1n(n1)anxn2(3)=n=2n(n1)anxn2

Substituting (1,2,3) into the given ODE givesn=2n(n1)anxn2+x2n=1nanxn1+n=0anxn=0(3A)n=2n(n1)anxn2+n=1nanxn+1+n=0anxn=0

Now we make all powers of x the same by rewriting n=2n(n1)anxn2=n=0(n+2)(n+1)an+2xn and n=1nanxn+1=n=2(n1)an1xn. The way the above is done is by using the rule: When adding a value to the summation index n inside the sum, then we must at same time subtract the same value from the starting index n.

Hence (3A) now becomesn=0(n+2)(n+1)an+2xn+n=2(n1)an1xn+n=0anxn=0 To be able to compare coefficients of x, we expand up to n=1 the sums in order to make all sums start from n=2. This gives(2)(1)a2+(1+2)(1+1)a3x+n=2(n+2)(n+1)an+2xn+n=2(n1)an1xn+(a0+a1x)+n=2anxn=0(2a2+a0)+x(6a3+a1)+(n=2(n+2)(n+1)an+2+n=2(n1)an1+n=2an)xn=0

Now we are to compare coefficients on each power of x. The above gives the three equations(4)2a2+a0=06a3+a1=0(n+2)(n+1)an+2+(n1)an1+an=0n2

First equation above givesa2=12a0 Second equation in (4) givesa3=a16 And the third equation in (4) gives the recursive equation which allows us to find all an after these(n+2)(n+1)an+2+(n1)an1+an=0n2 Or(5)an+2=an(n1)an1(n+2)(n+1)n2 Therefore, for n=2 the above givesa4=a2a1(2+2)(2+1)=a2a112 But a2=12a0, therefore the above becomes(6)a4=(12a0)a1(4)(3)=12a0a112=a02a124 And for n=3 (5) givesa5=a3(31)a2(3+2)(3+1)=a32a220 But a2=12a0,a3=a1(2)(3), the above becomes(7)a5=(a1(2)(3))2(12a0)20=a1(2)(3)+a020=a1+6a0120 And for n=4 (5) givesa6=a4(41)a3(4+2)(4+1)=a43a330 But a4=a02a124 and a3=a16. The above becomesa6=(a02a124)3(a16)30=a0+2a124+a1230=a0+2a1+12a1(30)(24)=a0+14a1(30)(24)=a0+14a1720 And so on. Therefore, from (1)y(x)=n=0anxn=a0+a1x+a2x2+a3x3+a4x4+a5x5+a6x6+=a0+a1x12a0x2a16x3+(a02a124)x4+(a1+6a0120)x5+(a0+14a1720)x6+=a0+a1x12a0x2a16x3+(a024a112)x4+(a1120+a020)x5+(a0720+7a1360)x6+

Therefore(8)y(x)=a0(112x2+124x4+120x51720x6+)+a1(x16x3112x4+1120x5+7360x6+) The series solution above contains two unknowns a0,a1. There are the same as the constant of integrations. Since this is a second order ODE, then there will be two unknowns in the general solutions. These can be found from initial conditions. For example, assuming y(0)=y0,y(0)=y0. Then from (8) at x=0, it gives y(0)=a. Taking one derivative of (8) gives(9)y(x)=a0(x+424x3+)+a1(136x2+) At x=0 the above becomes y0=a1. Therefore (8) can be written as(10)y(x)=y(0)(112x2+124x4+120x51720x6+)+y(0)(x16x3112x4+1120x5+7360x6+)

4.2.10 Problem 9

   4.2.10.1 Part a
   4.2.10.2 Part b

a) Determine all the equilibrium points of the given system. b) Select two equilibrium points and classify them as saddle, node, spiral or center and whether they are stable or unstable.x=2xx2xyy=3y3xy2y2

solution

4.2.10.1 Part a

equilibrium points are the solutions in x,y of (1)2xx2xy=0(2)3y3xy2y2=0

Which can be written as(1)x(2xy)=0(2)y(33x2y)=0

From (1), we see that (3)x=0 is a solution and 2xy=0 or (4)x=2y Is another solution.  For each x value in (3,4), now we solve for y from (2). When x=0 then (2) becomesy(32y)=0 Which has solutions y=0,y=32. Therefore {0,0} and {0,32} are two solutions found so far. And when x=2y then (2) becomesy(33(2y)2y)=0y(36+3y2y)=0y(36+y)=0

Which has solutions y=0 and y=3. When y=0 then x=2y gives x=2. Therefore {2,0} is a solution, and when y=3 then x=2y gives x=23=1. Hence {1,3} is another solution. Putting all these together gives the solutions as{0,0},{0,32},{2,0},{1,3}

4.2.10.2 Part b

The given system is matrix form is[xy]=F=[2xx2xy3y3xy2y2]=[f(x,y)g(x,y)]

The Jacobian matrix for the system is given by the gradient of FJ=FJ=[fxfygxgy]=[(2xx2xy)x(2xx2xy)y(3y3xy2y2)x(3y3xy2y2)y]=[22xyx3y33x4y]

Selecting points {0,0} and {0,32} for analysis.

At Point {0,0} the linearized system matrix A is the Jacobian matrix evaluated at this equilibrium point. Hence A=[22xyx3y33x4y]x=0y=0=[2003]

The eigenvalues are found by solving det(AλI)=0 or|2λ003λ|=0(2λ)(3λ)=0

Hence λ1=2 and λ2=3. Since both eigenvalues are positive, then this is unstable critical point. It is a negative attractor also called a node.

At Point {0,32} the linearized system matrix A is the Jacobian matrix evaluated at this equilibrium point. Hence A=[22xyx3y33x4y]x=0y=32=[23203(32)34(32)]=[120923]

The eigenvalues are found by solving det(AλI)=0 or|12λ0923λ|=0(12λ)(3λ)=0

Hence λ1=12 and λ2=3. Since one eigenvalue is positive, and one eigenvalue is negative, then this is unstable critical point. It is a saddle point.