Maple 2015 Questions and Posts

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

Hello. I want to replicate the results of this table in Maple.

My "attempt" uses the inbuilt dictionary, but it will only find the words of length n, not the frequency of letters in each ....

word_length.mw

Let say I have an array

A := Array([[1], [2], [3]])
Normally if I would like to extend it I will write, eg:

Extend(A, [[4], [5]])

But because this is not 1D array Maple will rise an error. Is there a way to go over this limitation?

hi

in attached file below not answer found for dsolve?

please help me

thanks..

dsove.mw

restart

J := 1:

PDE := diff(T(z, t), z, z)-.2*(diff(T(z, t), t, t)) = int(.7703831837*(diff(T(z, tau), tau, tau))/(t-tau)^.3, tau = 0 .. t):

with(inttrans):

sol := laplace(PDE, t, s):

sol2 := subs([laplace(T(z, t), t, s) = U(z, t), T(z, 0) = sin(J*Pi*z), (D[2](T))(z, 0) = 0], sol)

diff(diff(U(z, t), z), z)-.2000000000*s^2*U(z, t)+.2000000000*s*sin(Pi*z) = -1.000000000*s^.3000000000*sin(Pi*z)+1.000000000*s^1.300000000*U(z, t)

(1)

sol3 := dsolve([sol2, (D[1](U))(0, t) = 0, (D[1](U))(1, t) = 0], U(z, t))

"sol3 := "

(2)

U(z, t) = invlaplace(rhs(sol3), s, t)

Error, invalid input: rhs expects 1 argument, but received 0

 

sol4 := simplify(subs(z = 0, rhs(sol3)))

Error, invalid input: rhs expects 1 argument, but received 0

 

``

``



Download dsove.mw

 

 

 

so I'm trying this:

restart;

sigma := 0.143e-18;

l_0 := 1.87;

l0 := 1.87;

roll := rand(0 .. 25.0);

f_gauss := proc (x) options operator, arrow; exp(-(1/2)*x^2/`σ_x`^2)/sqrt(2*Pi*`σ_x`^2) end proc;

f_norm := proc (dx) options operator, arrow; int(f_gauss(x), x = -(1/2)*dx .. (1/2)*dx) end proc;

sol_gauss := proc (mix) options operator, arrow; evalf(eval(-ln((int(f_gauss(x)*exp(-2*sigma*N2O*sqrt((1/4)*l_0^2-x^2)), x = -(1/2)*dx .. (1/2)*dx))/f_norm(dx))/(sigma*N2O), [N2O = 0.25e20*mix/100])) end proc;

for ii to 10 do

a := roll();

eval(sol_gauss(a), [dx = l_0, `σ_x` = l0])

end do

Already searched and browsed multiple different threads and still cannot find a solution.

Apologizing the noob nature of this question.

 

 

In this code below, Why is the factor command not working?

f := a^2+x^2-2*ax;

a^2+x^2-2*ax

(1)

factor(f);

a^2+x^2-2*ax

(2)

expand((x-a)*(x-a));

a^2-2*a*x+x^2

(3)

``

NULL

NULL


Download factor_polynomial_2_multivariable.mwfactor_polynomial_2_multivariable.mw

 

hi...how i can pdsolve this equation numerically or analyticlly?

this equation is time-fractional  equation with generalized Cattaneo model

where

 

 is the fractional derivative operator considered in the
Caputo sense.

 

FRACTION.mw

restart

k := 1; -1; rho := 1; -1; h := 1; -1; alpha := 2-Upsilon; -1; 0 < Upsilon and Upsilon <= 1

0 < Upsilon and Upsilon <= 1

(1)

k*(diff(T(z, t), z, z)) = rho*(diff(T(z, t), [`$`(t, alpha)]))

diff(diff(T(z, t), z), z) = diff(T(z, t), [`$`(t, 2-Upsilon)])

(2)

k*(diff(T((1/2)*h, t), z)) = 1:

k*(diff(T((-h)*(1/2), t), z)) = 0:

T(z, 0) = 0

T(z, 0) = 0

(3)

NULL



Download FRACTION.mw

 

Dear Pros, I'm a biginer so I have a question about my program.

