Let’s look at a calculus problem involving differentiation: finding the rate of change of a function at any given point.

For this example, we’ll consider a simple polynomial function and calculate its derivative.

Then, we’ll use $Python$ to graph the function along with its derivative, making it easy to understand how the slope changes across different points.

Problem: Differentiating a Polynomial Function

Given a polynomial function:
$$
f(x) = 3x^3 - 5x^2 + 2x - 4
$$
we want to:

1. Find its derivative, $ f’(x) $.
2. Plot both $ f(x) $ and $ f’(x) $ over a specified range to observe the behavior of the function and how the slope changes with $ x $.

Solution

1. Differentiate $ f(x) $ with respect to $ x $.
   • The derivative, $ f’(x) $, represents the rate of change or slope of $ f(x) $ at each point.
2. Use $Python$ to plot both $ f(x) $ and $ f’(x) $ on the same graph.

Derivative Calculation

The derivative of $ f(x) = 3x^3 - 5x^2 + 2x - 4 $ is:
$$
f’(x) = 9x^2 - 10x + 2
$$

Python Code

Here’s the $Python$ code to calculate and plot both $ f(x) $ and $ f’(x) $:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29

import numpy as np
import matplotlib.pyplot as plt

# Define the function and its derivative
def f(x):
    return 3 * x**3 - 5 * x**2 + 2 * x - 4

def f_prime(x):
    return 9 * x**2 - 10 * x + 2

# Generate x values for the plot
x_values = np.linspace(-2, 3, 100)
y_values = f(x_values)
y_prime_values = f_prime(x_values)

# Plot f(x) and f'(x)
plt.figure(figsize=(10, 6))
plt.plot(x_values, y_values, label="f(x) = 3x^3 - 5x^2 + 2x - 4", color="blue")
plt.plot(x_values, y_prime_values, label="f'(x) = 9x^2 - 10x + 2", color="red", linestyle="--")

# Add labels and title
plt.xlabel("x")
plt.ylabel("y")
plt.title("Function and its Derivative")
plt.legend()
plt.grid(True)
plt.axhline(0, color="black",linewidth=0.5)
plt.axvline(0, color="black",linewidth=0.5)
plt.show()

Explanation of the Code

1. Function Definitions:
   • f(x) is defined to represent $ 3x^3 - 5x^2 + 2x - 4 $.
   • f_prime(x) represents the derivative $ 9x^2 - 10x + 2 $.
2. Plotting:
   • We use a range of $ x $ values from $-2$ to $3$ to cover the region of interest.
   • Both $ f(x) $ and $ f’(x) $ are plotted on the same graph for easy comparison.
   • The function $ f(x) $ is shown in blue, and the derivative $ f’(x) $ is shown as a red dashed line.

Visualization

The plot shows:

• Blue Curve: Represents $ f(x) $, the original polynomial function.
• Red Dashed Curve: Represents $ f’(x) $, the derivative, showing the slope of $ f(x) $ at each point.

The intersection of $ f’(x) $ with the $x$-$axis$ indicates where $ f(x) $ has a local maximum or minimum (i.e., where the slope of $ f(x) $ is zero).

Interpretation

1. Function Behavior: $ f(x) $ varies according to the cubic shape, with local peaks and troughs.
2. Derivative Insights: $ f’(x) $ reveals where $ f(x) $ is increasing or decreasing:
   • When $ f’(x) > 0 $: $ f(x) $ is increasing.
   • When $ f’(x) < 0 $: $ f(x) $ is decreasing.
3. Applications: Understanding the derivative function is crucial in finding the rate of change, optimizing functions, and analyzing trends, especially in fields like $physics$, $economics$, and $engineering$.

Let’s tackle a basic econometrics problem often used to understand relationships between variables in economics: Simple Linear Regression.

This example can be applied to analyze the relationship between two variables, such as income and expenditure, or education and wage.

Problem: Relationship Between Education Level and Wage

The goal is to analyze the relationship between years of education and wage.

This example assumes we have a dataset where each data point represents an individual’s years of education and corresponding wage.

By fitting a simple linear regression model, we can determine if there is a statistically significant relationship between these variables and interpret the results.

Objective

Use Python to:

1. Fit a simple linear regression model where the independent variable $(X)$ is Years of Education and the dependent variable $(Y)$ is Wage.
2. Plot the regression line with the data points to visualize the relationship.

Example Dataset

For simplicity, we’ll create a synthetic dataset of years of education and wage, simulating a realistic scenario where wage generally increases with more years of education, though with some variability.

Python Code

