clear
disp('Application to lead poisoning in blood, tissue and bones')
format long
A = [ -.0361 .0124 .000035; .0111 -.0286  0; .0039 0 -.000035]
pause
f = [ 49.3 0 0]'
pause
xp_sol = -A\f
pause
[u d] = eig(A)
pause
d
pause
x0 = [0 0 0]'
pause
c = u\(x0 -  xp_sol)
pause
syms t
disp('x = c(1)*u(:,1)*exp(d(1,1)*t) + c(2)*u(:,2)*exp(d(2,2)*t) +c(3)*u(:,3)*exp(d(3,3)*t) + xp_sol')
pause
x = c(1)*u(:,1)*exp(d(1,1)*t) + c(2)*u(:,2)*exp(d(2,2)*t) +c(3)*u(:,3)*exp(d(3,3)*t) + xp_sol
pause
ezplot(x(1),[0 80000])
hold on
ezplot(x(2),[0 80000])
ezplot(x(3),[0 80000])