I have a lot of arrays which are result from 2 while loop. Now, I want creat a matrix from them but i can't. So, could you help me to do it.

For detail: 

V[1]:=[ 1 2 3]

V[2]:=[2 3 4]

V[3]:=[3 4 5]

V[4]:=[2 6 7]

V[5]:=[7 8 9]

...

V[n]

with type of V[i] is a array.

I searched and found a solution by manual to create a matrix as follow:

V_matrix:=<V[1],V[2],V[3]>

but in this case i can't but manual with n=100

please help me to have a Matrix.

Thank a lots.

Hi everybody, 

 

Could somebody just explain me why commands L[1..5] and L[3..6] below do not seem to return "the same type of thing" ?

 

L := Array(1..10, [$(1..10)]);

L := Array(1..10, {(1) = 1, (2) = 2, (3) = 3, (4) = 4, (5) = 5, (6) = 6, (7) = 7, (8) = 8, (9) = 9, (10) = 10})

(1)

L[1..5];
whattype(%);

Array([1, 2, 3, 4, 5])

 

Array

(2)

L[3..6];
whattype(%);

Array(3 .. 6, {3 = 3, 4 = 4, 5 = 5, 6 = 6}, datatype = anything, storage = rectangular, order = Fortran_order)

 

Array

(3)

 


Thank you for your attention


Download array.mw

how to specify color at the turning point of mountain shape in dualaxisplot

woud like to specify a point with another color at the top of peak shape line diagram

After using the Groebner and PolynomialIdeals packages, Maple goes into a long calculation when I make an entry of the form

name:=polynomial expression. This can take 10's of minutes for an expression of two lines. The only solution I have found is to save the sheet and restart it and enter the line name:= etc. before loading Groebner and PolynomialIdeals. This is most inconvenient. Is there a better workaround?

how to change rotation angle into slope of line?

with(MTM):
evalf(solve(atan(-1) - 45*Pi/180 = atan(x), x))*180/Pi;

above can not be solved

 

if angle is 30 degree

is following correct?

with(MTM):
angle := 30:
evalf((atan(-1) - 30*Pi/180)/Pi)*(-1.0);

I would like compute the vertices of the convex hull of finitely many points in some \mathbb{R}^{n}, where the coordinates of these points could be rational numbers or irrational numbers like \sqrt{2}.

I know that PolyhedralSets[ConvexHull] is a direct command, however it does not support irrational numbers. But if I transform all coordinates into floating point numbers, for some rational numbers the rounding error would be large enough to affect the result. For example, after the transform some collinear points are no longer collinear.

Thanks!

 

 

 

how to use Riemann matrix to output Riemann surface?

and plot this surface?

 

with(algcurves):
f:=y^3+2*x^7-x^3*y;
pm:=periodmatrix(f,x,y);
evalf(pm, 5);
rm:=evalf(periodmatrix(f,x,y,Riemann),10);
M := rm;
A := proc (x, y) options operator, arrow; RiemannTheta([x, y], M, [], 0.1e-1, output = list)[2] end proc;
plot3d(Re(A(x+I*y, 0)), x = 0 .. 1, y = 0 .. 4, grid = [40, 40]);

is this graph Riemann Surface?

if so, how to convert A into polynomials?

Hello fellow maple struglers!

This is my first question, so hopefully I did it in the correct manner. My problem (I think) is of a simple nature, but I can't seem to get it to work. The equation I have is the following:

All is nice and well and as you all can see it depends on two variables. Now what I want to do is to express the function in a shape where it depends on the ratio t/d only. Simple calculus shows that for this equation it is rather easy, where the result is:

(62,8*t/d+0,0945)/pi

Now I would like to know whether it is possible to automate this in maple, such that I can plot Wr_over_p as a function of t/d. I know I can simply use the form above and make the plot, but further down in my maple file I have an equation which is much more complex, and I want to plot thatone as a function of t/d as well.

I hope I have made my problem clear to some of you guys, but if you require some further explaination than please let me know. Thanks in advance for your help, I hope there is someone who can help me with this!

First 37 38 39 40 41 42 43 Last Page 39 of 71