How i convert the following mathematica code to Maple

f1[t^m_] := Gamma [m + 1] t^(m - a[t]) /Gamma[m + 1 + a[t]];
f1[t] := Gamma [2] t^(1 - a[t]) /Gamma[2 + a[t]];

f1[p_Plus] := Map[f1, p]
f1[c_*f_] := c*f1[f] /; FreeQ[c, t]
f1[c_] := 0 /; FreeQ[c, t]

f1[x^2*t^2 + t + x]

gives

t^(1 - a[t])/Gamma[2 + a[t]] + (2 t^(2 - a[t]) x^2)/Gamma[3 + a[t]]

Glyph2.maple

The above link is a maple workbook containing the package Glyph2 together with the install script and two worksheets illustrating the difficulty I am having with the code in the procedure assignaltbasis (near the end of the code in Glyph2.mpl).

The problematic code is shown standalone is basisdemo worksheet, but reproduced below for reference. Interestingly, even though the table should be evaluated in a seq command, outside the module it works without eval. Inside the module evaluation does not make any difference.

indxs := [1, 2, 3, 4];
altbasis := [seq(seq(e[i[]], i = combinat:-choose(indxs, 2*k)), k = 1 .. 2)];
gen := table([seq(op(gen[i]) = bivec[i], i = 1 .. nops(gen))]);
exprs := [seq(exteriorp(seq(gen[op(2*i - 1 .. 2*i, b)], i = 1 .. nops(b)/2)), b = altbasis)];

The worksheet Basischangedemo shows how the exported procedure “assignaltbasis” is used in the module and the error that results if a bivector basis is used. From experimentation I have found that the error results because “assignaltbasis” passes the table name to the procedure “exteriorp” instead of the bivectors that should have resulted from evaluation of the table. Including the eval command within the seq procedure does not make any difference to the result. I am unclear whether this is a problem with evaluation of the table or a problem that results from namespace conflict.

I have a further namespace problem with the local variables, AltRoot, AltToStd, and StdToAlt which store the results from assignaltbasis.  If I declare these variables as exports, I get an error message “attempt to assign to a protected variable”. Therefore I added an exported procedure “ShowAltBasis” to display the results.  As you can see, this procedure produces the fully-qualified names of each blade making the result difficult to read. I would prefer to get rid of the “Glyph2:-” qualification.

The waters are muddied even further because the other worksheet, “basisdemo”, contains the core code extracted from assignaltbasis showing that the code using the table actually works (outside the module). From this I conclude that a namespace conflict exists within the module, preventing the table from fully evaluating.

My understanding is that within the module local and exported names are more or less equivalent, so I am puzzled why local and exported names in assignaltbasis are so problematic.

I would appreciate any advice about these issues. Clearly my understanding of maple’s namespaces is incomplete.

Plot([x^2, subs(x=0, 2*x)], axis[1]=[color=white]);

can not find the y=0 this red line

how to plot y=x^2 and x = 0 on the same graph to show the slope at x=0 ?

How to derive from n*x^(n-1) to limit definition ?

I am trying to solve a simple two-equation linear system with solve, but I keep getting this weird result where the magnitude order of the numerator and denominator don't cancel out. Whenever I have a linear system with floating-point numbers, more unknowns than equations and try to solve it for a specific set of variables, this happens.

Any suggestions on how to get around this? Or do I need to solve it manually?

How to proof exp(x) is transcendental function without using calculus?

How to write this equation in maple?

g(k-r-1);

conditions:=

1.r from 0 to k, 

2.g(0)=1;

3.g(n)=0, where n not equal to 0.

How to convert diff(f(t),t) + diff(g(t),t)

to diff(t(f,g), f)+ diff(t(f,g), g)

?

How to cmaple or mw file maple dualaxisplot output to jpeg file after press and run?

How to cmaple print int(f(x),x) after append to text command instead of /

                                             |

How can I make Maple preserve old outcomes?

I can display an animated plot of the Normal distribution pdf and was wondering if anybody has a routine to animate the Normal cdf. If this is so - is it possible to generalise for any continuous statistical distribution?

Thanks for reading! 

Dear Community!

I'm struggleing with a problem long since. I highly appreciate any help with this theme.

The problem is the following:

- I have a set of points, what comes from a numerical solution of a complex, but periodical function. Therefore I have a set of points (X;Y). The points are doesn't matter, but in my case, it look like this:

I want to use the Fourier-method to approximate this points with a function.

The best result I could get is this:

But it is not acceptable due to the high inaccuracy at the starting, and finishing points (there is a diagram inaccuracy %):

I'm feel like, I'm doing something wrong. Unfortunately, I had no time to look deep into the math here. Can somebody tell me, that how can I get a better result, using this method?

Thank you very much for the help in advance.

Best regards

Dávid

Hi guys,

I'm trying to apply boundary condition on the seris i generated using Adomian decomposition. Could some one please check it for me?
 

Download second_problem.mw

Hi everyone,

Please, I really need your expertise advice(s) on what i am not doing right in the code attached below. I was actually writing the code on Multi-step DTM, but, instead of continuing from the last point, it starting all over again. Please, your expertise will save a soul here. Regard

Thanks

MsDTM3.mw