Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Hello,

Is it possible to create animation of convolution of two functions?

For example f(t)=u(t)-u(t-2) and g(t) = tu(t)-(t-4)u(t-4), where u(t) is a step function.

I would like to generate animation for this convolution.

Any help would be appreciated.

 

Thanks.

 

     Maple is seriously used in my article Approximation of subharmonic functions in the half-plane by the logarithm of the modulus of an analytic function. Math. Notes 78, No 4, 447-455 in two places. The purpose of this post is to present these applications.                                                                                                 First, I needed to prove the elementary inequality (related to the properties of the minimal harmonic majorant of the function 1/Im z in a certain strip)                                                                                                    2R+sqrt(R)-R(R+sqrt(R))y - 1/y   1/4                                                                                                  for    y ≥ 1/(R+sqrt(R)) and  y ≤ 1/R, the parameter R is greater than or equal to 1.   The artless attemt                                                                          
restart; `assuming`([maximize(2*R+sqrt(R)-R*(R+sqrt(R))*y-1/y, y = 1/(R+sqrt(R)) .. 1/R)], [R >= 1])

maximize(2*R+R^(1/2)-R*(R+R^(1/2))*y-1/y, y = 1/(R+R^(1/2)) .. 1/R)

(1)

fails. The second (and successful) try consists in the use of optimizers:

F := proc (R) options operator, arrow; evalf(maximize(2*R+sqrt(R)-R*(R+sqrt(R))*y-1/y, y = 1/(R+sqrt(R)) .. 1/R)) end proc:

F(1)

.171572876

(2)

 

Optimization:-Minimize('F(R)', {R >= 1})

[.171572875253809986, [R = HFloat(1.0)]]

(3)

To be sure ,
DirectSearch:-Search(proc (R) options operator, arrow; F(R) end proc, {R >= 1})
;

[.171572875745665, Vector(1, {(1) = 1.0000000195752754}, datatype = float[8]), 11]

(4)

Because 0.17
"158 < 0.25, the inequality is  proved.   "
Now we establish this  by the use of the derivative. 

solve(diff(2*R+sqrt(R)-R*(R+sqrt(R))*y-1/y, y) = 0, y, explicit)

1/(R^(3/2)+R^2)^(1/2), -1/(R^(3/2)+R^2)^(1/2)

(5)

maximize(1/sqrt(R^(3/2)+R^2)-1/(R+sqrt(R)), R = 1 .. infinity, location)

(1/2)*2^(1/2)-1/2, {[{R = 1}, (1/2)*2^(1/2)-1/2]}

(6)

minimize(eval(2*R+sqrt(R)-R*(R+sqrt(R))*y-1/y, y = 1/sqrt(R^(3/2)+R^2)), R = 1 .. infinity, location)

3-2*2^(1/2), {[{R = 1}, 3-2*2^(1/2)]}

(7)

evalf(3-2*sqrt(2))

.171572876

(8)

The second use of Maple was the calculation of the asymptotics of the following integral (This is the double integral of the Laplacian of 1/Im z over the domain {z: |z-iR/2| < R/2} \ {z: |z| ≤ 1}.). That place is the key point of the proof. Its direct calculation in the polar coordinates fails.

`assuming`([(int(int(2/(r^2*sin(phi)^3), r = 1 .. R*sin(phi)), phi = arcsin(1/R) .. Pi-arcsin(1/R)))/(2*Pi)], [R >= 1])

(1/2)*(int(int(2/(r^2*sin(phi)^3), r = 1 .. R*sin(phi)), phi = arcsin(1/R) .. Pi-arcsin(1/R)))/Pi

(9)

In order to overcome the difficulty, we find the inner integral

`assuming`([(int(2/(r^2*sin(phi)^3), r = 1 .. R*sin(phi)))/(2*Pi)], [R*sin(phi) >= 1])

(R*sin(phi)-1)/(sin(phi)^4*R*Pi)

(10)

and then we find the outer integral. Because
`assuming`([int((R*sin(phi)-1)/(sin(phi)^4*R*Pi), phi = arcsin(1/R) .. Pi-arcsin(1/R))], [R >= 1])

