Maple 15 Questions and Posts Maple 15 Questions and Posts Feed

These are Posts and Questions associated with the product, Maple 15

cerebrrin

 

One thing to remember is that a lot of children, even thought have been not born or aware when these attacks happened, are no stranger to your cause and affect we have now using on an authentic bases. When family has ever traveled by airplane they may already be conscious of the extream meassures airports security take to ensure safety. If you are a military family they may all ready be use to frequent hightened security levels on base, bomb sniffing dogs, and limiting non-military personnel on assiette. Military kids may already know of family or friends who were personally involved the intrusions. They may have knowledge of online marketing missions inside of Middle East and why dad recently been sent over there (or soon). Here just what you in a position to educate your children about 9/11.

 

 

http://rhinorx90eveningblog.com/cerebrrin/

     Parallel curves on surfaces. The distance between the points of the curves is measured along the curves of intersection of the surface and perpendicular planes.
     (According to tradition, it also does not make sense.)

equidistant_curve_surface_MP.mw

 









 

 

I'm running into a very simple problem with the way that Maple integrates Heaviside functions. Naively, it should act like a step function, but it is not integrating properly. See the attached document.

int(int(Heaviside(-x^2-y^2+1), x = -1 .. 1), y = -1 .. 1)

0

(1)

evalf(Int(Heaviside(-x^2-y^2+1), [x = -1 .. 1, y = -1 .. 1]))

3.141592654

(2)

int(piecewise(-x^2-y^2+1 > 0, 1, 0), [x = -1 .. 1, y = -1 .. 1])

Pi

(3)

``


Note that the symbolic integration of the Heaviside function (defined to be 1 inside the unit circle and 0 outside) gives zero, whereas it should clearly give the area of the unit circle, which the numerical integration does. I even checked that the (suposedly equivalent) piecewise definition symbolically evaluates to the area, and it, too, gets the right answer.

Anyone have any clue as to why the symbolic integration of this Heaviside function is so wrong? My understanding is that if we do the integral as two nested 1D integrals, the returned function (as a function of y) is zero everywhere except at y=0, but that result cannot be right either.

Thoughts?

 

Download Heaviside-error.mw

      General description of the method of solving underdetermined systems of equations. As a particular application of the idea proposed a universal method of calculation for all kinds of link (lever) mechanisms. With the description and examples.
      The method can be used for powerful CAD linkages.

Description: Calculation_method_of_linkages.pdf

Attachment:
figure_1.mw
figure_2.mw

Or all in one
Calculation_method_of_linkages_(with_attach.).pdf


        Some examples of a much larger number calculated by the proposed method. Examples gathered here not to look for them on the forum and opportunity to demonstrate the method.  Among the examples, I think, there are very complicated.

https://vk.com/doc242471809_408704758
https://vk.com/doc242471809_408704572
https://vk.com/doc242471809_376439263
https://vk.com/doc242471809_402619761
https://vk.com/doc242471809_402610228
https://vk.com/doc242471809_401188803
https://vk.com/doc242471809_400465891
https://vk.com/doc242471809_400711315
https://vk.com/doc242471809_387358164
https://vk.com/doc242471809_380837279
https://vk.com/doc242471809_379935473
https://vk.com/doc242471809_380217387
https://vk.com/doc242471809_363266817
https://vk.com/doc242471809_353980472
https://vk.com/doc242471809_375452868
https://vk.com/doc242471809_353988163 
https://vk.com/doc242471809_353986884 
https://vk.com/doc242471809_353987119
https://vk.com/doc242471809_324249241
https://vk.com/doc242471809_324102889
https://vk.com/doc242471809_322219275
https://vk.com/doc242471809_437298137
https://vk.com/doc242471809_437308238
https://vk.com/doc242471809_437308241
https://vk.com/doc242471809_437308243
https://vk.com/doc242471809_437308245
https://vk.com/doc242471809_437308246
https://vk.com/doc242471809_437401651
https://vk.com/doc242471809_437664558

 

 

Round := proc(x,n::integer:=1)
parse(sprintf(cat("%.",n,"f"),x));
end proc:

roundcoeffs1:=proc(p,x,n:=1) local t,c;
c:=map(Round, [coeffs(p,x,t)],n);
add(i, i = zip(`*`, c, [t]));
end:

ggg:=.9940413618*y^3-1.785839107*c*A*y^3-2.357517322*c*A*y^2+.375393240*c*y*B-.3575173222*c*A*y-.2082022533*c*B-0.1787591445e-1*y^2-0.1787591445e-1*y-0.5958638151e-2+.2141608926*c*A+.7917977467*c*B*y^3+2.375393240*c*B*y^2;

roundcoeffs1(ggg, [y^3, c*A*y^3, c*A*y^2, c*y*B, c*A*y, c*B, y^2, y, c*A, c*B*y^3, c*B*y^2], 4);


