Semester 6
Complex Analysis
Practical File
NAME : SHREEJI
COURSE: MATHS HONS (3rd year)
SECTION:A
ROLL NO. :15032563057
INDEX
[Link]. Name of the Practical
1. To compute the poles and corresponding residues of complex
functions.
2. To compute the integral over a straight line path.
3. Declaring complex numbers and graphical representation.
4. Algebra of complex numbers.
5. To find modulus, conjugate and face angle of an array of
complex numbers.
6. Contour Integration.
7. Plot the complex function and analyze the graph.
8. Taylor’s Series expansion of a given function at a given point.
9. To perform Laurent Series expansion.
Q. 1 To compute the poles and corresponding residues of
complex functions.
1. 4Z+1/Z2-Z-2
>>num=[4 1]
num =
4 1
>>den=[1 -1 -2]
den =
1 -1 -2
>> [r,p,k]=residue(num,den)
r=
p=
-1
k=
[]
2. 4Z+3/2Z3-3.4Z2+1.98Z-0.408
>>num=[4 3]
num =
4 3
>>den=[2 -3.4 1.98 -0.408]
den =
2.0000 -3.4000 1.9800 -0.4080
>> [r,p,k]=residue(num,den)
r=
32.5402
-16.2701 +12.3344i
-16.2701 -12.3344i
p=
0.7118
0.4941 + 0.2061i
0.4941 - 0.2061i
k=
[]
3. 2Z+1/Z3+5Z2+8Z+4
>>num=[2 1]
num =
2 1
>>den=[1 5 8 4]
den =
1 5 8 4
>> [r,p,k]=residue(num,den)
r=
1.0000
3.0000
-1.0000
p=
-2.0000
-2.0000
-1.0000
k=
[]
Q.2. To compute the integral over a straight line path
%AIM- TO PERFORM LINE INTEGRATION
% Compute integral over a straight line path between two
specified points
%f=z+1, where C is the straight line path from -1+i to 2-i
function line_integration
syms t;
x=t;
y=(1-2*t)/3;
z=x+i*y;
f=z+1;
integrand=f*diff(z,t);
F=int(integrand,'t',-1,2)
Q3) Declaring complex numbers and graphical representation
Function plot_complex_no
z1=3+2*i;
z2=7-4*i;
plot(z1,'*');
hold on;
plot(z2,'*');
hold off;
end
Q4) Algebra of complex numbers
Function algebra_complexno
z1=3+2*i;
z2=7-4*i;
z3=z1+z2;
z4=z1-z2;
z5=z1*z2;
z6=z1/z2;
str1=['Sum of the complex numbers is:' num2str(z3)];
str2=['Difference of the complex numbers is:' num2str(z4)];
str3=['Product of the complex numbers:' num2str(z5)];
str4=['Division in the complex numbers gives:' num2str(z6)];
disp(str1);
disp(str2);
disp(str3);
disp(str4);
end
>>algebra_complexno
Sum of the complex numbers is:10-2i
Difference of the complex numbers is:-4+6i
Product of the complex numbers:29+2i
Division in the complex numbers gives:0.2+0.4i
>>
Q5) To find modulus, conjugate and face angle of an array of
complex numbers
Function complex_array
z=input('enter the complex numbers: ');
absolute = abs(z);
conjugate = conj(z);
theta = angle(z);
str1 = ['the absolute value of the complex number is: ' num2str(absolute)];
str2 = ['the conjugate of the complex number is: ' num2str(conjugate)];
str3 = ['the phase angle of the complex number is: ' num2str(theta)];
disp(str1);
disp(str2);
disp(str3);
end
output
>>complex_array
enter the complex numbers: [1+2i, 3+5i, -2-7i]
the absolute value of the complex number is: 2.2361 5.831 7.2801
the conjugate of the complex number is: 1-2i 3-5i -2+7i
the phase angle of the complex number is: 1.1071 1.0304 -1.8491
Q.6 Contour Integration
function contour_integration1
symst;
y=t;
x=t^2+1;
z=x+i*y;
f=z.^2-2.*z+1;
integrand=f*diff(z,t);
answer=int(integrand,'t',-2,2)
>> contour_integration1
answer =
176/3*i
function contour_integration2
symst;
y=sin(t);
x=cos(t);
z=x+i*y;
f=z.^3+2.*z.^2+1;
integrand=f*diff(z,t);
answer=int(integrand,'t',0,2*pi)
>> contour_integration2
answer =
function contour_integration3
symst;
x=t;
y=(-2*t+1)/3;
z=x+i*y;
f=z+1;
integrand=f*diff(z,t);
answer=int(integrand,'t',1,-1)
>> contour_integration3
answer =
-22/9+2/3*i
Q. 7 Plot the complex function and analyze the graph
Function complexgrid
z=cplxgrid(40);
cplxmap(z,real(z))
figure
cplxmap(z,imag(z))
1
1
0.5
0.5
0
0
-0.5
-0.5
-1
-1
1 1
0.5 1 0.5 1
0 0.5 0.5
0
0 0
-0.5 -0.5 -0.5 -0.5
-1 -1 -1 -1
function complexgrid2
z=cplxgrid(40);
cplxmap(z,real(z.^2))
figure
cplxmap(z,imag(z.^2))
1
0.5
0.5
0
0
-0.5
-0.5
-1
-1 1
1
0.5 1
0.5 1
0 0.5
0 0.5
0 0
-0.5 -0.5 -0.5
-0.5
-1 -1 -1 -1
function complexgrid3
z=cplxgrid(40);
cplxmap(z,real(z.^3))
figure
cplxmap(z,imag(z.^3))
1
0.5
0.5
0
0
-0.5
-0.5
-1
-1 1
1
0.5 1
0.5 1 0.5
0.5 0
0 0
0 -0.5
-0.5 -0.5
-0.5
-1 -1 -1 -1
function complexgrid4
z=cplxgrid(40);
cplxmap(z,real((z.^4-1).^(1/4)))
figure
cplxmap(z,imag((z.^4-1).^(1/4)))
1
0.8
0.6 0.5
0.4 0
0.2 -0.5
0
1 -1
1
0.5 1
0.5 1
0 0.5
0 0.5
0 0
-0.5 -0.5 -0.5 -0.5
-1 -1 -1 -1
function complexgrid5
z=cplxgrid(40);
cplxmap(z,real(conj(z)))
figure
cplxmap(z,imag(conj(z)))
1
0.5
0.5
0
0
-0.5
-0.5
-1
1 -1
1
0.5 1
0.5 0.5 1
0 0.5
0 0
-0.5 0
-0.5 -0.5 -0.5
-1 -1 -1 -1
function complexgrid6
z=cplxgrid(40);
cplxmap(z,real(exp(z)))
figure
cplxmap(z,imag(exp(z)))
1.5
3
1
2.5
0.5
2
0
1.5
-0.5
1
-1
0.5
0 -1.5
1 1
0.5 1 0.5 1
0.5 0 0.5
0
0 0
-0.5 -0.5 -0.5
-0.5
-1 -1 -1 -1
function complexgrid7
z=cplxgrid(40);
cplxmap(z,real(z+1))
figure
cplxmap(z,imag(z+1))
1
2
0.5
1.5
0
1
-0.5
0.5
0 -1
1 1
0.5 1 0.5 1
0.5 0 0.5
0
0 0
-0.5 -0.5 -0.5
-0.5
-1 -1 -1 -1
Q.8. Taylor’s Series expansion of a given function at a given
point
function complexgrid6
z=cplxgrid(40);
cplxmap(z,real(exp(z)))
figure
cplxmap(z,imag(exp(z)))
1.5
3
1
2.5
0.5
2
0
1.5
-0.5
1
-1
0.5
0 -1.5
1 1
0.5 1 0.5 1
0.5 0 0.5
0
0 0
-0.5 -0.5 -0.5
-0.5
-1 -1 -1 -1
function complexgrid7
z=cplxgrid(40);
cplxmap(z,real(z+1))
figure
cplxmap(z,imag(z+1))
1
2
0.5
1.5
0
1
-0.5
0.5
0 -1
1 1
0.5 1 0.5 1
0.5 0 0.5
0
0 0
-0.5 -0.5 -0.5
-0.5
-1 -1 -1 -1
Taylor series expansion of a given function at a given point
function taylor1
syms xy;
z=x+i*y
f=exp(z)
taylor(f,6,0)
z=
x+i*y
f=
exp(x+i*y)
ans =
exp(i*y)+exp(i*y)*x+1/2*exp(i*y)*x^2+1/6*exp(i*y)*x^3+1/24*exp(i*y)*x^4+1/120*exp(i*y)*x^5
function taylor2
symsxy;
z=x+i*y
f=exp(z^2)
taylor(f,6,0)
z =
x+i*y
f =
exp((x+i*y)^2)
ans =
exp(-y^2)+2*i*exp(-y^2)*y*x+exp(-y^2)*(1-2*y^2)*x^2+exp(-
y^2)*(2*i*y-4/3*i*y^3)*x^3+exp(-y^2)*(1/2-
2*y^2+2/3*y^4)*x^4+exp(-y^2)*(i*y-4/3*i*y^3+4/15*i*y^5)*x^5
Q. 9 To perform Laurent Series expansion
%AIM:to perform laurents series expansion of a given function f(z) around a
%given point
functionlaurent_series
symsfz;
f=(sin(z)-1)/(z^4);
disp('the laurent series expansion of the function at z=0')
laurent_s=maple('series',f,'z=0',8)
end
laurent_series
thelaurent series expansion of the function at z=0
laurent_s =
series(-1*z^(-4)+1*z^(-3)-1/6*z^(-1)+1/120*z-1/5040*z^3+O(z^4),z,4)
function laurent_series1
symsfz;
f=cot(z)/(z^4);
disp('the laurent series expansion of the function at z=0')
laurent_s=maple('series',f,'z=0',8)
end
laurent_series1
thelaurent series expansion of the function at z=0
laurent_s =
series(1*z^(-5)-1/3*z^(-3)-1/45*z^(-1)-2/945*z-1/4725*z^3+O(z^4),z,4)
function laurent_series2
symsfz;
f=1/z*(z+5);
disp('the laurent series expansion of the function at z=0')
laurent_s=maple('series',f,'z=0',8)
end
laurent_series2
thelaurent series expansion of the function at z=0
laurent_s =
series(5*z^(-1)+1,z)