Here’s the $Python$ code to create the dataset, fit the linear regression model, and visualize the results:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41

import numpy as np
import matplotlib.pyplot as plt
from sklearn.linear_model import LinearRegression
from sklearn.metrics import r2_score

# Generate synthetic data
np.random.seed(0)  # For reproducibility
years_of_education = np.random.randint(10, 20, 50)  # Education years from 10 to 20
wage = 2 * years_of_education + np.random.normal(0, 5, 50)  # Linear relation with noise

# Reshape the data for sklearn
X = years_of_education.reshape(-1, 1)
y = wage

# Fit the Linear Regression model
model = LinearRegression()
model.fit(X, y)

# Predict wages based on the model
predicted_wage = model.predict(X)

# Calculate R-squared value for model performance
r_squared = r2_score(y, predicted_wage)
print(f"R-squared: {r_squared:.2f}")

# Plot the data and regression line
plt.figure(figsize=(10, 6))
plt.scatter(years_of_education, wage, color="blue", label="Actual Data", alpha=0.6)
plt.plot(years_of_education, predicted_wage, color="red", label="Regression Line", linewidth=2)
plt.xlabel("Years of Education")
plt.ylabel("Wage")
plt.title("Relationship Between Education Level and Wage")
plt.legend()
plt.grid(True)
plt.show()

# Output model parameters
slope = model.coef_[0]
intercept = model.intercept_
print(f"Slope (Effect of Education on Wage): {slope:.2f}")
print(f"Intercept: {intercept:.2f}")

Explanation of the Code

1. Data Generation:
   • We generate random years of education for $50$ individuals, ranging from $10$ to $20$ years.
   • Wage is simulated as a linear function of education with added noise to introduce variability.
2. Model Fitting:
   • The LinearRegression model is used to fit the data, with Years of Education as the predictor and Wage as the outcome.
3. Visualization:
   • A scatter plot displays the actual data points.
   • The regression line represents the predicted wage based on years of education.
4. Model Evaluation:
   • The R-squared value measures how well the model explains the variance in wage based on years of education.
   • The slope of the regression line indicates the increase in wage for each additional year of education.

Visualization

The plot shows:

• Data Points: Each point represents an individual’s education level and wage.
• Regression Line: The red line shows the model’s predictions, helping visualize the general trend that wage increases with education level.

Interpretation

1. Slope: If the slope is, for example, $2.5$, this suggests that each additional year of education is associated with an increase of ＄$2.50$ in wage.
2. Intercept: This is the estimated wage for an individual with $0$ years of education.
3. R-squared: A higher R-squared value (close to $1$) indicates that the model effectively explains the variance in wage, though it’s normal to have a moderate value in real-world data with noise.

Applications in Econometrics

Simple linear regression is a foundational tool in econometrics used to explore and quantify relationships between economic variables.

By understanding such relationships, economists can make informed predictions and policy recommendations.

Let’s explore cryptography problem: Modular Exponentiation.

This operation is fundamental in many cryptographic algorithms, especially public-key cryptography methods like $RSA$ encryption.

Modular exponentiation efficiently computes large powers of a number under a modulus, which is essential for encryption and decryption processes.

Problem: Modular Exponentiation

Given three integers $( a )$, $( b )$, and $( n )$, calculate $( a^b \mod n )$.

For large $( b )$, directly computing $( a^b )$ and then taking the modulus $( n )$ would be computationally intensive.

Instead, we can use an efficient algorithm called Exponentiation by Squaring to handle large powers.

Example

Let’s calculate $( 3^{45} \mod 13 )$:

• $( a = 3 )$
• $( b = 45 )$
• $( n = 13 )$

The solution involves repeatedly squaring and taking the modulus to avoid computing the large number $( 3^{45} )$ directly.

Objective

Use $Python$ to implement modular exponentiation and calculate $( a^b \mod n )$ efficiently.

Then, visualize the results of $( a^b \mod n )$ for various values of $( b )$ to see how the result changes as $( b )$ increases.

Python Code

Here’s the $Python$ code to solve the modular exponentiation problem and visualize the results:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33

import matplotlib.pyplot as plt

# Modular exponentiation function using Exponentiation by Squaring
def modular_exponentiation(a, b, n):
    result = 1
    base = a % n
    while b > 0:
        if b % 2 == 1:  # If b is odd, multiply the result by the base
            result = (result * base) % n
        b = b // 2
        base = (base * base) % n  # Square the base
    return result

# Parameters
a = 3  # Base
n = 13  # Modulus

