vv

12755 Reputation

19 Badges

8 years, 298 days

MaplePrimes Activity


These are answers submitted by vv

I have attached a worksheet.
An element of the algebra is represented as Vector(v1,...,vn,v0); also included a display procedure (optional)

e[i] &x e[j]  = 0 if not defined otherwise.

I have modified f because the indices in your definition exceed n; you also have provided two incompatible variants.

algebra-eiej-code-1.mw  (edit: updated)

 

 

You want Groebner[Basis], not LinearAlgebra[Basis] ! So, use:

G := Groebner[Basis](K, 'tord', degrevlex(r,u,v));

It is not very clear what are you trying to obtain.

Mathematically a polynomial map is a vector function (a,b,c) |--> (r,u,v) having polynomial components.

In Maple you may represent it by a list [r,u,v] or vector Vector([r,u,v]) provided that the indeterminates (i.e. variables) are a,b,c.

Or, you may consider the procedures:
p:=unapply( [r,u,v], a,b,c );
or
p:=unapply( Vector([r,u,v]), a,b,c );

The Maple Optimization package accepts Boolean variables for linear problems.
See ?Optimization and ?Optimization,Options help pages.

For nonlinear problems there is the DirectSearch package:

http://www.maplesoft.com/applications/view.aspx?SID=101333

 

 

Maple can find them, but must be asked politely! :-)

X := (x, k)-> int(cos(s^k), s = 0 .. x):

X(infinity,k) assuming k>1;

X(infinity,5/2); evalf(%);

     0.7178115070

Similarly for Y.

 

lagrange:=proc(f, g::list, vars::list, lambdas::list)
local L;
if nops(g)<>nops(lambdas) then error "lambdas and g must agree" fi;
L:=f + add(g[i]*lambdas[i],i=1..nops(g));
solve( {seq( diff(L,u),u=vars),op(g)},{op(vars),op(lambdas)})
end:

lagrange(a*b,[z-a-b],[a,b],[p]);

lagrange(a[2]*b[2], [z-a[2]-b[2]],[a[2],b[2]], [lambda]);

 

lagrange(x^2+y+z,[x-2*y-3*z, x+4*y+z-2 ],[x,y,z],[a,b]);

 

 Edit.
In your code use:
solve({eq1,eq2,eq3},{x,y,lambda});
and it will work.

 

f:=cos(t),sin(t);  # on [0,Pi]
T:=Pi; dx:=2;
fper:=cos(T*frac(t/T))+dx*floor(t/T),sin(T*frac(t/T));  # for t>=0 only

plot( [fper, t=0..4*T],scaling=constrained);




You can simply add p[i]+p[j] <> 0 etc in the system. E.g.

expand((x+y)*x*(x+y+5)):
solve({%, x+y<>0, x<>0}, {x,y});

 

The expresion

f:=3^(-(1/2)*n)*2^((1/6)*n)-2^((2/3)*n)*6^(-(1/2)*n);

should simplify to 0 even for a complex n !!

Now, convert(f,exp) should produce exp( ln(3) *...)...
but it gives

(IIRC Maple used to get this right.)

Without any assumption:

subs(n = ln(z), f):

simplify(%);

    0

f:=x->piecewise(x<1,x^2,(2-x)^2):

f_per:=x->f(x-2*floor(x/2)):

seq is faster but $ works in symbolic contexts e.g.

[x[i], y[i]] $ i=1..n;

diff( exp(a*x),  x$n);

where n is a symbol.

You may use

simplify(g,radical,symbolic);


Be aware that the validity depends on the domain of u in this case. 

Why don't you switch to 1D math and use simply:

diff(u, y, y, t);

 

 

It is a Maple evaluation peculiarity. In sum, the argument  N[i] is evaluated first (before the index i is defined).

You must use

sum('N[i]', i = 1..4);
    1

Or, better, use add, for which the evaluation rule is different.
add(N[i], i = 1..4);
   
1

restart;
p:=10:
dx_:=seq(x[i]=0..1,i=1..p-1):
sx_:=add(x[i],i=1..p-1):
H:=int( Heaviside(xp-sx_)*Heaviside(1-xp+sx_),  dx_ ):
int( exp(-xp^3)*H, xp=0..p,  numeric);
        0.000002546120390

# However the method is slow for higher accuracy and does not provide the exact result.
# It has the advantage that it works in other situations too.

First 105 106 107 108 109 110 111 Page 107 of 113