Advanced Search

Maple Questions and Posts

These are Posts and Questions associated with the product, Maple
View unanswered posts
Maple Questions and Posts Feed

Print "Hello" in Green?...

Asked by: mapleatha 160
print
April 07 2015
1  11

How do we print "Hello" in green in a .mw file in Maple 13?

Thanks!

mapleatha

Picking out Christoffel Symbols...

Asked by: shillis 20
April 07 2015
3  2

I would like to pick out specific christoffel symbols once maple calculates them. Right now I am using CoefficientList, but I would like a more direct way, as I am working with someone using Maple 17, and they are put in a different order in the CoefficientList. Here is the bit of my code I would like help with.

Error, invalid input: with expects...

Asked by: bfathi 10
April 07 2015
0  6

hi friends

After this cods i see very error

 > restart;

read(orbit.sav ): whit(plots):
ax := -G*Mz*x/(x^2+y^2)^(3/2);
ay := -G*Mz*y/(x^2+y^2)^(3/2);
i := 'i'; j := i+1;
for k from 0 to 3 do
x := 7*10^6; Vx := 0;
y := 0; Vy := 9000;
dt := evalf(1/2^k);
for i from 0 to 328 do
X[i] := evalf(x); Y[i] := evalf(y);
for n to 40*2^k do
x := evalf((1/2)*ax*dt^2+Vx*dt+x); y := evalf((1/2)*ay*dt^2+Vy*dt+y);
Vx := evalf(ax*dt+Vx); Vy := evalf(ay*dt+Vy)
od;
if i mod 41= 0 then
dX[k, i] := X[i]-XS[j]; dY[k, i] := Y[i]-YS[j]
fi
od;
p[k] := plot([seq([(X[i]-XS[j])*(1/1000), (Y[i]-YS[j])*(1/1000)], i = 0 .. 328)], color = green) end do;
p1 := display({seq(p[k], k = 0 .. 3)}, thickness = 3)
SI := [seq(41*i, i = 0 .. 8)]
p2 := plot({seq([seq([(1/1000)*dX[k, i], (1/1000)*dY[k, i]], k = 0 .. 1), [0, 0]], i = SI)}, color = black)
display({p1, p2}, scaling = constrained, labels = ['dx', 'dy'])
display({p1, p2}, view = [-.1 .. .5, -.4 .. .2], scaling = constrained, labels = ['dx', 'dy'])

can you help me Please?

Thank you

 

 

 

ColumnGraph in 3D...

Asked by: Alger 524
3d column plotting
April 07 2015
1  7

Dear all,

Is there a Maple command to get ColumnGraph in 3D as bar3 in Matlab.

In 2D, I use:

Statistics[ColumnGraph](A);  # equivalent to bar(A) in Matlab

Thanks

how do i solve this problem?...

Asked by: hooshmand 5
fsolve
April 07 2015
1  1

Hi,

we want to know what is the meaning of this statement?

We expect to have the following statement . But unfortunately we don,t get it

> with(difforms);
> sol := fsolve({diff(S, x) = 0, diff(S, y) = 0}, {x, y});


I do not take values above code.

Can't figure this out...

Asked by: Rouben Rostamian 8711
April 06 2015
1  7

Consider

> z := Int(f(t-s), s=0..1);

Int(f(t-s), s = 0 .. 1)

The result is clearly a function of t (and definitely not a function of s, which is a dummy variable).  Let's define

> F := unapply(z, t);

proc (t) options operator, arrow; Int(f(t-s), s = 0 .. 1) end proc

Then

> F(q);

Int(f(q-s), s = 0 .. 1)

as expected.  However:

> F(s);

Int(f(0), s = 0 .. 1)

which is not correct, as this confuses the argument s of F with the dummy integration variable s.  How would you salvage the situation?

 

Problem with integration....

Asked by: schzan 130
April 06 2015
0  11

I am performing the int operation on an experssion but it does not consider constants to be constants. For example, maple is considering this constant namely, f'(-1) or eval(diff(f(x), x), x=-1) as a function and henceforth does not give the desired output. Please tell me how to make eval(diff(f(x), x), x=-1) as a constant function or to say, how to convert any given function in to a constant function.

Thank you for your time and help!

Error, invalid left hand side in assignment...

Asked by: wenny 25
April 06 2015
1  2

Here is my code

ma[1] := [-885.880598855072776, [bh = 0., g0h = 0., g1h = 0.825946224722250e-4]]

ma[2] := [-877.957885609114328, [bh = 347.116836805625, g0h = 0., g1h = 0.164861392564e-3]]

ma[3] := [-863.445203144006655, [bh = 0., g0h = 0., g1h = 0.787090377403668e-4]]

avb := 90.5276611129000; avg0 := 0; avg1 := 0.92225359e-4;

for j from 1 to 3 do

assign(ma[j][2]);

A[j] := [(avb-bh)^2, (avg0-g0h)^2, (avg1-g1h)^2] 

end do;

 

Result is that

A[1] := [4763.19965962732, 0., 1.13103562500664*10^(-10)]

Error, invalid left hand side in assignment

 it only shows A[1], but no A[2] and A[3]

Please help!

Convolution animation...

Asked by: nMaple 85
animation
April 06 2015
1  21

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.

My math discovery done with Maple

by: Markiyan Hirnyk 7819 Maple
cases
April 06 2015
4  2

 

     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

how to using maple code convert to mathematica cod...

Asked by: candy898 60
convert mathematica
April 06 2015
0  1

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.

Installing with windows 8.1...

Asked by: pole99hode 5
windows
April 06 2015
1  1

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 ?

the series solution...

Asked by: Parham2016 70
dsolve series
April 06 2015
0  6

 

``

``

restart

``

Order := 6

ode := diff(f(eta), eta, eta)+(a/eta+(3*n+1)*`&Nu;u`*(1-(eta/a)^((n+1)/(2*n)))/(n+1))*f(eta)/(4*a*eta) = 0;

diff(diff(f(eta), eta), eta)+(1/4)*(a/eta+(3*n+1)*`&Nu;u`*(1-(eta/a)^((1/2)*(n+1)/n))/(n+1))*f(eta)/(a*eta) = 0

 (1)

ans := dsolve({ode}, f(eta), type = 'series')

Error, (in series/exp) unable to compute series

  

``

``

``

``

Order := 6;

6

 (2)

ode := diff(T(r), r, r)+(diff(T(r), r))/r+(3*n+1)*(1-r^((n+1)/n))*`&Nu;u`*T(r)/(n+1) = 0;

diff(diff(T(r), r), r)+(diff(T(r), r))/r+(3*n+1)*(1-r^((n+1)/n))*`&Nu;u`*T(r)/(n+1) = 0

 (3)

ans := dsolve({ode}, T(r), type = 'series')

NULL

``

 

Download series.mw

series.mw

 

A Simple Example of a maple.ini File and Its Locat...

Asked by: mapleatha 160
April 05 2015
1  3

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 Startup Code Menu Item in Maple 13...

Asked by: mapleatha 160
startup
April 05 2015
1  3

There is no menu item called "Startup Code"

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

Please help!

 

Thanks!

First 1282 1283 1284 1285 1286 1287 1288 Last Page 1284 of 2255