# Calculate modular exponentiation for different powers of b
b_values = list(range(1, 46))  # Powers from 1 to 45
mod_exp_results = [modular_exponentiation(a, b, n) for b in b_values]

# Plot the results
plt.figure(figsize=(10, 6))
plt.plot(b_values, mod_exp_results, marker='o', linestyle='-', color='teal')
plt.xlabel('Exponent (b)')
plt.ylabel(f'{a}^b mod {n}')
plt.title(f'Modular Exponentiation Results of {a}^b mod {n}')
plt.grid(True)
plt.show()

# Test output for a specific case
specific_result = modular_exponentiation(a, 45, n)
print(f"3^45 mod 13 = {specific_result}")

Explanation of the Code

1. Modular Exponentiation Calculation:
   • The modular_exponentiation function implements the Exponentiation by Squaring algorithm.
     This method allows efficient calculation of large powers under a modulus without computing the entire number.
   • Starting with result = 1 and reducing the exponent by half in each step, we only multiply and take modulus operations.
2. Plotting Results:
   • We calculate the modular exponentiation for various values of $( b )$ (from $1$ to $45$) and store the results.
   • The plot shows $( a^b \mod n )$ on the $y$-$axis$ against the exponent $( b )$ on the $x$-$axis$, illustrating how the result changes as $( b )$ increases.

Visualization

The resulting plot shows the values of $( 3^b \mod 13 )$ for $( b = 1, 2, \ldots, 45 )$:

• The $y$-$axis$ represents the results of $( 3^b \mod 13 )$.
• Notice the cyclic behavior in the values, which is a property of modular arithmetic.

Explanation

This cyclic pattern is essential in cryptography because it enables efficient computations and helps keep the results within manageable bounds.
This property of modular arithmetic forms the basis for cryptographic systems like $RSA$, which rely on modular exponentiation for encryption and decryption.

Applications

Modular exponentiation is widely used in cryptographic systems:

• Public-Key Cryptography: Algorithms like $RSA$ and $Diffie$-$Hellman$ rely on modular exponentiation to ensure secure communication.
• Digital Signatures: Modular exponentiation is also part of signing and verifying digital messages.
• Blockchain: Modular exponentiation is employed in blockchains for secure transaction validation.

Let’s explore a discrete mathematics problem: the Fibonacci sequence, a classic example in combinatorics and discrete mathematics.

Problem: Fibonacci Sequence

The Fibonacci sequence is a series of numbers in which each number (after the first two) is the sum of the two preceding ones.

The sequence starts with 0 and 1, and is defined as follows:

$$
F(0) = 0, \quad F(1) = 1
$$
$$
F(n) = F(n-1) + F(n-2) \quad \text{for } n \geq 2
$$

The sequence goes: $0, 1, 1, 2, 3, 5, 8, 13, 21, 34, …$

The Fibonacci sequence appears in many areas, including biology (e.g., patterns of leaves on a stem) and computer science (e.g., algorithms and data structures).

In this example, we’ll calculate the first $20$ terms of the Fibonacci sequence and visualize the results with a bar graph to illustrate the growth pattern.

Objective

Generate the first $20$ terms of the Fibonacci sequence, then plot the results to observe how the values grow as the sequence progresses.

Python Code

Here’s the Python code to calculate the Fibonacci sequence and visualize it:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24

import matplotlib.pyplot as plt

# Function to generate the Fibonacci sequence
def fibonacci_sequence(n):
    fib_sequence = [0, 1]
    for i in range(2, n):
        fib_sequence.append(fib_sequence[i-1] + fib_sequence[i-2])
    return fib_sequence

# Parameters
n_terms = 20  # Number of terms to generate

# Generate the Fibonacci sequence
fib_sequence = fibonacci_sequence(n_terms)

# Plot the Fibonacci sequence
plt.figure(figsize=(10, 6))
plt.bar(range(n_terms), fib_sequence, color='skyblue')
plt.xlabel('Term Index (n)')
plt.ylabel('Fibonacci Number F(n)')
plt.title('Fibonacci Sequence')
plt.xticks(range(n_terms))
plt.grid(axis='y')
plt.show()

Explanation of the Code

1. Fibonacci Sequence Generation:
   • We define a function fibonacci_sequence that generates the first n terms of the Fibonacci sequence using a loop.
   • Starting with [0, 1], we append each new term as the sum of the two previous terms.
2. Plotting the Sequence:
   • We use a bar graph to visualize the terms in the Fibonacci sequence.
   • The $x$-$axis$ represents the term index $(n)$, and the $y$-$axis$ represents the value of the Fibonacci number at that index.

