MaplePrimes Questions

How do I calculate run of 0s or run of 1s using maple?


Q[1] := (e^(-n*T*s)-e^(-(n+1)*T*s)+(-exp(-Z[1]*n*T)*(s-Z[1])*exp(-n*T*(s-Z[1]))+exp(-Z[2]*n*T)*(s-Z[2])*exp(-n*T*(s-Z[2])))/(Z[1]-Z[2])+2*exp(-n*T*s)*(-1+Heaviside(-n*T)))/c+Z[1]*Z[2]*exp(-n*T*(s-Z[1]))/((s-Z[1])*exp(Z[1]*n*T)*(Z[1]-Z[2])*c)-Z[2]*Z[1]*exp(-n*T*(s-Z[2]))/((s-Z[2])*exp(Z[2]*n*T)*(Z[1]-Z[2])*c); 1; Q[2] := ((s-Z[1])*exp(-n*T*(s-Z[1]))*exp(-Z[1]*n*T)/(Z[2]-Z[1])+exp(-n*T*s)*(-1+Heaviside(-n*T)))/c+Z[1]*Z[2]*exp(-n*T*(s-Z[1]))/((s-Z[1])*exp(Z[1]*n*T)*(Z[1]-Z[2])*c)+(exp(-Z[2]*n*T)*(s-Z[2])*exp(-n*T*(s-Z[2]))/(Z[1]-Z[2])+exp(-n*T*s)*(-1+Heaviside(-n*T)))/c-Z[2]*Z[1]*exp(-n*T*(s-Z[2]))/((s-Z[2])*exp(Z[2]*n*T)*(Z[1]-Z[2])*c)+(e^(-n*T*s)-e^(-(n+1)*T*s))/c

((s-Z[1])*exp(-n*T*(s-Z[1]))*exp(-Z[1]*n*T)/(Z[2]-Z[1])+exp(-n*T*s)*(-1+Heaviside(-n*T)))/c+Z[1]*Z[2]*exp(-n*T*(s-Z[1]))/((s-Z[1])*exp(Z[1]*n*T)*(Z[1]-Z[2])*c)+(exp(-Z[2]*n*T)*(s-Z[2])*exp(-n*T*(s-Z[2]))/(Z[1]-Z[2])+exp(-n*T*s)*(-1+Heaviside(-n*T)))/c-Z[2]*Z[1]*exp(-n*T*(s-Z[2]))/((s-Z[2])*exp(Z[2]*n*T)*(Z[1]-Z[2])*c)+(e^(-n*T*s)-e^(-(n+1)*T*s))/c

(1)

Q[1] = Q[2]"(->)"true"(->)"true"(->)"true"(->)"true

Q[2] = ((s-Z[1])*exp(-n*T*(s-Z[1]))*exp(-Z[1]*n*T)/(Z[2]-Z[1])+exp(-n*T*s)*(-1+Heaviside(-n*T)))/c+Z[1]*Z[2]*exp(-n*T*s)/((s-Z[1])*(Z[1]-Z[2])*c)+(exp(-Z[2]*n*T)*(s-Z[2])*exp(-n*T*(s-Z[2]))/(Z[1]-Z[2])+exp(-n*T*s)*(-1+Heaviside(-n*T)))/c-Z[2]*Z[1]*exp(-n*T*s)/((s-Z[2])*(Z[1]-Z[2])*c)+(e^(-n*T*s)-e^(-T*s)*e^(-n*T*s))/c
"(->)"true
Q[2] = ((s-Z[1])*exp(-n*T*(s-Z[1]))*exp(-Z[1]*n*T)/(Z[2]-Z[1])+exp(-n*T*s)*(-1+Heaviside(-n*T)))/c+(exp(-Z[2]*n*T)*(s-Z[2])*exp(-n*T*(s-Z[2]))/(Z[1]-Z[2])+exp(-n*T*s)*(-1+Heaviside(-n*T)))/c+Z[1]*Z[2]*exp(-n*T*s)/((s-Z[1])*(Z[1]-Z[2])*c)-Z[2]*Z[1]*exp(-n*T*s)/((s-Z[2])*(Z[1]-Z[2])*c)+(e^(-n*T*s)-e^(-T*s)*e^(-n*T*s))/c
"(->)"true"(->)"true``

