Maple 2016 Questions and Posts

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

hi.

please help me for remover this problem.

''''

-Float(infinity)*signum((5.*A3*A1-24.*A2^2)*A1/A2^2)''''

ReducedCantiler.mw
 

restart

f := -(2/3)*eta^3+(1/2)*eta^2+eta; -1; g := -eta^2+1; -1; h := -eta^2+1; 1; F := proc (eta) options operator, arrow; A1*f end proc; 1; G := proc (eta) options operator, arrow; A2*g end proc; 1; H := proc (eta) options operator, arrow; A3*h end proc

proc (eta) options operator, arrow; A3*h end proc

(1)

Q1 := diff(F(eta), eta, eta, eta)+.5*H(eta)*((diff(F(eta), eta))^2+F(eta)*(diff(F(eta), eta, eta)))/G(eta)^2+2*(diff(G(eta), eta))*(diff(F(eta), eta, eta))/G(eta)-(diff(H(eta), eta))*(diff(F(eta), eta, eta))/H(eta); 1; Q2 := diff(G(eta), eta, eta)+H(eta)*((diff(F(eta), eta))*G(eta)+.5*F(eta)*(diff(eta, eta)))/G(eta)^2+2*(diff(G(eta), eta))^2/G(eta)-((diff(H(eta), eta))*(diff(H(eta), eta)))/H(eta)+(diff(F(eta), eta, eta))^2-(H(eta)/G(eta))^2; 1; Q3 := diff(H(eta), eta, eta)+(.5*1.3)*H(eta)*(5*(diff(F(eta), eta))*H(eta)+F(eta)*(diff(H(eta), eta)))/G(eta)^2+2*(diff(G(eta), eta))*(diff(H(eta), eta))/G(eta)-(diff(H(eta), eta))^2/H(eta)+(1.3*1.44)*H(eta)*(diff(F(eta), eta, eta))/G(eta)-(1.3*1.92)*(H(eta)/G(eta))^3

-2*A3+.65*A3*(5*A1*(-2*eta^2+eta+1)*A3*(-eta^2+1)-2*A1*(-(2/3)*eta^3+(1/2)*eta^2+eta)*A3*eta)/((-eta^2+1)*A2^2)+4*A3*eta^2/(-eta^2+1)+1.872*A3*A1*(-4*eta+1)/A2-2.496*A3^3/A2^3

(2)

Eq1 := int(Q1*f, eta = 0 .. 1);

-0.2600000000e-1*A3*(24.*A1*A2^2-65.*A1*A2*A3+64.*A3^2)/A2^3

(3)

sol := solve({Eq1 = 0, Eq2 = 0, Eq3 = 0}, {A1, A2, A3}); J := min(select(`>`, sol, 0))

Error, invalid input: `>` expects 2 arguments, but received 1

 

A11 := evalf(simplify(sol[1, 1])); A22 := evalf(simplify(sol[1, 2])); A33 := evalf(simplify(sol[1, 3]))

Error, invalid subscript selector

 

``


 

Download ReducedCantiler.mw

 

How to find the integral of (x+y)/(x+y+z) over the part of the unit ball  centered at the origin which lies in the positive octant { x>=0 , y>=0, z>=0 } ? Numeric calculations suggest Pi/9.

How to find the double integral of sin(x^2)*cos(y^2) over the disk of radius R which is centered at the origin? 

Here is my try 

restart; evalf(VectorCalculus:-int(sin(x^2)*cos(y^2), [x, y] = Circle(`<,>`(0, 0), 1), inert), 15);
                       0.722091449378409
identify(%);
                       0.722091449378409

 

Hello people in mapleprimes,

I want to modify the following expression a to b. For it, I used c code.
But, as seen in d, just using c returns a. Originally, I intended phi being between 0 and 1, so
that 1-phi >0. Isn't any better way to modify a to b? and it is better to have maple sort expression remembering 1-phi to be positive, and to have maple give any priority to it in modifying a to b.

a:=(-2*theta+1)*phi/(phi-1)+theta;
b:=(2*theta-1)*phi/(-phi+1)+theta;
c:=simplify(-1*numer(op(1,aa)))/simplify((-1)*denom(op(1,aa)))+theta;
d:=simplify(c)
(-2*theta+1)*phi/(phi-1) + theta

I will be very glad if you borrow me your idea.

Thank you in advance.

Hi All,

