Maple 2016 Questions and Posts

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

The angle unit of the SI system is rad and its dimension is 1. In engineering angles are often scaled in degrees. Maple supports the unit "deg" or "arcdeg". However when calculating with these units it does not convert deg to rad automatically and replace rad by 1. Here is an example:

The result of evalf(B) should be 15.708 (= 5*pi), i.e. a value without dimension, neither arcdeg nor rad. How do I get the correct result?

Dear all!.       

I have an expression with ramdom variables. Can i use NLP to optimize it?. In this case, does Maple take the histogram of theses ramdom variables?. Thank!!!

What is the difference between the Expression Palette and the Layout Palette? First I thought that the Layout Palette is supporting the entering of names with e.g. sub- and superscripts, but I learned (the hard way) that the templates of the Layout Palette are "mathematically active".

What is the meaning of the colours green and violett in both palettes?

Can somebody of Maple users execute the following command

restart; pdsolve({diff(w(x, y, z), x)+diff(w(x, y, z), y, y)+2*(diff(v(x, y, z), x)-(diff(u(x, y, z), y))-2*w(x, y, z)) = diff(w(x, y, z), z, z), 3*(diff(u(x, y, z), x, x))+2*(diff(u(x, y, z), y, y))+2*(diff(v(x, y, z), x, y))+2*(diff(w(x, y, z), y)) = diff(u(x, y, z), z, z), 3*(diff(v(x, y, z), y, y))+2*(diff(v(x, y, z), x, x))+2*(diff(u(x, y, z), x, y))-2*(diff(w(x, y, z), x)) = diff(v(x, y, z), z, z)}, {u(x, y, z), v(x, y, z), w(x, y, z)})

restart; pdsolve({diff(w(x, y, z), x)+diff(w(x, y, z), y, y)+2*(diff(v(x, y, z), x)-(diff(u(x, y, z), y))-2*w(x, y, z)) = diff(w(x, y, z), z, z), 3*(diff(u(x, y, z), x, x))+2*(diff(u(x, y, z), y, y))+2*(diff(v(x, y, z), x, y))+2*(diff(w(x, y, z), y)) = diff(u(x, y, z), z, z), 3*(diff(v(x, y, z), y, y))+2*(diff(v(x, y, z), x, x))+2*(diff(u(x, y, z), x, y))-2*(diff(w(x, y, z), x)) = diff(v(x, y, z), z, z)}, {u(x, y, z), v(x, y, z), w(x, y, z)})

Error, (in simplify/normal) Maple was unable to allocate enough memory to complete this computation.  Please see ?alloc

 

 

 

in Maple 2016.1.1 on a powerful comp and report the obtained result as an answer to the question?
 That would be very kind of her/him. Thanks in advance. 

Download pdsolve.mw

I am learning how to use Maple with boundary value ODE. Given this ODE

y''''(x)+ lam* y(x) =0

with some B.C., say  y(0)=0,y'(0)=0,y''(L)=0,y'''(L)=0, where L is length.

I can't figure the correct syntax to use. It seems Maple do not like the syntax I am using, but it works on a second order ODE?

Here is my attempt:

restart;
assume(lam>0); assume(L>0);
bc:=y(0)=0,D[1](y)(0)=0,D[2](y)(L)=0,D[3](y)(L)=0;
dsolve({diff(y(x),x$4)+lam*y(x)=0,bc},y(x));

Error is 

Error, (in evalapply) too few variables for the derivative with respect to the 2nd variable
Error, (in dsolve) found the following equations not depending on the unknowns of the input system: {bc}

But on a simpler second order ODE, the syntax works

restart;
assume(lam>0); assume(L>0);
bc:=y(0)=0,D[1](y)(0)=0;
dsolve([diff(y(x),x$2)+lam*y(x)=0,bc],y(x));

No error. 

Is the syntax I am using in first example wrong? what would be the correct syntax? I googled for long time, and can't find one example that shows how to use BVP with higher order ODE. I am Maple newbie.

 

Hello people in mapleprimes,

I put a function h into my original package named z5.

And, I used it as

> h(x+y);

then, the output is

x*`z5:-h`(x)/(x+y)+y*`z5:-h`(y)/(x+y)

I want to have

x*h(x)/(x+y)+y*h(y)/(x+y),

not z5:-h for h.

How should I do for "h" to be kept from changing  to "z5:-h"?

Best wishes.

taro

 

part other than z5 needs designating a directory for library.**********

>new_lib_dir3 :=Please put an appropriate directory

>z5:=module() option package;
>export numer_expand,h;
>numer_expand := a->expand(numer(a))/denom(a); h:= F-> expand
>(evalindets(D(F)/F, specfunc(D), d-> op(d)*'h'(op(d)))); end module;

>LibraryTools:-Save('z5',new_lib_dir3);
>restart;

>libname:=new_lib_dir3,libname;