Error, (in sprintf) number expected for floating point format

with(Optimization):

theta := Complex(1,1);
Minimize(theta^3-3*(A*theta^2+B), {0 <= theta^3-3*(A*theta^2+B)}, assume = nonnegative)

Error, (in Optimization:-NLPSolve) complex value encountered

from mathematica,

 

n = 5;
CalabiYau[z_, k1_, k2_] := Module[{z1 = Exp[2Pi I k1/n]Cosh[z]^(2/n), z2 = Exp[2Pi I k2/n]Sinh[z]^(2/n)}, {Re[z1], Re[z2], Cos[alpha]Im[z1] + Sin[alpha]Im[z2]}];
Do[alpha = (0.25 + t)Pi; Show[Graphics3D[Table[ParametricPlot3D[CalabiYau[x + I y, k1, k2], {x, -1, 1}, {y, 0, Pi/2}, DisplayFunction -> Identity, Compiled ->False][[1]], {k1, 0, n - 1}, {k2, 0, n - 1}], PlotRange -> 1.5{{-1, 1}, {-1, 1}, {-1, 1}}, ViewPoint -> {1, 1, 0}]], {t, 0, 1, 0.1}];

 

n := 5;

z1 := exp(2*3.14*I*k1/n)*cosh(z)^(2/n);
z2 := exp(2*3.14*I*k2/n)*sinh(z)^(2/n);

alpha = (0.25 + t)Pi;

xx := Re(z1);
yy := Re(z2);
uu := cos(alpha)*Im(z1) + sin(alpha)*Im(z2);

 

where k1, k2, alpha are variables

print([xx,yy,uu]);

i find algcurve has implicitize

how to use this implicitize to find 3d surface?

is there any other method to find?

 

i searched groebner basis can do this, but in mathematica is different from maple example

i got this error in window 8 in surface 2  then follow this post and install again still error

https://www.maplesoft.com/support/faqs/detail.aspx?sid=139020

then follow

http://www.maplesoft.com/support/faqs/detail.aspx?sid=32607

then follow and install again same error

http://www.maplesoft.com/support/faqs/detail.aspx?sid=32631

and install again same eror

then i add option -f c:\Program File (x86)\MapleXX in cmd and then no error any more 

but no install succeed 

where it go, it still not install

then i try again, there is no room enough to install,  hard disk do not have enough space, then i go to c:\Windows\Temp, after deleted file in it, still not enough space

 

i find 

https://www.maplesoft.com/support/install/maple15_install.html

but template do not state how to activate later

how to write this template and how to clear the temp file created by previous failed cmd install method

Hi,

 

Let, fixed an integer i and 1<=j<=2^{i}-1, for each x and y in [0,1] let the following mapping

Then, with the above procedure we can obtained, for a fixed i, all the mappings for j=1,...,2^{i}-1

 

However, How can I to evalute the "components" of the above procedure? For instance, I can not to compute CreaF(2)[1](0.35,0.465) (i.e., the first function in the "vector" CreaF(2), in x=0.35, y=0.465). 

 

Thanks very much for your time.

 

Hi, I am new to maple, but I think that my question should be simple.

I have a matrix where each element is an expression. I want to compute the matrix for different constant and to save it without crushing the previous matrix. 

If the file that I joined, I have a first part where the constant are defined. In the second part the expression of the matrix is defined. Finally, I compute each matrix with different constant. Each results is called C_p0, C_s0, C_g0. When I called them back, only the last matrix computed remains.

I would like to be able to save each matrix to performed operation on them later.

Thank you. 

 

Forum_Question1.mw

Homogénéisation

 

restart; with(plots); with(DifferentialGeometry); with(LinearAlgebra); with(Physics)

  NULL

Paramètre des matériaux

 

p[p] := [34.68, 34.82]:
NULL

 

NULLNULL

Tenseurs Élémentaires

 

NULL

Tenseur de rigidité

 

