Artifact 9860ceeb94b5f4faf0fc613544ce99bdc9d9f8d6f3bc0fc833b2062c2646f97e:
- Executable file
r38/packages/odesolve/extend.tst
— part of check-in
[f2fda60abd]
at
2011-09-02 18:13:33
on branch master
— Some historical releases purely for archival purposes
git-svn-id: https://svn.code.sf.net/p/reduce-algebra/code/trunk/historical@1375 2bfe0521-f11c-4a00-b80e-6202646ff360 (user: arthurcnorman@users.sourceforge.net, size: 2293) [annotate] [blame] [check-ins using] [more...]
% Test and demonstration of the ODESolve extension interface % F.J.Wright@Maths.QMW.ac.uk, Time-stamp: <17 July 2000> % Load odesolve before inputting this file if not using test interface: % load_package odesolve; % Hook into the general ODE solver: algebraic procedure ODESolve_Hook_Demo (ode, y, x); %% For any ODE, if the dependent variable is z then this hook %% procedure returns a solution corresponding to ODESolve failing %% to find any solution; otherwise it returns nil (nothing) and so %% is ignored. if y=z then {ode=0}; % Set the hook: symbolic(ODESolve_Before_Hook := '(ODESolve_Hook_Demo)); % Hook into the nonlinear ODE solver: algebraic procedure ODESolve_Non_Hook_Demo (ode, y, x, n); %% If the ODE is nontrivially nonlinear and the order is 3 then %% this hook procedure returns a solution corresponding to ODESolve %% failing to find any solution; otherwise it returns nil (nothing) %% and so is ignored. if n=3 then {ode=0}; % Set the hook: symbolic(ODESolve_Before_Non_Hook := '(ODESolve_Non_Hook_Demo)); % Hook into the general linear ODE solver: algebraic procedure ODESolve_Lin_Hook_Demo (odecoeffs, driver, y, x, n, m); %% If the ODE is linear and the order is 3 then this hook procedure %% returns a solution corresponding to ODESolve failing to find any %% solution; otherwise it returns nil (nothing) and so is ignored. %% (NB: Algebraic-mode lists are indexed from 1 in REDUCE!) if n=3 then {(for i := m : n sum part(odecoeffs,i+1)*df(y,x,i)) = driver}; % Set the hook: symbolic(ODESolve_Before_Lin_Hook := '(ODESolve_Lin_Hook_Demo)); % Test all the hooks: % The general ODE solver: odesolve(df(y,x)); % hook ignored odesolve(df(z,x)); % hook operates % The nonlinear ODE solver: odesolve(y*df(y,x,2)+1); % hook ignored odesolve(y*df(y,x,3)+1); % hook operates % The general linear ODE solver: odesolve(df(y,x,2)+1); % hook ignored odesolve(df(y,x,3)+1); % hook operates % Clear the hooks: symbolic(ODESolve_Before_Hook := nil); symbolic(ODESolve_Before_Non_Hook := nil); symbolic(ODESolve_Before_Lin_Hook := nil); end;