Scribd Logo
Upload

Welcome to Scribd!

  • 5 views

    Assign 6 SEM2

    Uploaded by

    Vadanya Sehgal
    The document discusses solving differential equations using Mathematica. It presents several examples of second and third order differential equations, along with their initial conditions. The solutions are plotted to visualize the behavior of the functions.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd

    Assign 6 SEM2

    Uploaded by

    Vadanya Sehgal
    5 views2 pages
    The document discusses solving differential equations using Mathematica. It presents several examples of second and third order differential equations, along with their initial conditions. The solutions are plotted to visualize the behavior of the functions.

    Original Title

    Assign6SEM2

    Copyright

    © © All Rights Reserved

    Available Formats

    PDF, TXT or read online from Scribd

    Did you find this document useful?

    Is this content inappropriate?

    The document discusses solving differential equations using Mathematica. It presents several examples of second and third order differential equations, along with their initial conditions. The solutions are plotted to visualize the behavior of the functions.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    5 views2 pages

    Assign 6 SEM2

    Uploaded by

    Vadanya Sehgal
    The document discusses solving differential equations using Mathematica. It presents several examples of second and third order differential equations, along with their initial conditions. The solutions are plotted to visualize the behavior of the functions.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    of 2

    In[1]:= EQ1 = y '''[x] - y ''[x] + 100 * y '[x] - 100 y[x] ⩵ 0

    Out[1]= - 100 y[x] + 100 y′ [x] - y′′ [x] + y(3) [x] ⩵ 0

    In[2]:= Sol1 = DSolve[EQ1, y[x], x]


    y[x] → ⅇ 3 + 1 Cos10 x + 2 Sin10 x
    x
    Out[2]=

    In[4]:= PlotEvaluate[Sol1〚1, 1, 2〛 /. {C[1] → {1, 2, 3}, C[2] → {1, 2, 3}, C[3] → {1, 2, 3}}],
    {x, 0, 10}, PlotLegends → "Expressions",
    PlotStyle → Thickness[0.01], Thickness[0.002], Thickness[0.003]

    15 000

    ⅇx + cos10 x + sin10 x
    10 000
    Out[4]= 2 ⅇx + 2 cos10 x + 2 sin10 x

    3 ⅇx + 3 cos10 x + 3 sin10 x
    5000

    2 4 6 8 10

    In[24]:= Quit

    In[10]:= EQ2 = y '''[x] + 7.5 * y ''[x] + 14.25 y '[x] - 9.125 * y[x] ⩵ 0


    Out[10]=
    - 9.125 y[x] + 14.25 y′ [x] + 7.5 y′′ [x] + y(3) [x] ⩵ 0

    In[11]:= IC1 = y[0] ⩵ 10.05


    Out[11]=

    y[0] ⩵ 10.05

    In[12]:= IC2 = y '[0] ⩵ - 54.975

    Out[12]=

    y′ [0] ⩵ - 54.975

    In[13]:= IC3 = y ''[0] ⩵ 257.5125


    Out[13]=

    y′′ [0] ⩵ 257.513

    In[14]:= Sol2 = DSolve[{EQ2, IC1, IC2, IC3}, y[x], x]


    Out[14]=

    y[x] → 0.05 ⅇ-4. x 1. ⅇ4.5 x + 200. Cos1.5 x - 200. Sin1.5 x 


    2

    In[15]:= Plot[Evaluate[Sol2〚1, 1, 2〛], {x, - 10, 10}]


    Out[15]=

    9
    1 × 10

    8
    5 × 10

    -10 -5 5 10

    8
    -5 × 10

    In[1]:= Quit

    In[5]:= EQ1B = y '''[x] + 0.55 y ''[x] + 4.3 y '[x] ⩵ 0


    Out[5]= 4.3 y′ [x] + 0.55 y′′ [x] + y(3) [x] ⩵ 0

    In[6]:= Sol1B = DSolve[EQ1B, y[x], x]

    y[x] → 3 + ⅇ - 0.477983 1 - 0.0639535 2  Cos2.05533 x +


    -0.275 x
    Out[6]=

    - 0.0639535 1 + 0.477983 2  Sin2.05533 x

    In[8]:= PlotEvaluate[Sol1B〚1, 1, 2〛 /. {C[1] → {1, 2, 3}, C[2] → {1, 2, 3}, C[3] → {1, 2, 3}}],
    {x, 0, 10}, PlotLegends → "Expressions",
    PlotStyle → Thickness[0.001], Thickness[0.002], Thickness[0.003]

    3 1 + ⅇ-0.275 x -0.541937 cos2.05533 x +

    Out[8]= 2 + ⅇ-0.275 x -1.08387 cos2.05533 x + 0.82806

    3 + ⅇ-0.275 x -1.62581 cos2.05533 x + 1.24209


    2

    2 4 6 8 10

    You might also like