Scribd
Upload
18 views

MAT120 Lab Assignment

The document provides instructions for a mathematics assignment involving calculus concepts such as derivatives, integrals, and special functions. It outlines 6 problems to solve related to derivatives and plots, Simpson's rule, Fresnel integrals, optimization of functions, the beta function, and the gamma function. Students are instructed to show work and properly display solutions, and plagiarism will result in zero points. The assignment is due by November 25th, 2023.

Uploaded by

shudipto.siam
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
18 views

MAT120 Lab Assignment

The document provides instructions for a mathematics assignment involving calculus concepts such as derivatives, integrals, and special functions. It outlines 6 problems to solve related to derivatives and plots, Simpson's rule, Fresnel integrals, optimization of functions, the beta function, and the gamma function. Students are instructed to show work and properly display solutions, and plagiarism will result in zero points. The assignment is due by November 25th, 2023.

Uploaded by

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

Department of Mathematics and Natural Sciences

MAT120 Lab - Integral Calculus and Differential Equations


Assignment

Copying of codes from any websites or other classmates will not be tolerated.
You are ONLY ALLOWED to take help from the colab files that I have provided
but be careful NOT to copy any codes.
The outputs must be organized and readable. It must show necessary steps of
the calculation and properly display the solutions/results.
Be creative, use your intuition. Answer the questions by yourself and Do not
Cheat. Plagiarized copies will receive zero points regardless of the circum-
stances.
Total marks is 40
Deadline: 23:59 25th November, 2023

Attempt All Questions

1. Derivatives and Plots: (5+5 = 10)


For the following functions:

1. f (x) = 2x(x − 3)(sin x − 5), −2π ≤ x ≤ 2π


1 2
−2t
e√
2. h(t) = 2π
, −4 ≤ t ≤ 4

a Plot these functions and their derivatives. Include enough points so that the
curve you plot appears smooth. Use different colors and label to represent the
function and its derivative respectively. Use gray lines, one horizontal at y = 0
and the other vertical at x = 0. Limit the values of y so that the plot is
understandable. Label the axes x and y or t and y respectively.
b Use Trapezoidal rule and quad to integrate each of these functions in the given
domain. Print the results such that it shows the integral and the respective
solutions. For example,
Z 4 − 1 t2
e 2
√ dt
−4 2π
Result from Trapezoidal rule:
Result from quad:
2. Simpson’s Rule: (5+2+3 = 10)
Simpson’s 1/3 Rule is a numerical method used to approximate definite integrals.
It’s based on approximating the area under a curve using quadratic approximations.
The formula is given as
Z b
h 
f (x) dx = f (a) + 4f (x1 ) + 2f (x2 ) + ... + 4f (xn−1 ) + 2f (xn−2 ) + f (b) (1)
a 3
The algorithm:

• Inputs: A function f (x), specified interval [a, b] and a value for the number of
sub-intervals N.
• Compute the width of each subinterval h : h = (b−a)
N
• Set new variables ’x i = a’ and ’sum = f(a) + f(b)’.
• Begin a loop that runs from i = 1 till i = N − 1. In each iterations, update the
value of x : x i += i*h check if i is even or odd.
• If odd then sum += 4*f(x i) otherwise if even then sum += 2*f(x i)
• Finally update sum = h/3 * sum and return the output.

For the following integrals:


Z 4  
2
I1 = sin(x) − dx
1 x
Z 3
I2 = (x2 + 2x + 1) dx
1

1 rd
a) Implement the Simpson’s 3
rule to evaluate both of the integrals.
b) Evaluate both integrals using the package scipy.integrate and show both result
and the error separately as outputs keeping atleast six decimal places.
c) Use sympy to solve both the integrals symbolically without the limits. Plot
the resulting functions in a single graph and label each function correctly using
legends to show it on the plot.

3. Fresnel Integral: (5)


If one looks at the diffraction of a plane wave, the diffracted wave (in the Fresnel
approximation) is given by the expression
Z z
exp ix2 = C(z) + iS(z)

Ψ(z) = (2)
0

where
Z z Z z
2
sin x2 dx
 
C(z) = cos x dx S(z) =
0 0
Carry out these integrals numercially (Using any preferred algorithm e.g. Trapezoidal
Rule or Monte Carlo) for values of z ∈ [0, 20] and generate a parametric plot S vs C.
The graph plotted is known as the Cornu Spiral.

4. Highschool Calculus: (5 + 5 = 10)

a) Plot the function


x3 − 100 cos(x)
in the interval [-1,4] with axis labelling and grid. Find the root of this function
in the same interval using Newton’s Method or Bisection Method.
b) Find the maximum value of

x2
f (x) = exp(x) sin(x) − +5
2
for −1 ≤ x ≤ 3.
(EVALUATE THE FIRST AND SECOND DERIVATIVES TO FIND THE
DESIRED EXTREMA)

Special Functions:

5. The Beta function β(x, y) is defined for all x, y > 0 by


Z π
2
β(x, y) = 2 (sin t)2x−1 (cos t)2y−1 dt
0

Write a function named Beta(x, y) that takes x and y as arguments and returns the
result of the integral.
(a) Evaluate β 5, 4). (2)

6. The gamma function Γ(n) is defined for all n > 0 by


Z ∞
Γ(n) = xn−1 e−x dx
0

Write a function named Gamma(n) that takes n as an argument and returns the
result of the integral. (Use mathinf by using import math)
 √
(a) Show that Γ 21 = π. (2)
(b) Find Γ 72 . Provide both rational and decimal output.

(1)

You might also like