Visualization

The resulting bar graph displays the first $20$ terms of the Fibonacci sequence:

• The values grow slowly at first but increase rapidly as the index gets higher, showing the characteristic exponential growth pattern of the Fibonacci sequence.
• This growth pattern illustrates how the Fibonacci sequence can approximate exponential growth, even though each term is the sum of only two preceding terms.

Applications

The Fibonacci sequence has various applications:

• Computer Science: Recursive algorithms, dynamic programming.
• Biology: Natural patterns, such as branching in trees and arrangement of leaves.
• Finance: Technical analysis in stock trading.

This example demonstrates how $Python$ can be used to solve discrete math problems and visualize sequences that exhibit interesting growth patterns.

Let’s solve a $pharmacokinetics$ problem using $Python$ and visualize the results.

$Pharmacokinetics$ is the study of how drugs move through the body, and a common problem is calculating drug concentration over time.

Problem: Drug Concentration Over Time (One-Compartment Model)

We’ll assume a single dose of a drug is administered intravenously, and it distributes uniformly throughout a single compartment in the body.

Our goal is to calculate the concentration of the drug in the bloodstream over time and plot it.

Model Assumptions

1. One-compartment model: The drug distributes instantly throughout the body.
2. First-order elimination: The drug is metabolized and eliminated at a rate proportional to its concentration.

Equations and Parameters

The concentration $(C(t))$ of the drug in the bloodstream at time $(t)$ can be calculated with:
$$
C(t) = \frac{D}{V} e^{-k \cdot t}
$$
where:

• $(D)$: The dose administered (in $mg$),
• $(V)$: The volume of distribution (in $L$),
• $(k)$: The elimination rate constant (in $1/h$), calculated as$ (k = \frac{\ln(2)}{t_{1/2}})$,
• $(t_{1/2})$: The drug’s half-life (in hours),
• $(t)$: Time after administration (in hours).

Python Code

Let’s write Python code to calculate and plot the concentration $(C(t))$ over a specified time period.

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25

import numpy as np
import matplotlib.pyplot as plt

# Define parameters
D = 500          # Dose in mg
V = 50           # Volume of distribution in L
t_half = 4       # Half-life of the drug in hours

# Calculate the elimination rate constant
k = np.log(2) / t_half

# Time range (0 to 24 hours)
time = np.linspace(0, 24, 100)

# Calculate concentration over time
concentration = (D / V) * np.exp(-k * time)

# Plot the concentration-time curve
plt.plot(time, concentration, color='blue', label='Drug Concentration')
plt.xlabel('Time (hours)')
plt.ylabel('Concentration (mg/L)')
plt.title('Drug Concentration Over Time')
plt.grid(True)
plt.legend()
plt.show()

Explanation of the Code

1. Parameter Definitions:
   • The dose $(D)$, volume of distribution $(V)$, and half-life $(t_{1/2})$ are defined based on typical values.
2. Elimination Rate Constant:
   • We calculate $(k)$ using the half-life formula $(k = \frac{\ln(2)}{t_{1/2}})$.
3. Concentration Calculation:
   • For each time point in the defined range ($0$ to $24$ hours), we calculate the concentration $(C(t))$ using the exponential decay formula.
4. Plotting:
   • The concentration over time is plotted on a graph with time on the $x$-$axis$ and concentration on the $y$-$axis$.

Visualization

The output graph displays:

• A blue curve representing the drug concentration in the bloodstream over time.
• The concentration decreases exponentially, showing how the drug is metabolized and eliminated from the body.

This model is useful for understanding the $pharmacokinetics$ of a drug after an intravenous dose, helping in dosage planning and ensuring therapeutic levels are maintained without reaching toxic levels.

Let’s solve s geometry problem with $Python$ and visualize the result with a graph.

Here’s the problem:

Problem: Line-Circle Intersection

Given a circle and a line in a 2D plane, find the intersection points of the line and the circle, and then plot them.

Setup

1. Circle: Defined by its center $(x_c, y_c)$ and radius $(r)$.
2. Line: Defined by the linear equation $(y = mx + b)$, where $(m)$ is the slope, and $(b)$ is the y-intercept.

Solution Outline

To find the intersection points:

1. Substitute the line equation $(y = mx + b)$ into the circle equation $(x - x_c)^2 + (y - y_c)^2 = r^2$.
2. Simplify to a quadratic equation in terms of $(x)$.
3. Solve the quadratic equation:
   • If there are no real solutions, the line does not intersect the circle.
   • One real solution indicates a tangent (one intersection point).
   • Two real solutions indicate two intersection points.