To get EOM of a specific type of vibrating shell structure, I constructed the Hamilton equation of the system with messy and complex integrand (so many derivatives and variables are included). Now I have to calculate the variation of the Hamilton equation (del(int((T+U-V),t=t0..t1)=0) to get the EOMs and BCs of the system.

Is there anybody to know how to take a variation of such bulky integral in Maple?

Is there anyone to have experience in finding the EOM of shell structure using variational calculus in Maple?

Any comment will be useful.

Regards,

 

For some strange reason my matrices are shown in a very bizarre format (see imgur link). All I want is to simply see the 6x6 matrix with the numerical values - any help would be greatly appreciated. http://imgur.com/a/qaQmb

How to put number 1,2,3...12  around a circle, like a clock face using Maple code ?

Right now I can only do this manually by plotting a circle, then save the file and use Windos paint to

mannual put in the numbers

Let us consider 

with(Statistics);
U := RandomVariable(DiscreteUniform(-10, 10)):
V := RandomVariable(DiscreteUniform(-10, 10)):
Probability(U^2-V^2 <= 1/9, numeric);
  0.

, whereas a positive number greater than 1/21 is expected. 

 

I have not been using Maple for the last couple of years because I have Maple installed on my big and heavy laptop. Such laptop forces me to sit at a desk when I work which I dont like. This has resulted in that I have only (for the most part) used my ipad for the last couple of years hence dropping my laptop and Maple. The problem with such approach has been that is more difficult to write scientific papers because you cant run a full version of word on the ipad hence you cant use a citation manager. There does not exist an app on the ipad that offer one. This leads to a lot of extra work!

I have therefore now bought a Microsoft Surface Pro 4 which is a tablet (like the ipad) but its a not a "toy" tablet like the ipad or android tablets because it has Windows 10 Pro (64-bit) which means that you can run full versions of pc software like word 2016 and Maple. The version of Surface Pro 4 I bought has an Intel Core m3-6Y30 processor with 2 cores, 4GB RAM, 128 GB of storage, no fan and a processor base frequency of 0.9 GHz. The m3-6Y30 processor has something called dynamic frequency scaling which means (as I understand it) that the processor can go up to 2.2 GHz when it is under heavy load. Since Maple recommend that you have a system with a cpu with at least 1 GHz to run Maple 2016 I am not anticipating any problems here.

http://www.maplesoft.com/products/system_requirements.aspx

However, I am not completly sure if maplesoft means that the base frequency needs to be 1 Ghz or not? Since my tablet has not been delivered yet I cant comment on how Maple runs on it but as I said I am not anticipating any problems.

Now to my question: When will Maple 2017 come out? Since it is already 2017 I dont want to buy an old version of Maple now that I want to update my Maple software.

I can't find in the Linear Algebra package the command for creating the duplication matrix.

The duplication matrix is the n * n rows and n * ( n + 1 ) / 2 columns that transforms vech(A)
to vec(A) where A is a symmetric n by n matrix (Magnus & Neudecker).

If it not yet in Maple itself, has someone already made a proc for it?

Harry Garst

 

I want  to find the parametrs that verify the below inequalities at the same time:

beta*(v+u)>phi*(alpha+lamda+u+delta)

phi>beta

beta> (alpha+lamda+u+delta)

phi* (alpha+lamda+u+delta)>beta*(u+v+alpha),

question : find the values of u,v,delta,lamda,phi,beta,alpha such as all these parameters are in the intervall ]0,1]

can anyone help me ?

 

hi

why this equation does not any answer?

thanks

s-s.mw
 

restart

``

eq:={-J*g[1]*(diff(w(x), x, x, x, x, x, x))+J*c[1]*(diff(w(x), x, x, x, x))+A*g[113113]*(diff(w(x), x, x, x, x))+(beta[11]*A*0)*`&Delta;T`*(diff(w(x), x, x))+2*b*f[1133]*(Pi/L)^2*(d[33]*lambda[3]*`&Delta;T`*L/mu[33]-2*f[1133]*a*Pi/L-P[3]*`&Delta;T`*L)*sin(Pi*x/L)*sinh(h*Pi/(2*L))/(2*cosh(h*Pi/(2*L))*(-a33+d[33]^2/mu[33])) = 0, w(0) = 0, w(L) = 0, (D(w))(0) = 0, (D(w))(L) = 0, ((D@@3)(w))(0) = 0, ((D@@3)(w))(L) = 0}

{-J*g[1]*(diff(diff(diff(diff(diff(diff(w(x), x), x), x), x), x), x))+J*c[1]*(diff(diff(diff(diff(w(x), x), x), x), x))+A*g[113113]*(diff(diff(diff(diff(w(x), x), x), x), x))+b*f[1133]*Pi^2*(d[33]*lambda[3]*`&Delta;T`*L/mu[33]-2*f[1133]*a*Pi/L-P[3]*`&Delta;T`*L)*sin(Pi*x/L)*sinh((1/2)*h*Pi/L)/(L^2*cosh((1/2)*h*Pi/L)*(-a33+d[33]^2/mu[33])) = 0, w(0) = 0, w(L) = 0, (D(w))(0) = 0, (D(w))(L) = 0, ((D@@3)(w))(0) = 0, ((D@@3)(w))(L) = 0}

(1)

dsolve(eq, w(x))

``


 

Download s-s.mw

 

Graphical Programming with MapleSim in Vector Mechanics to Structures 2D

