Consider the boundary-value problem y ″ = y 3 + x , y ( a ) = α , y ( b ) = β , a ≤ x ≤ b To use the shooting method to solve this problem, one needs a starting guess for the initial slope y ′ ( a ). One way to obtain such a starting guess for the initial slope is, in effect, to do a ”preliminary shooting” in which we take a single step of Euler’s method with h = b − a.(a) Using this approach, write out the resulting algebraic equation for the initial slope.(b) What starting value for the initial slope results from this approach?I'm really not sure how to begin here; my idea was to write the 2nd order ODE using a new variable, such as: y ′ = u, u ′ = y ″ = y 3 + x. Then, maybe I could use Euler's method to solve for u = y ′ ; however, I'm not sure how to go about doing this. Any help appreciated.

Question

Consider the boundary-value problem                              y          ″                =                  y                      3                          +        x        ,                y        (        a        )        =        α        ,                y        (        b        )        =        β        ,                a        ≤        x        ≤        b            To use the shooting method to solve this problem, one needs a starting guess for the initial slope   y      ′    (  a  ). One way to obtain such a starting guess for the initial slope is, in effect, to do a ”preliminary shooting” in which we take a single step of Euler’s method with   h  =  b  −  a.(a) Using this approach, write out the resulting algebraic equation for the initial slope.(b) What starting value for the initial slope results from this approach?I&#039;m really not sure how to begin here; my idea was to write the 2nd order ODE using a new variable, such as:       y    ′    =  u,       u    ′    =      y    ″    =      y    3    +  x. Then, maybe I could use Euler&#039;s method to solve for   u  =      y    ′  ; however, I&#039;m not sure how to go about doing this. Any help appreciated.

Abbigail Vaughn · Accepted Answer

Yes, that is correct. Now the Euler step says that                              y          1                                    =                  y          0                +        h                  u          0                                              u          1                                    =                  u          0                +        h        (                  y          0          3                +                  x          0                )            As you can see, the differential equation does not really enter the equation for the value of       y    1  , so that one simply gets      u    0    =                    β            −              α                    b      −      a        ,the slope of the line connecting   (  a  ,      α    ) to   (  b  ,      β    ).Now, to get an initial guess for a multiple shooting or simultaneous collocation method, you could integrate from   (  a  ,      α    ,      u    0    ) forward and from   (  b  ,      β    ,      u    0    ) backwards with an ODE solver and blend the two solutions via cross-fading,   (  1  −  s  )      y    0    (  a  +  s  h  )  +  s      y    1    (  a  +  s  h  ). (Be careful in its use, this formula is or at least was patented).

Answered question

Answer & Explanation

New Questions in Integral Calculus