I have sets of names with numeric parts like X:={x12, x1, x3, x15, x2, x9}.

How can I sort X based on the numeric part so Xsort:={x1, x2, x3, x9, x12, x15}.

Thank you.

Trying to use the jets package from http://jets.math.slu.cz/ (file jets.s) for calculus of variations.

Test example -

coordinates([x,y],[f],3);
parameter(p);
g:=(f_x)^2+(f_y)^2-2*p*f;
vg:=variation(g,f);                  # first variation of g (from line 4354 of jets.s)

Now I want to find g from vg - inverse problem. The code for this starts on line 4367 of jets.s, command is lagrangian.

Tried -

Lg:=lagrangian(vg,f);
Lg:=lagrangian([vg,f]);
Lg:=lagrangian(vg,[f]);
Lg:=lagrangian(f=vg);

but all give syntax errors.

The authors of the package don't respond to query and the manual does't help.

Does anyone know what the correct syntax of lagrangian is. Or is there a better package to use for inverse problem.

Thank you.

I want to change the output of the mtaylor command, eg for 3 series terms

ftaylor:=mtaylor(f(x,y),[x,y],3);

gives

ftaylor := f(0,0)+D[1](f)(0,0)*x+D[2](f)(0,0)*y+1/2*D[1,1](f)(0,0)*x^2+D[1,2](f)(0,0)*x*y+1/2*D[2,2](f)(0,0)*y^2;

I want the output as ftaylor := g(x,y) + D[1](g(x,y)) *x + D[2](g(x,y))*y + .... etc for all the terms in ftaylor.

Later on g(x,y) can be defined and terms like D[1](g(x,y)) evaluated.

I tried

subs(f(0,0)=g(x,y),ftaylor);

but it doesn't work. Using eval doesn't work either.

Is there a way to do this?

Thank you in advance.