int((R*sin(phi)-1)/(sin(phi)^4*R*Pi), phi = arcsin(1/R) .. Pi-arcsin(1/R))

(11)

is not successful, we find the indefinite integral  

J := int((R*sin(phi)-1)/(sin(phi)^4*R*Pi), phi)

-(1/2)*cos(phi)/(Pi*sin(phi)^2)+(1/2)*ln(csc(phi)-cot(phi))/Pi+(1/3)*cos(phi)/(R*Pi*sin(phi)^3)+(2/3)*cos(phi)/(R*Pi*sin(phi))

(12)

We verify that  the domain of the antiderivative includes the range of the integration.
plot(-cos(phi)/sin(phi)^2+ln(csc(phi)-cot(phi)), phi = 0 .. Pi)

 

plot((2/3)*cos(phi)/sin(phi)^3+(4/3)*cos(phi)/sin(phi), phi = 0 .. Pi)

 

    That's all right. By the Newton-Leibnitz formula,

``
eval(J, phi = Pi-arcsin(1/R))-(eval(J, phi = arcsin(1/R)));

(1/3)*(1-1/R^2)^(1/2)*R^2/Pi+(1/2)*ln((1-1/R^2)^(1/2)*R+R)/Pi-(4/3)*(1-1/R^2)^(1/2)/Pi-(1/2)*ln(R-(1-1/R^2)^(1/2)*R)/Pi

(13)

Finally, the*asymptotics*is found by

asympt(eval(J, phi = Pi-arcsin(1/R))-(eval(J, phi = arcsin(1/R))), R, 3)

(1/3)*R^2/Pi-(3/2)/Pi+(1/2)*(ln(2)+ln(R))/Pi-(1/2)*(-ln(2)-ln(R))/Pi+O(1/R^2)

(14)

      It should be noted that a somewhat different expression is written in the article. My inaccuracy, as far as I remember it, consisted in the integration over the whole disk {z: |z-iR/2| < R/2} instead of {z: |z-iR/2| < R/2} \ {z: |z| ≤ 1}. Because only the form of the asymptotics const*R^2 + remainder is used in the article, the exact value of this non-zero constant is of no importance.

       It would be nice if somebody else presents similar examples here or elsewhere.

 

Download Discovery_with_Maple.mw

Does anyone know how to using some softward convert ?  for example  I have maple code , but I want to using mathematica code .I need fast way.

I have 16 running on 32 bit xp computer, but I cannot install on 64 bit m/c with windows 8.1.

Have tried -r ,-f, -i options as in install.html, but install log says nothing installed.

How do I install ?

I cannot make Maple 13 read the maple.ini file no matter where it is located.

I am using Windows 7, 64 bit.

I thought that ...maple 13\bin.win would do it.

The file contains a few commands of mine. Three of these commands

"read" three procedure files A.m, B.m C.m. respectively, and the other two

are

with(LinearAlgebra):
with(DEtools):

Could somebody up here direct me properly about this?

I would appreciate it.

Thank you!

mapleatha

 

There is no menu item called "Startup Code"

in the "Edit" menu item of my Maple 13.

Please help!

 

Thanks!

I would like to create a variable name that consists of a greek letter and an english letter. For example, δx.   If I type

>delta = 5;

Maple displays the output as δ = 5

However, I would like to have variable names that are output as δx.

Is there a way to do this in worksheet mode? I tried cat(delta,x) but that just gave me deltax.

 

Thanks for your help,

Harish

hi friends

i have a problem in maple with an error

dsnumsort:=proc(numpr::list,coor::list)
local i,j,n;
global C1,C2,C3,V1,V2,V3;
n:=nops(coor):
print("Order of the variables:");
for i from 2 to 2*n +1 do
for j from 1 to n do
if[numpr[i]]=
select(has,numpr,diff(coor[j](t),t)) then
C[j]:=i-1;V[j]:=i;
print(coor[j],C[j]," ",diff(coor[j](t),t),V[j]);
fi;
od;
od;
end :