((s-Z[1])*exp(-n*T*(s-Z[1]))*exp(-Z[1]*n*T)/(Z[2]-Z[1])+exp(-n*T*s)*(-1+Heaviside(-n*T)))/c+(exp(-Z[2]*n*T)*(s-Z[2])*exp(-n*T*(s-Z[2]))/(Z[1]-Z[2])+exp(-n*T*s)*(-1+Heaviside(-n*T)))/c = (2*Heaviside(-n*T)-1)*exp(-n*T*s)/c
"(->)"true"(->)"true"(->)"true

Q[2] = (2*Heaviside(-n*T)-1)*exp(-n*T*s)/c+Z[1]*Z[2]*exp(-n*T*s)/((s-Z[1])*(Z[1]-Z[2])*c)-Z[2]*Z[1]*exp(-n*T*s)/((s-Z[2])*(Z[1]-Z[2])*c)+(e^(-n*T*s)-e^(-T*s)*e^(-n*T*s))/c
"(->)"false"(->)"false"(->)"false

``


Download inexplicable.mw

Can anyone explain the false return on the last line?  MAPLE seems to recognize the simplified expression on the next to last line, but when substituted into the expression for Q2 MAPLE does not seem to recognize the simplification.

How can I plot a volume (many surfaces) in X,Y, Z axis where X,Y,Z are functions in 4 variables (a,b,c,d), and the domain for the 4 variables are 

-90<=a>=90, -10<=b>=10, -12<=c>=12, -90<=d>=0,

 

X := proc (a, b, c, d) options operator, arrow; 324.*cos(b)*sin(c)*cos(d)+324.*sin(b)*sin(d)+323.5*cos(b)*sin(c) end procX := proc (a, b, c, d) options operator, arrow; 324.*cos(b)*sin(c)*cos(d)+324.*sin(b)*sin(d)+323.5*cos(b)*sin(c) end proc

Y := proc (a, b, c, d) options operator, arrow; (324*1.*sin(a)*sin(b)*sin(c)+324*1.*cos(a)*cos(c))*cos(d)+(-1)*324.*sin(a)*cos(b)*sin(d)+323.5*sin(a)*sin(b)*sin(c)+323.5*cos(a)*cos(c)-100 end proc

Y := proc (a, b, c, d) options operator, arrow; (324*1.*sin(a)*sin(b)*sin(c)+324*1.*cos(a)*cos(c))*cos(d)+(-1)*324.*sin(a)*cos(b)*sin(d)+323.5*sin(a)*sin(b)*sin(c)+323.5*cos(a)*cos(c)-100 end proc

Z := proc (a, b, c, d) options operator, arrow; (324*cos(a)*sin(b)*sin(c)-324*sin(a)*cos(c))*cos(d)-324*cos(a)*cos(b)*sin(d)+323.5*cos(a)*sin(b)*sin(c)+(-1)*323.5*sin(a)*cos(c)+150 end proc

Z := proc (a, b, c, d) options operator, arrow; (324*cos(a)*sin(b)*sin(c)-324*sin(a)*cos(c))*cos(d)-324*cos(a)*cos(b)*sin(d)+323.5*cos(a)*sin(b)*sin(c)+(-1)*323.5*sin(a)*cos(c)+150 end proc

how to compute the ideal mapping from ideal A to ideal B

I want to cancell the common factors in Maple using some simple code. I use both factor and normal command to cancell two factors which are equal upto certain places of decimal and then after that the factors differ in their digits.

e.g.,

1.520790243780735576142660664279654952049482710775236871695628058010

1.520790243780735576142660664279654952049482702492695064053058915917

which match upto more than 40 digits.

 Using fnormal(f(x),40), both factors are cancelled automatically.

 

Is there any such code to canell the above factors as well as the following one,

e.g.,

0 and

4.9482702492695064053058915917266066427965495204 10^(-101)

Let Poly2 denote the vector space of polynomials

(with real coefficients) of degree less than 3.

Poly2 = {a1t^2+ a2 t+ a3 |a1; a2; a3 €R}

You may assume that {1,t; t^2}is a basis for Poly2.

(1) Show that L1 = {t^2 + 1; t-2 ; t + 3}and L2 = {2 t^2 + t; t^2 + 3; t}

are bases for Poly2.

(2) Let = 8t^2- 4+ 6 and = 7t^2- t + 9. Find the coordinates for

and with respect to the basis L1 and with respect to the basis L2

(3) find the coordinate change matrix P from the basis L1 to the basis L2.find P^-1

Just I answer part (1) can you help me to answer 2 and 3 