V := 1/((1+upsilon[23])*(-2*upsilon[12]*upsilon[21]-upsilon[23]+1)); G[12] := E/(2*(1+upsilon[12])); C[11] := (-upsilon[23]^2+1)*V*E[1]; C[22] := (-upsilon[12]*upsilon[21]+1)*V*E[2]; C[12] := upsilon[21]*(1+upsilon[23])*V*E[2]; C[23] := (upsilon[12]*upsilon[21]+upsilon[23])*V*E[2]; C[44] := (1/2)*(-2*upsilon[12]*upsilon[21]-upsilon[23]+1)*V*E[2]; C[55] := E[6]; C[33] := C[22]; C[13] := C[12]; C[66] := C[55]; C[21] := C[12]; C[32] := C[23]; C[iso] := Matrix(6, 6, {(1, 1) = C[11], (1, 2) = C[12], (1, 3) = C[12], (1, 4) = 0, (1, 5) = 0, (1, 6) = 0, (2, 1) = C[21], (2, 2) = C[22], (2, 3) = C[23], (2, 4) = 0, (2, 5) = 0, (2, 6) = 0, (3, 1) = C[21], (3, 2) = C[32], (3, 3) = C[22], (3, 4) = 0, (3, 5) = 0, (3, 6) = 0, (4, 1) = 0, (4, 2) = 0, (4, 3) = 0, (4, 4) = C[44], (4, 5) = 0, (4, 6) = 0, (5, 1) = 0, (5, 2) = 0, (5, 3) = 0, (5, 4) = 0, (5, 5) = C[66], (5, 6) = 0, (6, 1) = 0, (6, 2) = 0, (6, 3) = 0, (6, 4) = 0, (6, 5) = 0, (6, 6) = C[66]})

Matrice de rigidité

 

upsilon[23] := upsilon[p]:

Matrix([[C[11], C[12], C[12], 0, 0, 0], [C[21], C[22], C[23], 0, 0, 0], [C[21], C[32], C[22], 0, 0, 0], [0, 0, 0, C[44], 0, 0], [0, 0, 0, 0, C[66], 0], [0, 0, 0, 0, 0, C[66]]])

(1.2.1.1.1)

upsilon[23] := upsilon[s]:

Matrix([[C[11], C[12], C[12], 0, 0, 0], [C[21], C[22], C[23], 0, 0, 0], [C[21], C[32], C[22], 0, 0, 0], [0, 0, 0, C[44], 0, 0], [0, 0, 0, 0, C[66], 0], [0, 0, 0, 0, 0, C[66]]])

(1.2.1.1.2)

upsilon[23] := upsilon[g]:

Matrix([[C[11], C[12], C[12], 0, 0, 0], [C[21], C[22], C[23], 0, 0, 0], [C[21], C[32], C[22], 0, 0, 0], [0, 0, 0, C[44], 0, 0], [0, 0, 0, 0, C[66], 0], [0, 0, 0, 0, 0, C[66]]])

(1.2.1.1.3)

``

C[p0];

Matrix([[C[11], C[12], C[12], 0, 0, 0], [C[21], C[22], C[23], 0, 0, 0], [C[21], C[32], C[22], 0, 0, 0], [0, 0, 0, C[44], 0, 0], [0, 0, 0, 0, C[66], 0], [0, 0, 0, 0, 0, C[66]]])

(1)

``

 

equidistant_curve_MP.mw  Equidistant curves to the curves on the surface. (Without any sense, but real.)







Hi,

For each i and 1<=j<=2^{i}-1, define

 

 

Then, with the below procedure we can to generate, for a given i, all functions for j=1,...,2^{i}-1

 

 

However, I need to evaluate  the functions in the vector returned by this procedure. Then, how can I to defined sich functions? For instance, with

 

 

the instruction h(0,0.4) returns

 

With h:=(x,y=-->[CreaF(3)][2](x,y) do not works. Somebody know how to "evaluate" these functios (i.e., the components of the vector) above generated by the procedure CreaF?

 

Thank you very much for your time.

 

 

 

 

I have the following expression (generated by some other procedure):

This does not have a taylor expansion in pV[6] in the general case because the square roots can become negative:

taylor(xpr,pV[6]);
Error, does not have a taylor expansion, try series()

But I can get an expansion by restrictig the range of pV[6]:

taylor(xpr,pV[6]) assuming -0.01<pV[6],pV[6]<0.01;

So far things are perfectly fine. But when I try mtaylor:

mtaylor(xpr,pV[6]) assuming -0.01<pV[6],pV[6]<0.01;
Error, (in assuming) when calling 'mtaylor'. Received: 'does not have a taylor expansion, try series()'

So the assumption seems to be ignored. I can work around this by expanding in pV[6] first, using taylor, and then expanding the result from that using mtaylor (I really also want the expansions in the other pV components; 6 in total although in this example some do not show up). I'll have to convince myself that this work-around gives the correct result but I think it does. However, I don't particularly like it.

I consider this a bug and am tempted to submit an SCR. But before I do that; is there anything obvious I am missing here?

Thanks,

M.D.

PS: This was done using Maple 15. I'll check newer versions later.

mtaylor_assuming.mw

     Example of the equidistant surface at a distance of 0.25 to the surface
x3
-0.1 * (sin (4 * x1) + sin (3 * x2 + x3) + sin (2 * x2)) = 0
Constructed on the basis of universal parameterization of surfaces.

equidistant_surface.mw 


When I was editing the head of the question (? instead of .), its body disappeared. Please, insert it again.

Regard,

Markiyan Hirnyk

1 2 3 4 5 6 7 Last Page 1 of 44