At the present time before constructing or starting to develop a mechanical structures project it is necessary to model it using graphic programming; In this opportunity and used MapleSim as a computational tool belonging to the company Maplesoft. The modern approach to modeling and simulation makes the fabrication of complex designs easy to solve. We will cover some examples taken from the engineering being implemented in Maplesim with insertion of physical objects; To be seen in real time through video output; Then integrates with Maple to analyze the equations and data through the static and dynamic behavior of the fabricated. Solved methods of physical block components include functionality for many domains: rotational and translational mechanics, multi-body dynamics, logic, and structural blocks; With techniques like: Drag-and-Drop Physical Modeling Environment and Create Custom Components Directly From Their Equations, thus the systems that would take hours or days to build from equations; In principle they can be created in a fraction of time using MapleSim, so it can incorporate significantly more complex graphical algorithms. In MapleSim, I use the revolutionary multibody technology that perfectly combines advanced multi-domain modeling tools to provide all the functionality you need in one environment.

FAST_UNT_2017.pdf

Lenin Araujo Castillo

Ambassador Maple - Perú

 

 

Let:

f:=x->1/sqrt(2*Pi)*exp(-x^2/2);

I.e. f is a standard Gaussian PDF.

Then (in Maple 2016.1):

Int(convert(f(x)*f(y)*x*x*abs(x+y),piecewise,x),x=-infinity..infinity,y=-infinity..infinity):
evalf(%);

Returns:

1.692568751

However (again in Maple 2016.1):

int(convert(f(x)*f(y)*x*x*abs(x+y),piecewise,x),x=-infinity..infinity,y=-infinity..infinity):
evalf(%);

Returns:

-0.5641895835

This is clearly incorrect, as the integral of a positive function must be positive.

This also seems to be a problem in which ever version of Maple is used behind the scenes on this forum.

int(convert(1/sqrt(2*Pi)*exp(-x^2/2)*1/sqrt(2*Pi)*exp(-y^2/2)*x*x*abs(x+y),piecewise,x),x=-infinity..infinity,y=-infinity..infinity)

gives:

int(convert(1/sqrt(2*Pi)*exp(-x^2/2)*1/sqrt(2*Pi)*exp(-y^2/2)*x*x*abs(x+y),piecewise,x),x=-infinity..infinity,y=-infinity..infinity)

It is suggested  

hypergeom([1/3, 2/3], [3/2], (27/4)*z^2*(1-z)) = 1/z

if z > 1. Here is my try to prove that with Maple:


 

a := `assuming`([convert(hypergeom([1/3, 2/3], [3/2], (27/4)*z^2*(1-z)), elementary)], [z > 1])

-(1/((1/2)*(27*z^3-27*z^2+4)^(1/2)+(3/2)*z*(3*z-3)^(1/2))^(1/3)-1/((1/2)*(27*z^3-27*z^2+4)^(1/2)-(3/2)*z*(3*z-3)^(1/2))^(1/3))/(z*(3*z-3)^(1/2))

(1)

b := `assuming`([simplify(a, symbolic)], [z >= 1])

2*(-(12*(3*z+1)^(1/2)*z-12*z*(3*z-3)^(1/2)-8*(3*z+1)^(1/2))^(1/3)+(12*(3*z+1)^(1/2)*z+12*z*(3*z-3)^(1/2)-8*(3*z+1)^(1/2))^(1/3))/((3*z-3)^(1/2)*(12*(3*z+1)^(1/2)*z+12*z*(3*z-3)^(1/2)-8*(3*z+1)^(1/2))^(1/3)*(12*(3*z+1)^(1/2)*z-12*z*(3*z-3)^(1/2)-8*(3*z+1)^(1/2))^(1/3)*z)

(2)

plot(1/b, z = 1 .. 10)

 

simplify(diff(1/b, z), symbolic)

-48*(((3*z-2)*(3*z+1)^(1/2)+z*(3*z-3)^(1/2))*((12*z-8)*(3*z+1)^(1/2)-12*z*(3*z-3)^(1/2))^(1/3)+((12*z-8)*(3*z+1)^(1/2)+12*z*(3*z-3)^(1/2))^(1/3)*((-3*z+2)*(3*z+1)^(1/2)+z*(3*z-3)^(1/2)))/((3*z+1)^(1/2)*(3*z-3)^(1/2)*((12*z-8)*(3*z+1)^(1/2)+12*z*(3*z-3)^(1/2))^(2/3)*((12*z-8)*(3*z+1)^(1/2)-12*z*(3*z-3)^(1/2))^(2/3)*(((12*z-8)*(3*z+1)^(1/2)-12*z*(3*z-3)^(1/2))^(1/3)-((12*z-8)*(3*z+1)^(1/2)+12*z*(3*z-3)^(1/2))^(1/3))^2)

(3)

``


 

Download simplification.mw

5 6 7 8 9 10 11 Last Page 7 of 30