Scribd
Upload
4 views3 pages

MLL703 Minor Tutorial222222

The document outlines the instructions and questions for a laboratory exam in Mathematical and Computational Methods in Materials, scheduled for 23-Sep-2024. It includes five questions requiring programming tasks in Python, such as input handling, solving a quadratic equation, plotting a decay curve, and implementing Euler's method. Students must submit their answers in a Jupyter notebook format and are warned against using AI-generated content.

Uploaded by

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

MLL703 Minor Tutorial222222

The document outlines the instructions and questions for a laboratory exam in Mathematical and Computational Methods in Materials, scheduled for 23-Sep-2024. It includes five questions requiring programming tasks in Python, such as input handling, solving a quadratic equation, plotting a decay curve, and implementing Euler's method. Students must submit their answers in a Jupyter notebook format and are warned against using AI-generated content.

Uploaded by

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

MLL703

Mathematical and Computational Methods in Materials


Laboratory exam

Time: 1.5 hours DATE: 23-Sep-2024

Instructions:

1. Attempt all questions.


2. Do not copy from ChatGPT or Gemini. Your answers will be
run through Turnitin, which can detect AI-generated content.
3. Submit all your answers in a single .ipynb file (Jupyter notebook).
Make sure each question is clearly labeled with comments and
markdown. Upload your .ipynb file on Moodle with the format
keybrosid_name.ipynb.
4. No extra time will be given.

Question 1: Write a program that takes the input of a student's name and their marks, ensure
that the marks are entered as a floating-point number. Print the student's name and their marks,
along with the type of both variables. (3 marks)

Question 2: Write a Python code to solve a quadratic equation (10 marks)


𝑎𝑥 2 + 𝑏𝑥 + 𝑐 = 0

• Take input of a, b, and c from the user using the input function.
• Use the if else statement.
o Syntax of if else is
if condition:
# block of code executed if the condition is true
elif another_condition:
# block of code executed if the previous condition is false
and this one is true
else:
# block of code executed if all the above conditions are false

Question 3: Write a print(f"{}) statement to display the values of a = 5 and b = 10 in the


format: "The values of a and b are 5 and 10." (2 marks)
Question 4: A radioactive substance decays according to a decay constant of λ=0.1 (per unit
time), with an initial quantity of N0=1000. (10 marks)

1. Formulate the ordinary differential equation (ODE) that describes this decay
process.
2. Utilize the solve_ivp function from the scipy.integrate module to numerically solve
for the remaining quantity of the substance over time.
3. Plot the decay curve over the time interval from t=0 to t=50 using
the matplotlib.pyplot module.
4. Ensure that both axes are labeled using plt.xlabel and plt.ylabel, and provide a title for
the graph.

Question 5: Fill in the blanks: (5 marks)

import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import solve_ivp
# Function representing the ODE: dy/dx = f(x, y)
def f(x, y):
return _____________ # ODE: y’+2y=x3*e-2x

# Euler's method implementation


def euler_method(f, x0, y0, h, x_end):
# Create arrays to store the results
x_values = np.arange(_________) # create values from x0 to x_end with interval of h
y_values = np.zeros(len(x_values))
# Initial condition
y_values[0] = y0
# Euler's method iterative process
for i in range(1, len(x_values)):
y_values[i] = ________________________ # formula of euler method
return x_values, y_values
# Solve the ODE using solve_ivp
def solve_with_solve_ivp(f, x0, y0, x_end):
sol = solve_ivp(f, [x0, x_end], [y0], t_eval=np.arange(x0, x_end, 0.01))
return sol.t, sol.y[0]

# Initial conditions and parameters


x0 = 0
y0 = 1 # y(x0) = y0
h = 0.01 # Step size for Euler's method
x_end = 2 # We want to approximate up to x = 2

# Run Euler's method


x_euler, y_euler = euler_method(f, x0, y0, h, x_end)

# Solve using solve_ivp


x_solve_ivp, y_solve_ivp = ________________ # solve using function ivp

# Plot both methods


plt.plot(x_euler, y_euler, 'bo--', label="Euler's Method")
plt.plot(x_solve_ivp, y_solve_ivp, 'r-', label="solve_ivp (RK45)")

plt.xlabel('x')
plt.ylabel('y')
plt.title("Comparison of Euler's Method and solve_ivp for ODE")
plt.legend()
plt.grid(True)
plt.show()

You might also like