I'm working on a complex problem in Composite Materials. I've gotten to a near-result 6x6 matrix, with several cells containing polynomial denominators. I have an equation for simplifying these, which boils the polynomials down to a single variable, but I can't seem to get it to substitute in. Can anyone help me solve this? The problem is also time-sensitive.

 

I would like to plot this :

yz(x+y+2z)^3-8=0 such as x^2+x(y+2z)-1<0 and y^2+y(x+2z)-3<0

The difficulty is that bounds depend on other variables. I was trying to use Maple implicitplot3d function, but I don't see a way.

Is it possible to use Maple for such plot ? If so, how ?

 

Thank you for your answer

I am considering to write a wrapper for plot and related commands (could redefine the commands or introduce a new name) which facilitates export of Maple plots to postscript. The command should interpret some of the options and remove them from the options sequence before submitting the remaining ones to the original plot().
E.g. it should recognize a title="TITLE" parameter and process TITLE (e.g. write it to a specific file). Similarly I would want to be able to pass additional parameters, e.g. filename="FILE" in order to specify how the output file name should be set. Is this a sensible approach. How can I realize this detailed option 'parsing' in Maple?

Hi all,

I created a package and saved it already. I checked and found it at D:/CaoHoc/LuanVan/Code/workspace/HINHHOC9.m but when I load it using with command, it has error.

[.............Code to create package..............]

>url := currentdir();
                   url := "D:\CaoHoc\LuanVan\Code\workspace"

>
>save HINHHOC9, cat(url, "/HINHHOC9.m");
>
>libname := libname, url;
libname := "C:\Program Files\Maple 18\lib",

    "D:\CaoHoc\LuanVan\Code\workspace\LuanVan\binaries",

    "C:\Program Files\Maple 18\lib", "D:\CaoHoc\LuanVan\Code\workspace"

>
>         with(HINHHOC9);
Error, (in with) user level initialization for package `HINHHOC9' failed: file
or directory does not exist

Please help me to identify the problem. Thanks

how to calculate hlibert series as in maple with Gröbner Bases

would like to know the algorithm and try in another programming language such as F#

i find the algorithm in book Singular introduction to commutative algebra

page 320 and 322 

1. is it equal to the hilbert series function in maple?

eq1a := Homogenize(eq1, h);
eq2a := Homogenize(eq2, h);
eq3a := Homogenize(eq3, h);
T3:=lexdeg([a,b,c,h]);
GB := Basis([eq1a,eq2a,eq3a], T3); #a

MonomialHilbertPoincare(LeadingMonomial(GB[1],T3), LeadingMonomial(GB[2],T3), LeadingMonomial(GB[3],T3));

F:=[LeadingMonomial(GB[1],T3), LeadingMonomial(GB[2],T3), LeadingMonomial(GB[3],T3)];
InterReduce(F, ???);

 2. what is the maple function for degree reverse lex ordering ?

eq1a := Homogenize(eq1, h);
eq2a := Homogenize(eq2, h);
eq3a := Homogenize(eq3, h);
David Cox using Algebraic Geometry page 82 use resultant to eliminate variable h
eq1b := eq1a - x;
eq2b := eq2a - y;
eq3b := eq3a - z;
T2:=lexdeg([a,b,c],[x,y,z]);
GB := Basis([eq1b,eq2b,eq3b], T2);
r1 := resultant(eq1b, eq2b, h);
r2 := resultant(eq1b, eq3b, h);
r1 = r2

page 82 teach how to eliminate, after do question 2, discover r1 and r2 are the same.

how to eliminate the variable h with resultant after homogenize ideal with variable h

Why can't maple 15 solve this eqn. [n= 10]

solve(ithprime(n)=29,n);

 

 Do Hilbert series function classify all or only some type or some form of ideals?

 

Let Poly2 denote the vector space of polynomials

(with real coefficients) of degree less than 3.

Poly2 = {a1t^2+ a2 t+ a3 |a1; a2; a3 €R}

You may assume that {1,t; t^2}is a basis for Poly2.

(1) Show that L1 = {t^2 + 1; t-2 ; t + 3}and L2 = {2 t^2 + t; t^2 + 3; t}

are bases for Poly2.

(2) Let v = 8t^2- 4t + 6 and w = 7t^2- t + 9. Find the coordinates for

v and w with respect to the basis L1 and with respect to the basis L2

First 1384 1385 1386 1387 1388 1389 1390 Last Page 1386 of 2434