>with(z5);

>h(x+y);

               x*`z5:-h`(x)/(x+y)+y*`z5:-h`(y)/(x+y)

 

I have a problem with the right-click command „Unit Formatting“. Here is an example: I am trying to get in result in kW instead of W. This works as long as the value is assigned to the variable using := . If the value is calculated I cannot change the unit.

Right-clicking P and P2 and selecting „Unit Formatting“ I get the following selection menu

Right-clicking P3 offers no units to be selected. Entering a Custom Unit returns an error message.

Is that behaviour normal for Maple or am I making a mistake?

 

Hi,

I have reinstalled Maple 2016 on a new computer, the installation is fine. However, with time, on the old computer, I made several adjustments in the Tools->Options menu, I have added tasks, selected some palettes, removed some others, etc... Is there a way to transport those settings all at once, without performing all that again?

Thanks.

Hey, I have a problem loading one of my maple files as the headline states.

I have seen a sugestion on how to fix it:

http://www.mapleprimes.com/questions/129377-There-Were-Problems-During-The-Loading

but I do not know how to use it proberly, and can therefore not fix my problem.

What do I do?

Elektriske_felter.mw


 

Dear Community,

Would someone have a good and easy to understand/implement description of the Den Iseger algorithm for the numerical inversion of Laplace transform? Even better if someone would have a Maple script to do it, that would be superb.

Tx in advance,

best regards

Andras

hi...

how i can save results in text file?

for example in attached code maple i want save results for different ''x'' and fy(x), wherer x is between 0 to 1 i.e. (0,fy(0)) , (0.1,fy(0.1)) ,(0.2,fy(0.2)),....

in text file two column should be create.one of them for ''x'' and the onther for ,fy(x)

thanks...

Numerical_.mw

restart; E := 0.169e12; mu := 0.658e11; hl2 := 4; D1 := 1; n := 3; `αn` := 0.; beta := 16.474184; xi := 1.5; lambda := .1; dsys5 := {D1*(diff(y(x), x, x, x, x)) = `αn`/(1-y(x))^n-beta*(lambda*cosh(xi*(1-y(x)))-1/2*(lambda^2+1))/sinh(xi*(1-y(x)))^2, y(0) = 0, (D(y))(0) = 0, ((D@@2)(y))(1) = 0, ((D@@3)(y))(1) = 0}; dsol5 := dsolve(dsys5, 'maxmesh' = 900, numeric, output = listprocedure); fy := eval(y(x), dsol5); fy(1)

HFloat(0.3189919226404612)

(1)

NULL

NULL

NULL

NULL



Download Numerical_.mw

 

Hi

may every one help to me for dsolve this differentia1l equation?

error:

Error, (in dsolve/numeric/bvp) singularity encountered

Turbulent2-kw.mw

dsol1 := dsolve({diff(theta(eta), eta, eta)-3*Omega(eta)*(F(eta)*(diff(theta(eta), eta))-theta(eta)*(diff(F(eta), eta)))/(2*K(eta))+((diff(K(eta), eta))/K(eta)-(diff(Omega(eta), eta))/Omega(eta))*(diff(theta(eta), eta)) = 0, diff(F(eta), eta, eta, eta)+Omega(eta)*(3*F(eta)*(diff(F(eta), eta, eta))-(diff(F(eta), eta))^2)/(2*K(eta))+((diff(K(eta), eta))/K(eta)-(diff(Omega(eta), eta))/Omega(eta))*(diff(F(eta), eta, eta))+Omega(eta)/K(eta) = 0, diff(K(eta), eta, eta)+Omega(eta)*(1.5*F(eta)*(diff(K(eta), eta))-K(eta)*(diff(F(eta), eta)))/K(eta)+((diff(K(eta), eta))/K(eta)-(diff(Omega(eta), eta))/Omega(eta))*(diff(K(eta), eta))+(diff(F(eta), eta, eta))^2-Omega(eta)^2 = 0, diff(Omega(eta), eta, eta)+Omega(eta)*(3*F(eta)*(diff(Omega(eta), eta))+Omega(eta)*(diff(F(eta), eta)))/(2*K(eta))+((diff(K(eta), eta))/K(eta)-(diff(Omega(eta), eta))/Omega(eta))*(diff(Omega(eta), eta))+Omega(eta)*(diff(F(eta), eta, eta))^2/K(eta)-Omega(eta)^3/K(eta) = 0, F(0) = 0, K(0) = 0, Omega(0) = 0., theta(0) = 1, theta(1) = 0, (D(F))(0) = 0, (D(K))(1) = 0, (D(Omega))(1) = 0, ((D@@2)(F))(1) = 0}, numeric, method = bvp[middefer], output = listprocedure, initmesh = 512)

Error, (in dsolve/numeric/bvp) singularity encountered

 

NULL

plots[odeplot](dsol1, [(D(F))(eta), eta])

