Maple 2016 Questions and Posts

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

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

 

 

How does one obtain all solutions from dsolve? I see an option called Allsolutions, but this seems to only apply to solve and other functions. It does not work with dsolve.

For example, maple gives one solution for the following first order non-linear ODE. But the ODE has another solution y(x)=0 as well. How does one tell Maple to return all solutions? I am interested in this when using the 'implicit' option mainly. Here is an example

restart;
num:=-(exp(x)*sin(y(x))-2*y(x)*sin(x)):
den:=(exp(x)*cos(y(x))+2*cos(x)):
eq:=diff(y(x),x)=num/den;

r0:=dsolve(eq,y(x));

But when I tried y(x)=0, it turned out to also be a solution

odetest(y(x)=0,eq);
                          0

But dsolve did not return this solution on its own along with the first one.  But on another example, Maple did well, and returned all solutions. Here is the other example

eq:=(2*x*y(x)^2+2*y(x))+(2*x^2*y(x)+2*x)*diff(y(x),x);
dsolve(eq=0,y(x),'implicit');

In the above, Maple returned the two solutions. 

Is there a correct way to tell Maple dsolve to return all solutions all the time? Why did it return both solution in the above example, but not in the first example?

I am maple newbie. Thank you.

hi

how i can draw this equation in maple in 3D?

(x-y)^2+(x-z)^2+(z-x)^2=3

thanks

PLOT3D.mw


 

Hello there. I have to solve a simple line integral (3*y-x)dL from A(2;1) to B(3;-1). A and B are points in the line described by function y=5-2x. I found a function LineInt. Using paramether Line I wrote LineInt(3*y-x, Line(<2, 1>, <3, -1>)), it doesn't work with error below:

Error, (in Line:-ModuleCopy) invalid input: invalid specification of line

It seems that the problem is with first argument, that is the equation under the integral sign. Or no... How to write it properly? Thanks.

I mean the root of the equation

GAMMA(n-1/n)*GAMMA(1/n)/(n*GAMMA(n)) = 1

belonging to RealRange(Open(1),4). It should be noticed there are solutions outside this interval. Here is my try.

 

``

solve({n > 1, GAMMA(n-1/n)*GAMMA(1/n)/(n*GAMMA(n)) = 1, n < 4}, [n])``

[]

(1)

`in`(which*is*wrong, view*of)

simplify(eval(GAMMA(n-1/n)*GAMMA(1/n)/(n*GAMMA(n)), n = (1/2)*sqrt(5)+1/2))

1

(2)

Also

Student[Calculus1]:-Roots(A = 1, n = 1 .. 4)

[1.618033989]

(3)

There is a substitute

fsolve(GAMMA(n-1/n)*GAMMA(1/n)/(n*GAMMA(n)) = 1, n = 1 .. 4)

1.618033989

(4)

NULL

identify(%)

(1/2)*5^(1/2)+1/2

(5)

``

There is a shade of hope that GAMMA(n-1/n)*GAMMA(1/n)/(n*GAMMA(n))  can be simplified.

Download solution.mw

 PS. An SCR was submitted by me.

The following product

 

(product(mu^x[i]/factorial(x[i]), i = 1 .. n))

 

does not simplify to the most obvious form whatever I try

 

mu^(sum(x[i], i = 1 .. n))/(product(factorial(x[i]), i = 1 .. n))

 

What can it be?

 

 

First 45 46 47 48 49 50 51 Last Page 47 of 61