Difference equation my answer different from books

Click For Summary
SUMMARY

The discussion centers on solving a difference equation related to drug concentration in the bloodstream, specifically modeled by the equation y[n] = x[n] + αy[n-1], where α represents the drug's decay rate. The user attempted to derive a solution using methods like undetermined coefficients but arrived at a different expression than the textbook, which provides y[n] = (1 - α^(n+1)) / (1 - α) through iteration and geometric series recognition. The discrepancy arises from the user's approach, which does not fully account for the iterative nature of the problem and the need to determine the arbitrary constant K.

PREREQUISITES
  • Understanding of difference equations and their applications in modeling
  • Familiarity with geometric series and their summation formulas
  • Knowledge of methods for solving linear difference equations, including undetermined coefficients
  • Basic concepts of pharmacokinetics, particularly drug concentration dynamics
NEXT STEPS
  • Study the method of solving linear difference equations using iteration
  • Learn about geometric series and their applications in mathematical modeling
  • Explore pharmacokinetic models to understand drug absorption and elimination
  • Review the method of undetermined coefficients for solving difference equations
USEFUL FOR

This discussion is beneficial for students and professionals in mathematics, pharmacology, and biomedical engineering, particularly those interested in modeling drug dynamics and solving difference equations.

fahraynk
Messages
185
Reaction score
5

Homework Statement


When drugs are used to treat a medical condition, doctors often recommend starting with a higher dose on the first day than on subsequent days. In this problem, we consider a simple model to understand why. Assume that the human body is a tank of blood and that drugs instantly dissolve in the blood when ingested. Further assume that the drug vanishes from the body at a rate that is proportional to drug concentration.

Let X[n] represent the amount of drug taken on day n, and let y[n] represent the total amount of drug in the blood on day n, just after the dos x[n] has dissolved int eh blood, so that :
y[n]=x[n]+##\alpha##[n-1].

Assume that no drug is in the blood before day 0, and that one unit of drug is taken each day, starting with day 0. Determine an expression for the amount of drug in the blood immediately after the dose on day n has dissolved.

Homework Equations

The Attempt at a Solution



I plug ##y=P^n## into the difference equation to get ##y_h=K\alpha^n## and I plug ##y=C## into the equation for the particular solution and get for an answer :
## y = K\alpha^n - \frac{x(n)}{\alpha-1}##

The book says this:
Solve by iteration :
##n : y[n]##
##0 : 1##
##1:1+\alpha##
##2:1+\alpha+\alpha^2##
##3:1+\alpha+\alpha^2+\alpha^3##
##...## : ## ...##
##n : \frac{1-\alpha^{n+1}}{1-\alpha}##
Books answer is : ##\frac{1-\alpha^{n+1}}{1-\alpha}##

I think the book is using a power series representation or something like that. Is my answer wrong for x(n)=1 ? If so what am I doing wrong ? I rather not use a power series if I can answer the difference equation using methods I already know like undetermined coefficients.
 
Physics news on Phys.org
fahraynk said:
y[n]=x[n]+α[n-1].
Do you mean: ##y[n] = x[n] + \alpha y[n-1]## where ##\alpha \in (0,1)## determines the rate at which the drug disappears from the blood?

fahraynk said:
I think the book is using a power series representation or something like that.
It is calculating the first few terms of the solution, recognizing a geometric series and then using the standard formula for the partial sums of such a series.
 
fahraynk said:
I plug ##y=P^n## into the difference equation to get ##y_h=K\alpha^n## and I plug ##y=C## into the equation for the particular solution and get for an answer :
## y = K\alpha^n - \frac{x(n)}{\alpha-1}##
I'm not sure why you're writing the solution in terms of x(n) since you've assumed x(n) is a constant to find the particular solution. Just write the constant.

You didn't finish solving the problem. You must still determine the arbitrary constant K.
 

Similar threads

  • · Replies 1 ·
Replies
1
Views
2K
Replies
20
Views
4K
Replies
1
Views
1K
Replies
3
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
Replies
1
Views
2K
  • · Replies 11 ·
Replies
11
Views
3K
Replies
5
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
Replies
2
Views
2K