Us:=subs(G=1,m=1,L=1,U):
D2r:=[diff(z(t),t,t),diff(x(t),t,t),diff(y(t),t,t)]:
g:=subs(z=z(t),x=x(t),y=y(t),grad(Us,[z,x,y])):
IniC:=z(0)=0.75,D(z)(0)=0, x(0)=1,D(x)(0)=0,y(0)=0,D(y)(0)=1:

Ns:=dsolve({seq(D2r[i]=g[i],i=1..3),IniC},{z(t),x(t),y(t)},numeric);

dsnumsort(Ns(0),[z,x,y]):

for i from 0 to 1000 do ;
T:=i/25;
NsT:=Ns(T):

X[i]:=rhs(NsT[C1]); Vx[i]:=rhs(NsT[V1]);
Y[i]:=rhs(NsT[C2]); Vy[i]:=rhs(NsT[V2]);
Z[i]:=rhs(NsT[C3]); Vz[i]:=rhs(NsT[V3]);

KepVec[i]:=convert(crossprod([X[i],Y[i], Z[i]],[Vx[i],Vy[i], Vz[i]]),list);
KepAbs[i]:=norm(KepVec[i],2);
od:

but i see this error and I can't draw PLOT:

Error, invalid input: rhs received Ns(0)[C1], which is not valid for its 1st argument, expr

this cods is for draw plot:

spacecurve({[seq([X[i], Y[i], Z[i]], i = 0 .. 1000)], [[-1/2, 0, 0], [1/2, 0, 0]]}, labels = ['x', 'y', 'z']);

spacecurve([seq(KepVec[i], i = 0 .. 1000)], orientation = [0, 90], labels = ['x', 'y', 'z'])

plot([seq([(1/25)*i, KepAbs[i]], i = 0 .. 1000)], labels = ['t', 'MofI'])

can you helpe me?Thank you

I have a surfdata-plot which - by interpolation - goes through different circles with different inclinations with respect to each other. The surface is color-coded with zhue from violet to red. In order to see the valleys behind peaks I would like to have a filled 2D-plot of my zhue-surfdata colors with contours wherever the surface passes z-coordinate integers located somewhere below the surface (e.g. at z=-7). Unfortunately, I haven't found examples of contourplots for surfdata plots so far.

restart

with(plots);

[animate, animate3d, animatecurve, arrow, changecoords, complexplot, complexplot3d, conformal, conformal3d, contourplot, contourplot3d, coordplot, coordplot3d, densityplot, display, dualaxisplot, fieldplot, fieldplot3d, gradplot, gradplot3d, graphplot3d, implicitplot, implicitplot3d, inequal, interactive, interactiveparams, intersectplot, listcontplot, listcontplot3d, listdensityplot, listplot, listplot3d, loglogplot, logplot, matrixplot, multiple, odeplot, pareto, plotcompare, pointplot, pointplot3d, polarplot, polygonplot, polygonplot3d, polyhedra_supported, polyhedraplot, rootlocus, semilogplot, setcolors, setoptions, setoptions3d, spacecurve, sparsematrixplot, surfdata, textplot, textplot3d, tubeplot]

(1)

with(CurveFitting);

[ArrayInterpolation, BSpline, BSplineCurve, Interactive, LeastSquares, PolynomialInterpolation, RationalInterpolation, Spline, ThieleInterpolation]

(2)

R := [30, 32.5, 37.5, 42.5, 47.5, 50]:

incl := [0, .5, 1, -2, 5, 0]:

phases := [0, (1/12)*Pi, (1/8)*Pi, (1/4)*Pi, (1/2)*Pi, 0]:

colors := [grey, black, blue, red, green, grey]:

orbit := [R[j]*cos(t), R[j]*sin(t), incl[j]*cos(t-phases[j])]:

display(seq(spacecurve(orbit, t = 0 .. 2*Pi, color = colors[j], view = [-56 .. 56, -56 .. 56, -15 .. 15], labels = ["x [AU]", "y [AU]", "z [AU]"]), j = 1 .. 6));

 

NULL

R := [30, 32.5, 37.5, 42.5, 47.5, 50]:

incl := [0, .5, 1, -2, 2*evalf(Pi)-1.75, 0]:

pointplot(R, incl, color = [grey, red, blue, green, black, grey], labels = ["radius", "incl"]);

 

NewR := [seq(30+.4*i, i = 0 .. 50)]:

Newincl := ArrayInterpolation(R, incl, NewR, method = spline):

pointplot(NewR, Newincl, labels = ["radius", "incl"]);

 

phases := [0, (1/12)*evalf(Pi)+.3, 3*evalf(Pi)*(1/4), (1/4)*evalf(Pi)+.2, evalf(Pi)/(2.5)+.5, 0]:

pointplot(R, phases, labels = ["radius", "phase"]);

 

Newphases := ArrayInterpolation(R, phases, NewR, method = spline):

pointplot(NewR, Newphases, labels = ["radius", "phase"]);

 

t := [seq(0+i*(2*Pi*(1/50)), i = 0 .. 50)]:

f := proc (i, j) options operator, arrow; [NewR[i]*cos(t[j]), NewR[i]*sin(t[j]), Newincl[i]*cos(t[j]-Newphases[i])] end proc;

proc (i, j) options operator, arrow; [NewR[i]*cos(t[j]), NewR[i]*sin(t[j]), Newincl[i]*cos(t[j]-Newphases[i])] end proc

(3)

Surface := [seq([seq([NewR[i]*cos(t[j]), NewR[i]*sin(t[j]), Newincl[i]*cos(t[j]-Newphases[i])], i = 1 .. 51)], j = 1 .. 51)]:

plots[surfdata](Surface, labels = ["x [AU]", "y [AU]", "z [AU]"]);

 

``


Download Worksheet_Surfdata-contourplot.mwWorksheet_Surfdata-contourplot.mw

When trying to plot "Example1" from this page www.rhitt.com/courses/227/su99/Maple/html/surfaces2.html I get an error message that Maple cannot make a plot structure from an object with the transformation proc( (x,y) -> [x,y,-.3] end proc) . How can I plot this example?

Helo Dears, I need help to generate a vector whose detail given below

`for `K:=1 and m:=3 i need [[[psi[1,0]],[psi[1,1]],[psi[1,2]]]] and for K:=1 and m:=4 i need [[[psi[1,0]],[psi[1,1]],[psi[1,2]],[psi[1,3]]]]

 

`similarly for `K:=2 and m:=3 i need [[[psi[1,0]],[psi[1,1]],[psi[1,2]],[psi[2,0]],[psi[2,1]],[psi[2,2]]]] and for K:=1 and m:=4 i need [[[psi[1,0]],[psi[1,1]],[psi[1,2]],[psi[1,3]],[psi[2,0]],[psi[2,1]],[psi[2,2]],[psi[2,3]]]]

 

Simply I need to generate a general vector for arbitrary values of k and M as discus in above.

New_Wavelet_with_Operational_Matrix.mw

We're starting on indefinite integrals in my 1st year calculus class.

 

A quick example would be int(sin(x), x);=-cos(x)+C

 

Maple doesn't add the +C on the end of it's solution. Can someone explain or point me to a resource? I've tried searching but I can't find an answer. 

I have the following functions: and 

I want to differentiate e(a,A) wrt A and I keep getting the following error

Error, invalid input: ln expects its 1st argument, x, to be of type algebraic, but received [gamma/(c(A)^(sigma*phi)-1+gamma)].

Not sure how to proceed.

Thanks in advance for your help.

 

Hello there,

I have, of late, experienced some funny output when using the DifferentialGeometry package.

The command DGsetup retuns vector fields and 1-forms with stars enclosed in back primes attached to them. The same predicament occurs as a result of an operation (evalDG) but only when the vector or the form is not multiplied by a scalar.

This phenomenon does not occur with Maple 17!

Has anyone experienced such a phenomenon and how has they fixed the issue (providing it is an issue).

Many thanks in advance.

Kind regards,

Freddy Baudine

Here's the file:

DiffGeo_potential_issue.mw

First 1251 1252 1253 1254 1255 1256 1257 Last Page 1253 of 2224