## Maple 13 Questions and Posts

##### These are Posts and Questions associated with the product, Maple 13

Some time ago, I have used the Virtual 3D Solar System code to plot a few interesting figures with asteroids. It can be found here: http://www.maplesoft.com/applications/view.aspx?SID=4484&view=html . I plotted the asteroid belt as a closed, warped surface according to how the asteroids are inclined, with the arguments of perihelion and longitude of ascending node where they cross the equatorial plane to the other celestial hemisphere. The plot was only done with axes in units of distances, for convenience's sake in astronomical units of 1 AU = 150 million km.

But I am strongly wondering if it were somehow possible to alter or adjust the code such that there is a zylinder placed around such plot where one can see celestial coordinates (Rectascension, Declination) on the inside, so that if one wants to know what coordinates an orbit of a specific asteroid has that one can look in the chart and go out on the balcony and set the telescope to these coordinates. Of course the asteroid most likely is not there on his entire orbital path, but you have at least the location where the asteroid CAN be. Would this be possible to plot?

This probably doesn't work if a closed cylinder is placed around the Virtual Solar System. But how about scrolling the cylinder and the orbits flat on a 2D plot?

I need to calculate Weyl scalar for a metric using null tetrad using debever package in maple 13. However, I am stuck at the defination of h (representing the covariant complex null tetrad). Is it the product of covariant null tetrad? I have worked it out by using covariant, contravariant and both covariant and contravariant null tetrad (like l_a*l_a, l^a*l^a, l_a*l^a), however, i am not getting the right result not even for the example given in the maple help for plane wave. please help me out that how should i define this h.

thanks,
suresh

## The color option of plot...

Hi Maple People

# Some Maple code
restart
x:= Vector(10):
y:= Vector(10):

for z from -5 to 4 do
x[z+6]:=z^2 + 40:
y[z+6]:=z^2 + z + 41:
end do:

plot(x,y,style=point,symbol=asterisk)

Regards

Matt

## How to convert expression to sum?...

I get the maple result as a product and sum off terms

x= (a+b*c)/d

i want to convert it to

x= a/d+b*c/d

i try to use convert(x,?)

## How to optimize this code?...

a1:= f(x) :
> T1 :=simplify((taylor(a1,x=alpha,N+3))):
> E1:=subs([seq(((D@@i)(f))(alpha) = 0,i=1..m-1),f(alpha)=0,x=e[n]+alpha],T1):
> g1 :=(convert(simplify(series((E1,e[n]=0,N))),polynom));

## How to express this constraint into Maple?...

Hi everyone...!

Can somebody tell me how to express this equation in Maple?

xij <= zkl ; ∀ i ∈ I: S(i)=k, ∀ j ∈ B: R(j)=l;

Currently I'm dealing with containerization problem and have 4 indexes in the constraints (namely: i for item, j for container, k for shipment, l for route, S for Set of Shipment, and R for Set of Route) while x and z are binary variables. What I want to express is: (for example), item 1,2,3 are in shipment 1, item 4,5 are in shipment 2, etc etc. SO, if i = 1,2,3 then the value of k will be 1. If i = 4,5 then the value of k will be 2, etc. Same thing goes to j and l, (for example) if j = 1,2 then the value of l will be 1, etc etc. Further depcition is more or less like this:

S(i) = k

S(1) = 1

S(2) = 1

S(3) = 1

S(4) = 2

S(5) = 2

Thank you very much for the help.

## Two independent variables...

HI, I am trying to solve two PDEs but in boundry conditions there is arising an error plz help.
Nazi.mw

## lambda for higher order ODE...

hi, I want calculate lambda(lagrange multiplier) for lane emden type equation, I have made a programme for it but here is error at the end of programme,plz help

lambda.3rd.mw

## How to solve those ODEs ?...

plz help me, how do i solve singular ODEs of lane Emden type equation for homotopy analysis method in maple? there is arising an arror, invalid fraction

## DELAY DIFFERENTIAL EQN 2...

please write a program on matlab or maple for this dde

y'(x)=-y(x-pi/2)       0<x<=10

y(x)=sin(x)       x<=0

tau=pi/2

## Length of output exceed limit 1000000...

can anybody help me? i want to check the consistency of my scheme. My equation is too long if i check manually, so i used maple 13 to simplify my equation. But it cannot simplify it because of length of output exceed limit 1000000

 (1)

 (2)

.

## how to use mathematica code convert to maple code...

I find mapleform software   but I dont know how to use , and does anyone know another method to convert.  I know use

> with(MmaTranslator);
> MmaToMaple();

I can  automatically translated .nb files.  but this have

Error, (in readline) file or directory does not exist
Error, (in readline) file or directory does not exist
Error, (in readline) file or directory does not exist
Error, (in readline) file or directory does not exist

another Error, missing operator or `;`

also I try to use

> with(MmaTranslator);
> FromMmaNotebook(Mma_notebook_filename, options);

I still dont know how is works?  can you explains for me and show me some example  , here is my example files.plot.zip

## There is no Startup Code Menu Item in Maple 13...

There is no menu item called "Startup Code"

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

Thanks!

## Opening a maple file...

Hi, there

I work with maple 13 and made a file that its contet is plotting curves and its size is

151 MB.Now when I want to open it the file cannot be opened completely only about a quarter of its contents rises.

REGARDS

Yegan