Error, (in plots/odeplot) input is not a valid dsolve/numeric solution

 

``



Download Turbulent2-kw.mw

 

 

 

Hello:

I want to make a package or module. Not sure. I cant get my head around what to do. I have loads of formulas I would like to be able to use without going back through work sheets. At present just want to keep it simple.

Have loaded a sample work sheet to use

Intend to place package in Maple2016 library folder.

 

restart

Library Folder

C:\Program Files\Maple 2016\lib

``

``

with(LinearAlgebra):

``

``

Points

Point2 := proc (x1, y1) options operator, arrow; [x1, y1] end proc

proc (x1, y1) options operator, arrow; [x1, y1] end proc

(1)

Point3 := proc (x1, y1, z1) options operator, arrow; [x1, y1, z1] end proc

proc (x1, y1, z1) options operator, arrow; [x1, y1, z1] end proc

(2)

a1 := Point3(x1, y1, z1)

[x1, y1, z1]

(3)

a2 := Point3(x2, y2, z2)

[x2, y2, z2]

(4)

a3 := Point3(x3, y3, z3)

[x3, y3, z3]

(5)

Line through 2 points

LinePts := proc (p1, p2) options operator, arrow; (p2[2]-p1[2])*x+(p1[1]-p2[1])*y-p2[2]*p1[1]+p1[2]*p2[1] end proc

proc (p1, p2) options operator, arrow; (p2[2]-p1[2])*x+(p1[1]-p2[1])*y-p2[2]*p1[1]+p1[2]*p2[1] end proc

(6)

``

Line

Line := proc (a, b, c) options operator, arrow; a*x+b*y+c end proc

proc (a, b, c) options operator, arrow; a*x+b*y+c end proc

(7)

Lines s

LnPrll := proc (l1, l2) options operator, arrow; is(coeff(l1, x)*coeff(l2, y)-coeff(l2, x)*coeff(l1, y) = 0) end proc

proc (l1, l2) options operator, arrow; is(coeff(l1, x)*coeff(l2, y)-coeff(l2, x)*coeff(l1, y) = 0) end proc

(8)

Lines t

LnPrpnd := proc (l1, l2) options operator, arrow; is(coeff(l1, x)*coeff(l2, x)+coeff(l1, y)*coeff(l2, y) = 0) end proc

proc (l1, l2) options operator, arrow; is(coeff(l1, x)*coeff(l2, x)+coeff(l1, y)*coeff(l2, y) = 0) end proc

(9)

``

``

Quadrance between 2 points (2D)

Quadrance := proc (a1, a2) options operator, arrow; (a1[1]-a2[1])^2+(a1[2]-a2[2])^2 end proc

proc (a1, a2) options operator, arrow; (a1[1]-a2[1])^2+(a1[2]-a2[2])^2 end proc

(10)

``

Quadrance between 2 points (3D)

Qd3 := proc (a1, a2) options operator, arrow; (a1[1]-a2[1])^2+(a1[2]-a2[2])^2+(a1[3]-a2[3])^2 end proc

proc (a1, a2) options operator, arrow; (a1[1]-a2[1])^2+(a1[2]-a2[2])^2+(a1[3]-a2[3])^2 end proc

(11)

Triple Quad Formula (3 colinear points)

TQF := proc (Q1, Q2, Q3) options operator, arrow; (Q1+Q2+Q3)^2 = 2*Q1^2+2*Q2^2+2*Q3^2 end proc

proc (Q1, Q2, Q3) options operator, arrow; (Q1+Q2+Q3)^2 = 2*Q1^2+2*Q2^2+2*Q3^2 end proc

(12)

Spread*PolynomialsNULL

SpreadPoly := proc (n, s, rn) with(orthopoly); s-factor(1/2-(1/2)*T(n, 1-2*rn)) end proc

proc (n, s, rn) with(orthopoly); s-factor(1/2-(1/2)*T(n, 1-2*rn)) end proc

(13)

SpreadPoly(3, s, r)

s-r*(4*r-3)^2

(14)

``

 

Download Testpackage.mw

I updated the OSX from El Capitan 10.11 to Sierra 10.12.

After then, whenever I try to input Japanese characters, Internal Error notification appears, and it crashes.
As it is not realistic to go back to the old OS, as the time when I did back up was a little before, I wouldn't do so.
Then, I might end up torelating not using Japanese until the next version of maple appears next year.
Are the situations like this ?

I know there might not be any response as this is about inputting Japanese characters.

Best wishes.

taro

A student of mine has a problem, when trying to open a *.mw file directly from Finder, by double clicking og right-click and choose Open or Open with.  

Maple will prompt - the file does not exist.

 

If she uses Maple and and opening the same file thru file -> open etc. There is no problem. 

 

Any suggestions?

Kind regards 

Per Kirkegaard

 

 

First 44 45 46 47 48 49 50 Last Page 46 of 60