MaplePrimes Questions

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

Hi all
how can I solve an equation? Pleas, help!

I have a differential equation with one unknown.
I first solve the differential equation, and then use the boundary condition to find his unknown, however, get the error:
Error, (in sol_2) parameter 'X' must be assigned a numeric value before obtaining a solution.

How can I find X?

 F:=1:
 l:=1:
 a:=0.1:
 n:=4.5:

 B:=1.47*10^(-11):
 J_nc:=evalf(2*a*int(y^(1/n+1),y=0..a/2)):
 M_f:=z->piecewise(`and`(z>=0,z<=l),F*l,`and`(z>l,z<=2*l),F*l-F*(z-l),`and`(z>2*l,z<=3*l),F*l-F*(z-l)+F*(z-2*l));
 M_ed:=z->z-3*l:
 M_x:=z->M_f(z)+X*M_ed(z);
 eq := diff(V(z), `$`(z, 2)) = B*(abs(M_x(z))/J_nc)^n*signum:cond := V(0) = 0, (D(V))(0) = 0;
 sol_2 := dsolve({cond, eq}, numeric);
 sol_3:=solve({sol_2(3*l)=0},{X});



 

Please see attached document.

 

The 2 bits of data I am working with are a Constraint Function and Objective Function:

 

 

I need some pointers on how to do this. I keep getting solutions of

and I am ont receving an x value for some reason and leaving me unable to continue with the problem.

If anyone can send me some code I would be grateful. If anyone would like to send the correct full working code so I can see how to do this and review it I would be grateful.

 

 Code I had so far was:

 

 


 

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

hi friends

i have a problem in maple with an error

dx := diff(x(t), t, t) = -G*Mz*x(t)/(x(t)^2+y(t)^2)^(3/2):

dy := diff(y(t), t, t) = -G*Mz*y(t)/(x(t)^2+y(t)^2)^(3/2):

 G := 6.67*10^(-11); Mz := 6*10^24:
 IniC := x(0) = 7*10^6, (D(x))(0) = 0, y(0) = 0, (D(y))(0) = 9*10^3:
 Digits := 15:
 Ns := dsolve({dx, dy, IniC}, {x(t), y(t)}, numeric):
 dsnumsort(Ns(0), [x, y]);

>for i from 0 to 400 do T := 40*i; NsT := Ns(T); X[i] := rhs(NsT[C1]); Vx[i] := rhs(NsT[V1]); Y[i] := rhs(NsT[C2]); Vy[i] := rhs(NsT[V2]); MofI[i] := X[i]*Vy[i]-Y[i]*Vx[i] end do:

> with(plots);
> p1 := polarplot(6378*10^3, phi = 0 .. 2*Pi);
> p2 := plot([seq([X[i], Y[i]], i = 0 .. 327)], thickness = 2);

display({p1, p2}, labels = ['x', 'y'], scaling = constrained);

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

Error, invalid input: rhs received [t = HFloat(0.0), x(t) = HFloat(1.0), diff(x(t), t) = HFloat(0.0), y(t) = HFloat(0.0), diff(y(t), t) = HFloat(1.0), z(t) = HFloat(0.75), diff(z(t), t) = HFloat(0.0)][C1], which is not valid for its 1st argument, expr

 

can you helpe me?

 

I'm a student using Maple 18.

I'm still having some navigation issues with the Maple text editor. I realize that being in "math" mode or "text" mode causes Maple to interpret an <enter> key press differently.

However when I paste in some text, I get unexpected behavior. 

My professor will assign a list of problems, and I will paste them into Maple so it looks similar to this

1.) solve math problem 1

2.) solve math problem 2

....

 

Then I go to problem one, put my cursor at the end of the line and push <enter> hoping to make some room between problem 1 and 2 for my work. When I do this, the cursor jumps to the end of the list I had just pasted into the document.

How can I insert new blank space between the listed problems? I've had the same experience pasting in from MS word or just plain notepad.

Right now the only solution I've found is to copy the listed problems into maple 1 at a time.

 

Thanks in advance for any suggestions.

 

Brent

First 1294 1295 1296 1297 1298 1299 1300 Last Page 1296 of 2434