Scribd
Upload
52 views

Math Lab 5 - Jupyter Notebook

This document contains code from a Jupyter Notebook that is working through examples involving symbolic differentiation and solving differential equations symbolically. It imports sympy and defines symbols for variables. It then finds the general solution to a differential equation, computes the derivative of a function, finds the points where that derivative is zero to find maxima and minima, and computes the derivative of a complicated logarithmic function.

Uploaded by

redstonejazz5
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
52 views

Math Lab 5 - Jupyter Notebook

This document contains code from a Jupyter Notebook that is working through examples involving symbolic differentiation and solving differential equations symbolically. It imports sympy and defines symbols for variables. It then finds the general solution to a differential equation, computes the derivative of a function, finds the points where that derivative is zero to find maxima and minima, and computes the derivative of a complicated logarithmic function.

Uploaded by

redstonejazz5
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 1

3/4/24, 7:10 PM math lab 5 - Jupyter Notebook

In [2]: import sympy as sp


from sympy.plotting import (plot,plot_parametric)

In [3]: k, t = sp.symbols('k t')


y = sp.cos(k*t)
​
diff_eq = sp.diff(y, t, 2) + y
k_solution = sp.solve(diff_eq, k)
print(k_solution)

[-1, 1, pi/(2*t), 3*pi/(2*t)]

In [4]: A, B = sp.symbols('A B')


y_family = A*sp.sin(k*t) + B*sp.cos(k*t)
​
#substitute values of k solution to family
#check if the values of y satisfy the diff equation

In [5]: #compute dy/dx


x, y = sp.symbols('x y')
lhs = x**2-x*y+y**2
rhs = 3
dydx = sp.idiff(lhs-rhs, y, x)
dydx

Out[5]: 2𝑥 − 𝑦
𝑥 − 2𝑦
In [6]: #make eq1 and eq2, to solve it would be sp.solve then [eq1, eq2], [x,y]
horizontal_p = sp.solve([2*x-y, lhs-rhs],[x,y])
print(horizontal_p)
vertical_p = sp.solve([],[x,y])

[(-1, -2), (1, 2)]

In [7]: pcurve = sp.plot_implicit(lhs - rhs, (x,-2.5,2.5),(y,-2.5,2.5), show=False)

In [20]: #3a
x = sp.symbols('x')
y = (x**4) * ((7+x)**(1-(x**2)))
f = sp.expand_log(sp.log(y), force=True)
print(f)

(1 - x**2)*log(x + 7) + 4*log(x)

In [24]: #3b
logydiff=sp.diff((1-(x**2))*((sp.log(x+7))+4)*(sp.log(x)),x)
print(logydiff)

-2*x*(log(x + 7) + 4)*log(x) + (1 - x**2)*log(x)/(x + 7) + (1 - x**2)*(log(x + 7) + 4)/x

In [25]: #3c
ydiff = logydiff * y
print(ydiff)

x**4*(x + 7)**(1 - x**2)*(-2*x*(log(x + 7) + 4)*log(x) + (1 - x**2)*log(x)/(x + 7) + (1 - x**2)*(log(x + 7) + 4)/x)

In [28]:

ive of the derivative of the log is x**4*(x + 7)**(1 - x**2)*(-2*x*log(x + 7) + (1 - x**2)/(x + 7)) + 4*x**3*(x + 7)**(1 - x**2),

The derivative of y found with the derivative of the derivative of the log is x**4*(x + 7)**(1 - x**2)*(-2*x*log(x + 7) + (1 -
x**2)/(x + 7)) + 4*x**3*(x + 7)**(1 - x**2), the derivative of y is x**4*(x + 7)**(1 - x**2)*(-2*x*log(x + 7) + (1 - x**2)/(x +
7)) + 4*x**3*(x + 7)**(1 - x**2)

In [ ]: ​

localhost:8888/notebooks/math lab 5.ipynb 1/1

You might also like