% ODMP
% L.Linares, 2008
% For a second order differential equation with 
% double real Eigenvalues b1 b2, and knowing
% the initial conditions:
%   yo = initial value of the function of time f(0).
%   yo'= initial value of the derivative of f'(0).  
%   yf = final value of the function of time f(infin).
%   b1 = one Eigenvalue
%   b2 = the other Eigenvalue
% NOTE: In the code, * means multiply.
% It produces k1 and k2 so:
%  y(t) = k1.exp(b1.t)+k2.exp(b2.t)+yf

<< "yo Dyo yf b1 b2" PROMPT
   ->  a b K3 b1 b2
   << '(b2*K3-a*b2+b)/(b1-b2)'
      EVAL ->STR  "k1=" SWAP +
      '(a*b1-b-b1*K3)/(b1-b2)'
      EVAL ->STR  "k2=" SWAP +
   >>
>>