Hi all,

I'm trying to plot the regions of a cube that are defined by an inequality. At the moment I am able to plot only the surface defined by the corresponding equality. 

That's the code: 

with(plots)

implicitplot3d(x^2*y = z, z = 0 .. 1, y = 0 .. 1, x = 0 .. 1);

Now I would like to plot the two regions defined by the corresponding inequality. 

Have you got any idea on how to do this?

Bests

Manuele

Hi My main problem is that the new installer is ridiculously small that I cannot even press the buttons acuratley. I am running windows 8.1 on a Surface Pro 2. I had no problems installing maple 17.

I also ran into errors using the Bitrock installer 3 times, I am not even sure if it is installed correctly now.

is there anyway to get maple 18 without using the Bitrock installer?

Hello guys ...

I used a numerically method to solve couple differential equation that it has some boundary conditions. My problem is that some range of answers has 50% error . Do you know things for improving our answers in maple ?

my problem is :

a*Φ''''(x)+b*Φ''(x)+c*Φ(x)+d*Ψ''(x)+e*Ψ(x):=0

d*Φ''(x)+e*Φ(x)+j*Ψ''(x)+h*Ψ(x):=0

suggestion method by preben Alsholm:

a,b,c,d,e,j,h are constants.suppose some numbers for these constants . I used this code:


VR22:=0.1178*diff(phi(x),x,x,x,x)-0.2167*diff(phi(x),x,x)+0.0156*diff(psi(x),x,x)+0.2852*phi(x)+0.0804*psi(x);
VS22:=0.3668*diff(psi(x),x,x)-0.0156*diff(phi(x),x,x)-0.8043*psi(x)-0.80400*phi(x);
bok:=evalf(dsolve({VR22=0,VS22=0}));

PHI,PSI:=op(subs(bok,[phi(x),psi(x)]));
Eqs:={eval(PHI,x=1.366)=1,eval(diff(PHI,x),x=1.366)=0,eval(PHI,x=-1.366)=1,eval(diff(PHI,x),x=-1.366)=0,
eval(PSI,x=1.366)=1,eval(PSI,x=1.366)=1};
C:=fsolve(Eqs,indets(%,name));
eval(bok,C);
SOL:=fnormal(evalc(%));


I used digits for my code at the first of writting.

please help me ... what should i do?

X belongto A, eigenvector(X) = 0

from this statement , 

using linearalgebra package eigenvectors function

the eigenvector matrix [3][1],[3][2],[3][3] are 1 , contradict 1=0

so, need to find another kind of eigenvector in terms of algebra 

using original basic calculation solve, however got error

m := Matrix([[a1,a2,a3],[a4,a5,a6],[a7,a8,a9]]);
eigenvector1 := Eigenvectors(m);
solve(
[eigenvector1[2][1][1]=0, eigenvector1[2][1][2]=0, eigenvector1[2][1][3]=0,
eigenvector1[2][2][1]=0, eigenvector1[2][2][2]=0, eigenvector1[2][2][3]=0,
eigenvector1[2][3][1]=0, eigenvector1[2][3][2]=0, eigenvector1[2][3][3]=0]
);

solve(
[eigenvector1[2][1][1]=0, eigenvector1[2][1][2]=0, eigenvector1[2][1][3]=0,
eigenvector1[2][2][1]=0, eigenvector1[2][2][2]=0, eigenvector1[2][2][3]=0]
);

eigenvalue1 :=
(1/6)*(36*a7*a1*a3+108*a7*a2*a6+108*a8*a4*a3+36*a8*a5*a6

...

eigenvalue2 :=
-(1/12)*(36*a7*a1*a3+108*a7*a2*a6+108*a8*a4*a3+36*a8*a5*a6

...

eigenvalue3 :=
-(1/12)*(36*a7*a1*a3+108*a7*a2*a6+

...

solve(MatrixMatrixMultiply(Matrix([[a1,a2,a3],[a4,a5,a6],[a7,a8,a9]])-Matrix([[eigenvalue1,0,0],[0,eigenvalue1,0],[0,0,eigenvalue1]]), Matrix([[x],[y],[z]])),[x,y,z]);
solve(MatrixMatrixMultiply(Matrix([[a1,a2,a3],[a4,a5,a6],[a7,a8,a9]])-Matrix([[eigenvalue2,0,0],[0,eigenvalue2,0],[0,0,eigenvalue2]]), Matrix([[x],[y],[z]])),[x,y,z]);
solve(MatrixMatrixMultiply(Matrix([[a1,a2,a3],[a4,a5,a6],[a7,a8,a9]])-Matrix([[eigenvalue3,0,0],[0,eigenvalue3,0],[0,0,eigenvalue3]]), Matrix([[x],[y],[z]])),[x,y,z]);

got error when using solve

> solve(MatrixMatrixMultiply(Matrix([[a1, a2, a3], [a4, a5, a6], [a7, a8, a9]])-Matrix([[eigenvalue1, 0, 0], [0, eigenvalue1, 0], [0, 0, eigenvalue1]]), Matrix([[x], [y], [z]])), [x, y, z]);
Error, invalid input: solve expects its 1st argument, eqs, to be of type {`and`, `not`, `or`, algebraic, relation(algebraic), ({set, list})({`and`, `not`, `or`, algebraic, relation(algebraic)})}, but received Matrix(3, 1, {(1, 1) = ((2/3)*a1-(1/6)*(36*a7*a1*a3+108*a7*a2*a6+108*a8*a4*a3+36*a8*a6*a5-72*a7*a3*a5-72*a8*a6*a1-72*a9*a4*a2+48*a9*a5*a1-12*a9*a1^2-12*a5*a1^2+8*a1^3-12*a9^2*a1-12*a5^2*a1-12*a9^2*a5-12*a9*a5^2+36*a8*a6*a9+36*a7*a3*a9+36*a4*a2*a1+36*a4*a2*a5+8*a9^3+8*a5^3+12*(54*a7*a2^2*a6*a4*a1+114*a8*a6*a9*a1*a4*a2+6*a8*a6*a9*a1*a7*a3+54*a8*a4*a3^2*a7*a9-60*a9*a1^2*a8*a6*a5-60*a8*a6*a7*a3*a5^2-60*a8*a6*a4*a2*a9^2-24*a9*a1*a4^2*a2^2+6*...

MatrixOperation := module() option package;  export `+`, LinearAlgebra;
    `+` := proc(a::float, b::float) option overload;
 :-`+`(map(x->x^2,a),map(x->x^2,b));
    end proc;
end module;

with(MatrixOperation);
with(LinearAlgebra):
m := Matrix([[1,2],[3,4]]);
L := MatrixMatrixMultiply(m,m);

1 2  1 2
3 4  3 4

1*1+2*3 = 1 + 6 = 1 after overload + with a+b-a*b
1*2+2*4 = 2 + 8 = -6

L should be Matrix([[1, -6],[....]])

 http://en.wikibooks.org/wiki/Linear_Algebra/

Representing_Linear_Maps_with_Matrices

how to calculate the first step

(2,0) -> (1,1,1) and (1,4) -> (1,2,0)

how to use maple command to get (1,1,1) and (1,2,0)

how to use maple command to calculate rep(h)

to get (0,-1/2,1) and (1,-1,0)

http://en.wikibooks.org/wiki/Linear_Algebra

/Representing_Linear_Maps_with_Matrices

hello evreybody i have these Error :

Download fffffffff_INFINI.mw

please help 
thank you !

Hi Everybody,

I installed the Maple Toolbox for matlab and tried to run my old code.  For some reason, maple('restart') will cause Matlab to lock up.  Any ideas?

Windows 7 64-bit.  Matlab 2012b 32-bit, Maple 18 32-bit

Hi Mapleprimers,

I'm using CodeGeneration to convert a procedure I obtained with unapply() into a Matlab function.  I'm having problems getting the outputting function to run correctly in Matlab.  I'm going to dynamically generating equations, so directly editing the Matlab code won't work here.  I'm having problems getting any output in Matlab.  Here is the code I'm working on:  Series_addGear_codegen.mw

Ideally I would like output in a matrix.  I've tried putting the unapplied procedure in another procedure, but the CodeGeneration doesn't work.

This is the maple output from the function:

unapp(1,2,3,4,5);
{BAT_A = -2.267032891, BAT_V = 271, EM2_A = .4615464218, EM2_P = 125.0790803, EM2_T = -1, EM2_V = 271, EM2_W = 2, GBa_T = 12, GBa_W = 5/3, GBb_T = -4, GBb_W = 5, GEN_A = 1.805486469, GEN_P = 489.2868330, GEN_T = -12, GEN_V = 271, GEN_W = 5/3, ICE_mdot_g = 20}

Since I'll know the order of variables, I want the Matlab function to output:

[-2.267032891,  271,  .4615464218,  125.0790803,  -1, ...]

This is the output after putting the Matlab function in Matlab:

>> unapp(1, 1, 1, 1, 1)
Undefined function or variable 'BAT_A'.

Error in unapp (line 39)
unappreturn = unique([BAT_A == -t35 / 0.271e3 - t73 / 0.271e3 BAT_V == 271 EM2_A == t35 /
0.271e3 EM2_P == t35 EM2_T == -FD_T EM2_V == 271 EM2_W == FD_W GBa_T == t44 GBa_W == t41 GBb_T
== -ICE_T GB

This is the Matlab code that is generated by Maple 18:

MCode:=CodeGeneration[Matlab](unapp, declare = [FD_T::float, FD_W::float, GB_R::float, ICE_T::float, ICE_W::float],defaulttype=float,optimize,defaulttype = numeric);
Warning, could not preprocess. Found `abs` or similar in the 'While'/'For' conditions.
Warning, procedure/module options ignored
function unappreturn = unapp(FD_T, FD_W, GB_R, ICE_T, ICE_W)
t2 = 0.1e1 * FD_T * FD_W;
t5 = abs(FD_W);
t7 = abs(FD_T);
t9 = t5 ^ 2;
t13 = t7 ^ 2;
t15 = t9 * t5;
t21 = t13 * t7;
t23 = t9 ^ 2;
t31 = t13 ^ 2;
t33 = 0.1483000000e3 - 0.4267000000e1 * t5 - 0.1277000000e2 * t7 + 0.3640000000e-1 * t9 - 0.1160000000e1 * t5 * t7 + 0.2580000000e0 * t13 - 0.1181000000e-3 * t15 + 0.5994000000e-3 * t9 * t7 - 0.1171000000e-3 * t5 * t13 - 0.1739000000e-2 * t21 + 0.1245000000e-6 * t23 - 0.1200000000e-5 * t15 * t7 + 0.1584000000e-5 * t9 * t13 - 0.4383000000e-6 * t5 * t21 + 0.2947000000e-5 * t31;
if (-t2 == 0.0e0)
t35 = 0.0e0;
elseif (-t2 < 0.0e0)
t35 = t33;
else
t35 = -t33;
end
t36 = ICE_T * ICE_W;
t37 = 0.1e1 * t36;
t41 = ICE_W / GB_R;
t42 = abs(t41);
t44 = ICE_T * GB_R;
t45 = abs(t44);
t47 = t42 ^ 2;
t51 = t45 ^ 2;
t53 = t47 * t42;
t59 = t51 * t45;
t61 = t47 ^ 2;
t69 = t51 ^ 2;
t71 = 0.5280000000e-11 - 0.3849000000e-13 * t42 + 0.7190000000e2 * t45 + 0.1168000000e-15 * t47 - 0.1296000000e1 * t42 * t45 - 0.2489000000e1 * t51 - 0.1451000000e-18 * t53 - 0.1326000000e-3 * t47 * t45 + 0.8141000000e-2 * t42 * t51 + 0.4539000000e-2 * t59 + 0.6325000000e-22 * t61 + 0.2091000000e-6 * t53 * t45 - 0.3455000000e-5 * t47 * t51 - 0.2499000000e-4 * t42 * t59 + 0.5321000000e-4 * t69;
if (-t37 == 0.0e0)
t73 = 0.0e0;
elseif (-t37 < 0.0e0)
t73 = t71;
else
t73 = -t71;
end
unappreturn = unique([BAT_A == -t35 / 0.271e3 - t73 / 0.271e3 BAT_V == 271 EM2_A == t35 / 0.271e3 EM2_P == t35 EM2_T == -FD_T EM2_V == 271 EM2_W == FD_W GBa_T == t44 GBa_W == t41 GBb_T == -ICE_T GBb_W == ICE_W GEN_A == t73 / 0.271e3 GEN_P == t73 GEN_T == -t44 GEN_V == 271 GEN_W == t41 ICE_mdot_g == t36]);

Consider the following code:

with(LinearAlgebra):
with(Physics):
Setup(anticommutativeprefix = psi):
psiFermi := Vector(2,symbol = psi):
psiBose  := Vector(2,symbol = phi):
A := Matrix([[0,1],[1,0]]):
Transpose(psiFermi) . A;
Transpose(psiBose ) . A;

It produces the following output:

Why is the first line, for anticommuting components, not evaluated to the same form as the second line, for commuting components? The actual choice of the matrix A seems immaterial; the odd behaviour is present even if A is chosen to be the identity matrix!

In comparison, the 'contracted' (scalarly) expressions

Transpose(psiFermi) . A . psiFermi,
Transpose(psiBose ) . A . psiBose;

produce the following completely sensible output:

I am not able to simplify my equation, any help would be appreciated ! I want the  term to be eliminated

Download 1.mw

Hello,

I'm quite new to Maple and I have a serious problem when I'm trying to solve this system of equations for (a,b,p,v,u,g):

1) alpha[1]= v a u (1- b)
2) alpha[2]= v a ub
3) alpha[3]= v (1-a)=v-va
4) alpha[11]= 1/(2) auv (-1+b) (-b u p+b a u p+b g+p u-u p a) 
5) alpha[22]= 1/(2) a u v b (-b u p+b a u p+b g-g)
6) alpha[33]= 1/(2)(a- 1)apv
7) alpha[12]= -a u v b(pu-pua-g - bup +baup +bg)
8) alpha[13]= (a-1) (b-1)ap vu
9) alpha[23]= (1-a) a p v ub

I tried this command:

solve({v*(1-a) = alpha[3], a*u*v*b = alpha[2], v*a*u*(1-b) = alpha[1], (1/2*(-1+a))*a*p*v = alpha[33], (-1+a)*(-1+b)*a*p*v*u = alpha[13], (1-a)*a*p*v*u*b = alpha[23], (1/2)*a*u*v*b*(-b*u*p+b*a*u*p+b*g-g) = alpha[22], -a*u*v*b*(p*u-u*p*a-g-b*u*p+b*a*u*p+b*g) = alpha[12], (1/2)*a*u*v*(-1+b)*(-b*u*p+b*a*u*p+b*g+p*u-u*p*a) = alpha[11]}, {a, b, g, p, u, v}, 'parametric' = 'full', 'parameters' = {alpha[1], alpha[2], alpha[3], alpha[11], alpha[12], alpha[13], alpha[22], alpha[23], alpha[33]})

 but the Maple output is [ ]. I can find a solution manually but I don't understand why I cannot do it with Maple. It's very important that I find a solution as I have a much more complicated system to solve in a similar manner.

Thank you very much for your help!!

Elena

I have a long and complicated expression (say Eq1). I have to solve the two equations (obtained by taking real and imaginary parts of Eq1) for two unknowns. Normally it take around a day on my i5 (3.1 GHz, 4 cores), 8gb ram desktop. I looked into posts related to parallel programming but couldn't get much to reduce the computation time. Is there a way to reduce the computation time ?
Many thanks in advance.

Here's an example compound inequality I'm working on.

Working it out manually.... 

Compound Inequality
4477.25 <= 4477.25+.25*(t-32450) <= 16042.25;

Distribute the coefficient
4477.25 <= 4477.25+.25*t - 8112.50 <= 16042.25;

Combine like terms
4477.25 <= -3635.25+.25*t <= 16042.25;

Add 3635.25 to all sides
8112.50 <= .25*t <= 19677.50;

Divide all sides by .25
32450 <= t <= 78710;

How can I ask Maple to simplify this compound inequality? Obviously this is not the correct syntax, It seems Maple doesn't understand what I want it to do.

4477.25 <= 4477.25 + .25 * (t-32450) <= 16042.25;

                       0.00 <= 0.25 t - 8112.50 and 0.25 t <= 19677.50                (112)

Also is there a way to ask Maple to only perform one step? In the above example, is it possible to ask Maple to "Distribute the .25", then show the result, next ask it to combine like terms, etc?