4. Plot the circle, line, and intersection points.

Python Code

Here’s the $Python$ code for finding the intersection points and plotting them.

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61

import numpy as np
import matplotlib.pyplot as plt

def find_line_circle_intersections(xc, yc, r, m, b):
    # Solve for x in the quadratic equation
    A = 1 + m**2
    B = -2 * xc + 2 * m * (b - yc)
    C = xc**2 + (b - yc)**2 - r**2
    
    discriminant = B**2 - 4 * A * C
    
    # If discriminant is negative, no real intersection
    if discriminant < 0:
        return "No intersection"
    
    # Calculate the x-coordinates of intersection points
    x1 = (-B + np.sqrt(discriminant)) / (2 * A)
    x2 = (-B - np.sqrt(discriminant)) / (2 * A)
    
    # Calculate the corresponding y-coordinates
    y1 = m * x1 + b
    y2 = m * x2 + b
    
    return [(x1, y1), (x2, y2)] if discriminant > 0 else [(x1, y1)]

# Define the circle and line parameters
xc, yc, r = 0, 0, 5  # Circle with center (0, 0) and radius 5
m, b = 1, 1          # Line with slope 1 and y-intercept 1

# Find the intersection points
intersections = find_line_circle_intersections(xc, yc, r, m, b)

# Set up the plot
theta = np.linspace(0, 2 * np.pi, 100)
x_circle = xc + r * np.cos(theta)
y_circle = yc + r * np.sin(theta)

# Plot the circle
plt.plot(x_circle, y_circle, label='Circle', color='blue')

# Plot the line
x_vals = np.linspace(-r-2, r+2, 100)
y_vals = m * x_vals + b
plt.plot(x_vals, y_vals, label='Line', color='red')

# Plot the intersection points, if any
if intersections != "No intersection":
    for (x, y) in intersections:
        plt.plot(x, y, 'go')  # Green dots for intersection points
        plt.text(x, y, f'({x:.2f}, {y:.2f})', ha='right')

# Add labels and title
plt.xlabel('x')
plt.ylabel('y')
plt.axhline(0, color='black',linewidth=0.5)
plt.axvline(0, color='black',linewidth=0.5)
plt.title('Intersection of a Circle and a Line')
plt.legend()
plt.grid(True)
plt.axis('equal')
plt.show()

Explanation of the Code

1. Intersection Calculation:
   • The equation of the circle $(x - x_c)^2 + (y - y_c)^2 = r^2$ and the line equation $y = mx + b$ are combined.
   • Solving the resulting quadratic equation allows us to find the intersection points (if they exist).
2. Plotting:
   • The circle is drawn by parameterizing it with $(\theta)$.
   • The line is plotted using a range of $(x)$ values.
   • The intersection points (if any) are marked with green dots and labeled with their coordinates.

Visualization

The output graph displays:

• The circle and line in blue and red, respectively.
• Green dots indicate the intersection points with coordinates labeled for easy interpretation.

This graph illustrates the intersection of a circle and a line in a 2D plane.

Key Elements:

1. Circle (blue line): Centered at the origin $(0, 0)$ with a radius of $5$.
2. Line (red line): Represented by the equation $(y = x + 1)$, which has a slope of $1$ and a y-intercept of $1$.

Intersection Points:

The green dots mark the two intersection points where the line crosses the circle:

• Point 1: Located at approximately $(-4.00, -3.00)$.
• Point 2: Located at approximately $(3.00, 4.00)$.

Each point is labeled with its coordinates, making it easy to see the exact locations of intersection.

This visualization effectively shows where the line and circle meet, providing a clear solution to the problem of finding their intersection points.

To solve an optimization problem with a nonlinear objective function or nonlinear constraints in $Python$, we can use libraries like SciPy’s optimize module.

SciPy provides various functions to handle nonlinear optimization problems, including minimize with methods like 'SLSQP' (Sequential Least Squares Programming), which can handle nonlinear constraints.

In this example, we’ll set up a nonlinear objective function with a nonlinear constraint and solve it using SciPy.

Problem Setup

Suppose we want to minimize the following objective function:

$$
f(x, y) = (x - 2)^2 + (y - 1)^2
$$

This objective function represents the distance between $(x, y)$ and the point $(2, 1)$, which we aim to minimize.
We’ll also define a nonlinear constraint:

$$
g(x, y) = x^2 + y^2 - 1 \leq 0
$$

This constraint forces $(x, y)$ to lie within the unit circle centered at the origin.

Code Implementation

