Artifact 9bf68f04127d9c8248dae66c1ff47d0a60bf8133949a4ceb639297da4449bd33:
- Executable file
r36/XMPL/RANDPOLY.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: 1588) [annotate] [blame] [check-ins using] [more...]
% randpoly.tst % F.J.Wright@Maths.QMW.ac.uk, 14 July 1994 off allfac; on div, errcont; % Univariate: % ---------- randpoly x; % Equivalent to above: randpoly {x}; randpoly(x, dense); % univariate default already dense randpoly(x, degree=10, ord=5); % Bivariate: % --------- % Default is sparse randpoly {x,y}; randpoly({x,y}, dense); randpoly({x,y}, degree=10); % Lots of terms: randpoly({x,y}, dense, degree=10); randpoly({x,y}, dense, degree=10, ord=5); % Sparse: randpoly({x,y}, deg=10, ord=5); % Dense again: randpoly({x,y}, terms=1000, maxdeg=10, mindeg=5); % Exponent and coefficient functions: % ---------------------------------- randpoly({x,y}, expons = rand(-10 .. 10)); % Trivial example: randpoly({x,y}, expons = proc 5); randpoly({x,y}, expons = proc(2*random(0 .. 5))); randpoly({x,y}, coeffs = rand(-999 .. 999)); procedure coe; randpoly(a, terms=2)$ randpoly({x,y}, coeffs = coe); randpoly({x,y}, coeffs = coe, degree = 10); % Polynomials composed with general expressions: % --------------------------------------------- randpoly({x,y^2}); randpoly(x^2 - y^2); % This should give the constant term: sub(x=y, ws); randpoly({x^2 - a^2, y - b}); % This should give the constant term: sub(x=a, y=b, ws); % Polynomials with specified zeros: % -------------------------------- randpoly(x = a); % This should give 0: sub(x=a, ws); randpoly({x = a, y = b}); % This should give 0: sub(x=a, y=b, ws); % Invalid input detection: % ----------------------- randpoly({x,y}, degree=foo); randpoly({x,y}, foo); randpoly({x,y}, degree=-5); on allfac; off div, errcont; end;