Scribd
Upload
46 views4 pages

Lab

Uploaded by

nayanamirajkar14
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
46 views4 pages

Lab

Uploaded by

nayanamirajkar14
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 4

B ) Finding the angle between the radius vector and the tangent , angle

between two curves.


# Angle between the radius vector and the tangent for the curve
# r=a*(1+cost) at t=pi/3
from sympy import *
a,t=symbols('a,t')
R=a*(1+cos(t))
dRdt=diff(R,t)
R=R.subs(t,pi/3)
dRdt=dRdt.subs(t,pi/3)
PHI=atan(R/dRdt)
if PHI<0:
PHI=PHI+pi
print('The angle between the radius vector and the tangent =',PHI)

C ) Find the angle between the curves r =a(1−cos t) and r =2 a cos t at


1
t=acos ( )
3

from sympy import *


a,t=symbols('a,t')
R=a*(1-cos(t))
dRdt=diff(R,t)
R=R.subs(t,acos(1/3))
dRdt=dRdt.subs(t,acos(1/3))
PHI=atan(R/dRdt)
if PHI<0:
PHI=PHI+pi
r=2*a*cos(t)
drdt=diff(r,t)
r=r.subs(t,acos(1/3))
drdt=drdt.subs(t,acos(1/3))
phi=atan(r/drdt)
if phi <0:
phi=phi+pi
print('The angle of intersection =',abs(PHI-phi))

3. Finding partial derivatives, Jacobian and plotting the graph


−1 y
A) Findu x , u yx , u xx , u xy for u=tan ( x )

from sympy import *


import matplotlib.pyplot as plt
var('x y')
u=atan(y/x)
ux=simplify(diff(u,x))
uy=simplify(diff(u,y))
uxy=simplify(diff(ux,y))
uxx=simplify(diff(u,x,2))
uyx=simplify(diff(u,x,y))
print('The required partial derivatives:')
print('ux=',ux)
print('uy=',uy)
print('uxy=',uxy)
print('uyx=',uyx)
print('uxx=',uxx)
2 2
y x ∂f ∂f ∂ x ∂ f
B) If f ( x , y )=x + y , find ,
∂x ∂ y
∧also show that 2 = 2
∂x ∂ y

from sympy import *


var('x y')
f=x**y+y**x
fx=simplify(diff(f,x))
fy=simplify(diff(f,y))
fxx=simplify(diff(f,x,2))
fyy=simplify(diff(f,y,2))
print('The required partial derivatives:')
print('fx=',fx)
print('fy=',fy)
print('fxx=',fxx)
print('fyy=',fyy)

C) Jacobian of a given function


∂(u , v )
If u=2 x−3 y , v =5 x + 4 y find ∂(x , y )

from sympy import *


import numpy as np
import matplotlib.pyplot as plt
var('x y')
u=2*x-3*y
v=5*x+4*y
A=Matrix([u,v])
Jmatrix=A.jacobian([x,y])
print('Jacobian matrix-')
pprint(Jmatrix)
J=det(Jmatrix)
print(' Jacobian =', J)

4. Solution of first order differential equation and plotting the graphs


dy 2 3
A) Solve dx − x y =x with y(1)=0 and verify using Python code.

from sympy import *


import matplotlib.pyplot as plt
x=symbols('x')
y=Function('y')
deq=Derivative(y(x),x)-2*y(x)/x-x**3
ysoln=simplify(dsolve(deq,y(x),ics={y(1):0},hint='1st_linear'))
print("The solution of given differential equation is:")
pprint(simplify(ysoln))
plot(ysoln.rhs,(x,-2,2))

dy 2
B) Solve x dx + y= y log x with y(1)=1 and verify using Python code.

from sympy import *


import matplotlib.pyplot as plt
x=symbols('x')
y=Function('y')
deq=x*Derivative(y(x),x)+y(x)-log(x)*y(x)**2
ysoln=dsolve(deq,y(x),ics={y(1):1},hint='Bernoulli')
print("The solution of the given differential equation is:")
pprint(ysoln)
plot(ysoln.rhs,(x,1,5))

You might also like