The following $Python$ code demonstrates how to define and solve this nonlinear optimization problem using SciPy.

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35

from scipy.optimize import minimize, NonlinearConstraint
import numpy as np

# Define the nonlinear objective function
def objective(vars):
    x, y = vars
    return (x - 2)**2 + (y - 1)**2

# Define the nonlinear constraint function
def constraint(vars):
    x, y = vars
    return x**2 + y**2 - 1  # This should be <= 0 to satisfy the constraint

# Create a nonlinear constraint object
nonlinear_constraint = NonlinearConstraint(constraint, -np.inf, 0)

# Initial guess for the variables
initial_guess = [0, 0]

# Solve the optimization problem
result = minimize(
    objective,
    initial_guess,
    method='SLSQP',
    constraints=[nonlinear_constraint]
)

# Display the results
if result.success:
    print("Optimal solution found:")
    print("x =", result.x[0])
    print("y =", result.x[1])
    print("Minimum value of objective function =", result.fun)
else:
    print("Optimization failed:", result.message)

Explanation of Key Components

1. Objective Function:

   • The function objective(vars) takes an array vars with variables $( x )$ and $( y )$ and returns the value of the objective function $(x - 2)^2 + (y - 1)^2$.
     This is the function we want to minimize.

2. Constraint Function:

   • The function constraint(vars) defines the constraint $( x^2 + y^2 - 1 )$, which must be less than or equal to zero for a valid solution.
     In other words, the solution must lie within or on the unit circle.

3. Nonlinear Constraint:

   • We use NonlinearConstraint to define the constraint, setting the bounds to -np.inf (no lower bound) and 0 (upper bound), which enforces $( x^2 + y^2 - 1 \leq 0 )$.

4. Optimization Solver:

   • The minimize function from SciPy is used to find the minimum of the objective function.
   • We specify method='SLSQP' to handle the nonlinear constraint, and we pass the constraint as a list.
   • The initial guess is $([0, 0])$, which serves as a starting point for the algorithm.

5. Output:

   • If successful, result.x contains the optimal values for $( x )$ and $( y )$, and result.fun gives the minimum value of the objective function.
   • If the optimization fails, an error message is displayed.

Result

Optimal solution found:
x = 0.8944272288069395
y = 0.4472135198861174
Minimum value of objective function = 1.5278640450001992

The optimization was successful, yielding an optimal solution for $( x )$ and $( y )$ within the specified constraints. The optimal values are:

• $( x = 0.894 )$
• $( y = 0.447 )$

These values minimize the objective function, achieving a minimum value of approximately $1.528$.

This result satisfies the constraint that $( x^2 + y^2 \leq 1 )$, meaning the solution lies within the unit circle as required.

Summary

This example demonstrates how to solve a constrained nonlinear optimization problem in $Python$ using SciPy.

The SLSQP algorithm in SciPy‘s minimize function is well-suited for handling nonlinear constraints and is effective for finding solutions when either the objective function or constraints are nonlinear.

Creating a 3D bubble chart in $Python$ using plotly can effectively visualize data in three dimensions, with the addition of varying bubble sizes to represent a fourth data dimension.

In this example, we’ll use plotly.graph_objects to create a Scatter3d plot where each point (or bubble) has x, y, and z coordinates along with a size attribute for the bubble’s radius.

Here’s how to create a detailed 3D bubble chart with plotly.

Step-by-Step Code Explanation

1. Generate Data: We’ll generate synthetic data for x, y, and z coordinates.
   We’ll also create a size variable to determine the size of each bubble.
2. Create a Scatter3d Plot: Using plotly.graph_objects, we’ll create a Scatter3d object to represent each bubble, where:
   • x, y, and z correspond to the spatial coordinates.
   • marker.size controls the radius of each bubble.
   • marker.color is set to vary based on another variable for an additional data layer.
3. Customize Layout: We’ll customize the layout to make the chart more informative and visually appealing.

Code Example

Here’s how you can implement a 3D bubble chart using plotly.

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38

import numpy as np
import plotly.graph_objects as go

# Generate synthetic data
np.random.seed(0)
num_points = 100
x = np.random.normal(loc=0, scale=1, size=num_points)
y = np.random.normal(loc=1, scale=2, size=num_points)
z = np.random.normal(loc=2, scale=1.5, size=num_points)
size = np.random.randint(5, 20, size=num_points)  # Bubble sizes (can represent some other variable)

