## Total Profit

A company finds that the marginal profit, in dollars, from drilling a well that is $$x$$ feet deep is given by

$P^\prime(x) = \sqrt[5]{x}$ Find the profit when a well 250 ft deep is drilled.

import sympy as sp
from sympy.abc import x
import matplotlib.pyplot as plt
import numpy as np

marginal_profit = sp.root(x, 5)
marginal_profit
## x**(1/5)
profit = sp.integrate(marginal_profit, x)
profit
## 5*x**(6/5)/6
p = sp.lambdify(x, profit)
p(250)
## 628.5600350567877
sp.integrate(marginal_profit, (x, 0, 250)).evalf()
## 628.560035056788
x_values = np.linspace(0, 300)

plt.plot(x_values, p(x_values))
plt.fill_between(x_values, p(x_values), where=[x <= 250 for x in x_values], color="green", alpha=0.3)

## Increasing Total Profit

Laso Industries finds that the marginal profit, in dollars from the sale of $$x$$ digital control boards is given by

$P^\prime(x) = 2.6x^{0.1}$ A customer orders 1200 digital control boards, and later increases the order to 1500. Find the extra profit resulting from the increase in order size.

marginal_profit = (2.6 * x)**0.1

profit = sp.integrate(marginal_profit, x)
profit
## 1.00024099373274*x**1.1
second_purchase = sp.integrate(marginal_profit, (x, 0, 1500))
second_purchase
## 3117.48975261839
first_purchase = sp.integrate(marginal_profit, (x, 0, 1200))
first_purchase
## 2438.95630774585
extra_profit = second_purchase - first_purchase
extra_profit
## 678.533444872540
x_values = np.linspace(0, 2000)

p = sp.lambdify(x, profit)
plt.plot(x_values, p(x_values))
plt.fill_between(x_values, p(x_values), where=[x <= 1200 for x in x_values], color="green", alpha=0.3)
plt.fill_between(x_values, p(x_values), where=[(x <= 1500) and (x > 1200) for x in x_values], color="red", alpha=0.3)

## Accumulated Sales

Ragss, Ltd. estimates that its sales will grow continuously at a rate given by

$S^\prime(t) = 10e^t$ where $$S^\prime(t)$$ is the rate at which sales are increasing, in dollars per day, on day $$t$$.

• Find the accumulated sales for the first 5 days.
• Find the sales from the 2nd day through the 5th day.
from math import e
from sympy.abc import t

marginal_sales = (10 * e)**t
marginal_sales
## 27.1828182845905**t
sales = sp.integrate(marginal_sales, t)
sales
## 0.302793106564114*27.1828182845905**t
sp.integrate(marginal_sales, (t, 0, 5)).evalf()
## 4493847.84717322
sp.integrate(marginal_sales, (t, 2, 5)).evalf()
## 4493624.41444125