Thomas Richard

Mr. Thomas Richard

3407 Reputation

13 Badges

14 years, 76 days
Maplesoft Europe GmbH
Technical professional in industry or government
Aachen, North Rhine-Westphalia, Germany

MaplePrimes Activity


These are answers submitted by Thomas Richard

The r.h.s. of your solution (multiplied by a constant factor) can be obtained via

simplify(sol,symbolic);

which is also confirmed by odetest.

Your initial example cannot be correct; please plot the hypergeom and the elementary representations. Maybe some typo or copy&paste error (which I often make...)?

Three of the four in your list are covered by

convert(%,StandardFunctions);

 

The recommended tool for evaluating algebraic numbers is evala. However, some preprocessing is needed for this expression:

evala(factor(expr));

In theory, factoring should not be necessary. So that's a weakness, I'd say.

In the definition of C2, you need to remove all the spaces so that it's recognized as a floating-point number, just like the other parameters. Otherwise it will be interpreted as an expression involving the unbound symbol e, which is syntactically valid - thus no error message appears.

The solve call will provide the approximate roots then, without resorting to RootOf representations. And no need to call allvalues.

You almost had it correct; essentially only the curly braces around Cond had to be removed. Some more cleanup was done in the attached file.

exact_solution_error_TR.mw

Edit: two errors of mine fixed

You could strip off the first column by appending [..,2..] to the definition of TR. However, that is ignored by the summarize=embed option. To work around, omit the option and do the embedding in a separate step:

TR := FrequencyTable(Obs, weights = Eff, headers, tableweights = [4, 2, 2, 2, 2])[..,2..]:
DocumentTools:-Tabulate(TR):

By the way, some previous steps can be shortened considerably, e.g.

Obs := lhs~(L):
Eff := rhs~(L):

If we supply one more item in the first list, it will work right away.

And for the second list, did you mean to start with 0, or perhaps 2? If I choose the latter, a simple result is found:

restart:
with(gfun):
List1 := [0,4,16,36,64,100,144]:
List2 := [2,3,6,11,18,27,38]:

R1 := listtorec(List1,f1(n));
rsol1 := rsolve(op(1,R1),f1(n));
simplify(rsol1);

R2 := listtorec(List2,f2(n));
rsol2 := rsolve(op(1,R2),f2(n));
simplify(rsol2);

 

Try this function definition instead:

phi := (x, y) -> Transpose(Multiply(Miv,<x,y,1>));

 

You can also obtain a series solution of this IVP as follows:

sys_ivp := [sys_ode, ics]:
vars := [x(t),y(t)]:
sols := dsolve(sys_ivp, vars, 'series');
solp := convert(sols, 'polynom');
plot(map(rhs,solp),t=0..10);

 

Can your colleague install the free Maple Player?

It should display your Maple file(s) properly, and you wouldn‘t need to export to PDF.

The simplify command needs an option here (perhaps it was applied by default in that older version, but I doubt that):

simplify(f,'wronskian');

This will return 1 right away.

1. Check out the Maple Calculator App which is available for iOS and Android. It's fairly powerful, though I don't know the Sharp calculator you mentioned.

2. Please specify what you mean by "using" here. What do you want to emulate? Maple platform support is covering the three mainstream operating systems: Windows, Linux, macOS.

I have seen that error in a few runs on a Windows computer, but not easy to reproduce. On a Mac, it returns the generators quite fast (half a minute or so), and more reliably. I've uploaded my worksheet with saved output.
Both computers are running 2023.2.1, and I don't have any earlier patch levels installed.

When you submit an SCR, please include the version of Physics updates installed on your system.

symgen-2023.2.1.mw

By default, it should be seen when ordering. You should have received an email as well (so check your SPAM filter, please).

If it is still not found, please contact our Customer Service department.

Short explanation: it would break Maple's syntax rules.

1 2 3 4 5 6 7 Last Page 2 of 44