# Create a 3D scatter plot for the bubble chart
fig = go.Figure(data=[go.Scatter3d(
    x=x,
    y=y,
    z=z,
    mode='markers',
    marker=dict(
        size=size,             # Bubble sizes
        color=z,               # Color based on Z-axis values or any other metric
        colorscale='Viridis',  # Color scale for bubbles
        opacity=0.8,
        sizemode='diameter'    # Control size mode (use 'diameter' for relative sizing)
    )
)])

# Update layout for clarity
fig.update_layout(
    title="3D Bubble Chart",
    scene=dict(
        xaxis_title="X Axis",
        yaxis_title="Y Axis",
        zaxis_title="Z Axis"
    ),
    margin=dict(l=0, r=0, b=0, t=30)  # Adjust margins for better view
)

fig.show()

Explanation of Key Components

1. Data Generation: We generate random data for x, y, and z coordinates, along with a size array for bubble sizes.
   This size variable can represent any additional metric, such as population, magnitude, or frequency.

2. 3D Scatter Plot with Scatter3d:

   • x, y, and z represent the spatial coordinates.
   • marker.size controls the size of each bubble.
   • marker.color uses z values (or another variable) to set bubble colors, adding a layer of information.

3. Customization:

   • Color Scale: The Viridis color scale is used for a visually pleasing gradient effect based on z values.
   • Opacity: Setting opacity=0.8 allows overlapping bubbles to blend slightly, enhancing the 3D effect.
   • Layout: Axes are labeled, and margins are adjusted for an optimal view.

Summary

This 3D bubble chart effectively visualizes data across three spatial dimensions (x, y, and z), with the size of each bubble representing a fourth dimension.
This visualization is useful for showing multi-dimensional data, especially when analyzing clusters, correlations, or densities across multiple variables.

Visualizing complex mathematical surfaces in 3D can be incredibly insightful, especially when using interactive plotting libraries like Plotly.

In this tutorial, we’ll build an intricate 3D contour plot using plotly.graph_objects and Python’s numpy.

Although Plotly doesn’t directly support 3D contour plots, we can achieve a similar effect using the go.Surface plot type.

Let’s walk through the steps to create this visualization.

Step 1: Define a Complex Mathematical Function

To create a compelling 3D plot, we need a mathematical function that has varied, interesting features.

For this example, let’s combine trigonometric functions with an exponential decay factor.

This will produce wave-like patterns that decay outward, creating a complex and visually engaging surface.

Here’s the function we’ll use:

1
2
3
4

import numpy as np

def complex_function(x, y):
    return np.sin(x**2 + y**2) * np.exp(-0.1 * (x**2 + y**2)) + 0.5 * np.sin(3 * x) * np.cos(3 * y)

This function combines sine waves with exponential decay, resulting in a surface that has oscillations and a natural fade-out effect, adding to the complexity.

Step 2: Create a Grid of Points

We need to evaluate our function across a grid of $(x)$ and $(y)$ values.

Using numpy.linspace, we can create a finely spaced grid.

This allows us to generate smooth contours and ensures that the details of our surface are captured.

1
2
3
4

x = np.linspace(-5, 5, 100)
y = np.linspace(-5, 5, 100)
X, Y = np.meshgrid(x, y)
Z = complex_function(X, Y)

Here, X and Y form a 2D grid, and Z is the output of our function across this grid, representing height.

Step 3: Plot the Surface with Plotly’s go.Surface

To plot in 3D, we use Plotly’s go.Surface plot type.

go.Surface is flexible, allowing us to visualize a 3D surface with customizable color maps and contour lines.

We’ll set up contour lines along the z axis to give the appearance of a 3D contour plot.

1
2
3
4
5
6
7
8
9
10
11

import plotly.graph_objects as go

fig = go.Figure(data=go.Surface(
    x=X,
    y=Y,
    z=Z,
    colorscale='Viridis',
    contours={
        "z": {"show": True, "size": 0.1, "start": -1, "end": 1} # Define contour levels
    }
))

In this code:

• colorscale='Viridis' applies a perceptually uniform color scale, enhancing depth perception.
• contours settings add contour lines along the z axis, creating a layered effect that resembles contour lines on a topographic map.

Step 4: Customize the Layout

To make the plot more informative and visually appealing, we’ll add labels and adjust the camera view to highlight the surface’s details.

1
2
3
4
5
6
7
8
9
10
11
12

fig.update_layout(
    title="Complex 3D Contour Plot using Surface",
    scene=dict(
        xaxis_title="X Axis",
        yaxis_title="Y Axis",
        zaxis_title="Z Axis",
        camera=dict(
            eye=dict(x=1.25, y=1.25, z=0.75) # Adjust camera position for better view
        )
    )
)
fig.show()

