Skip to main content
1 of 2
Akku14
  • 17.4k
  • 15
  • 32

This differential equation is quite simple. So you can solve it analytically.

y[x_, \[Alpha]_] = 
     FullSimplify[
         y[x] /. First@
      DSolve[{\[Alpha]*y'''[x] + 8 y''[x] + 17 y'[x] + 10 y[x] == 0, 
        y[0] == 6*\[Alpha], y'[0] == -20, y''[0] == 84}, y, x], 
            0 >= x <= 1 && 1 <= \[Alpha] <= 2]

The is result not shown here.

And plot the DensityPlot

DensityPlot[y[x, \[Alpha]], {x, 0, 1}, {\[Alpha], 1, 2}]
Akku14
  • 17.4k
  • 15
  • 32