With these settings:

• xaxis_title, yaxis_title, and zaxis_title label each axis.
• camera positions the viewpoint, giving a balanced, 3D perspective on the plot.

Full Code

Here’s the complete code to generate this 3D contour-style plot:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38

import numpy as np
import plotly.graph_objects as go

# Define the function for complex surface
def complex_function(x, y):
    return np.sin(x**2 + y**2) * np.exp(-0.1 * (x**2 + y**2)) + 0.5 * np.sin(3 * x) * np.cos(3 * y)

# Generate a grid of points
x = np.linspace(-5, 5, 100)
y = np.linspace(-5, 5, 100)
X, Y = np.meshgrid(x, y)
Z = complex_function(X, Y)

# Create 3D surface plot with contours on z-axis
fig = go.Figure(data=go.Surface(
    x=X,
    y=Y,
    z=Z,
    colorscale='Viridis',
    contours={
        "z": {"show": True, "size": 0.1, "start": -1, "end": 1}
    }
))

# Customize plot
fig.update_layout(
    title="Complex 3D Contour Plot using Surface",
    scene=dict(
        xaxis_title="X Axis",
        yaxis_title="Y Axis",
        zaxis_title="Z Axis",
        camera=dict(
            eye=dict(x=1.25, y=1.25, z=0.75)
        )
    )
)

fig.show()

Output

Explanation of Key Components

• Complex Function: By combining trigonometric and exponential functions, the surface has both oscillations and decay, providing an intricate structure.
• Surface Plot and Contours: Although Plotly lacks a dedicated 3D contour plot type, adding contours on the z axis of a go.Surface plot effectively mimics this visualization style.
• Interactive Plotting: Plotly’s interactive nature allows users to explore the plot by rotating and zooming, revealing intricate details from different perspectives.

This plot is highly customizable.

Try experimenting with different functions, contour settings, and color scales to create unique 3D visualizations!

This plot simulates circular waves radiating outward from the origin.

The amplitude decays as the waves spread, creating a visually dynamic $3D$ $mesh$.

1. Import Required Libraries

First, import the necessary libraries:

1
2

import numpy as np
import plotly.graph_objects as go

2. Define the Circular Wave Grid

To create circular waves, we define a grid in polar coordinates (r for radial distance and theta for angle) and convert it to Cartesian coordinates (x and y).

This grid will form a circular, ripple-like structure.

1
2
3
4
5
6
7
8

# Define the polar grid
r = np.linspace(0, 10, 100)
theta = np.linspace(0, 2 * np.pi, 100)
r, theta = np.meshgrid(r, theta)

# Convert to Cartesian coordinates
X = r * np.cos(theta)
Y = r * np.sin(theta)

3. Define the Circular Wave Function for Z Coordinates

Define a wave function that decreases in amplitude as it moves outward.

We can use a sine function with a decaying exponential to achieve this effect:

1
2

# Circular wave function
Z = np.sin(3 * r - theta) * np.exp(-0.3 * r)

4. Reshape Data for the Mesh Plot

Flatten X, Y, and Z to 1D arrays for Plotly’s go.Mesh3d.

1
2
3

X = X.flatten()
Y = Y.flatten()
Z = Z.flatten()

5. Create the Mesh Plot

Now, create the $3D$ $mesh$ plot with customized color and transparency:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24

# Create the mesh plot for the circular wave
fig = go.Figure(
    data=[go.Mesh3d(
        x=X, y=Y, z=Z,
        color='blue',
        opacity=0.7,
        intensity=Z,
        colorscale="Cividis",
        alphahull=5
    )]
)

# Set layout options for a better view
fig.update_layout(
    title="3D Circular Wave Mesh Plot",
    scene=dict(
        xaxis_title='X-axis',
        yaxis_title='Y-axis',
        zaxis_title='Z-axis',
        camera=dict(
            eye=dict(x=1.5, y=1.5, z=1.5)
        )
    )
)

6. Display the Plot

To display the circular wave plot, use:

1

fig.show()

Explanation of Key Parameters

• Wave Equation: The function np.sin(3 * r - theta) * np.exp(-0.3 * r) simulates oscillations with amplitude decay, creating a wave pattern that radiates outward.
• Colorscale: The Cividis scale adds contrast, highlighting the wave peaks and valleys for a dramatic effect.
• Opacity: Lower opacity adds depth to the plot, making the wave’s radial decay easier to perceive.

This creates a visually engaging $3D$ $mesh$ that mimics ripples in water, showcasing the decaying wave amplitude as it moves outward from the origin.