How can I calculate GR tensors and geodesic equations for adS schwarzschild spacetime.

I've found the following project: http://www.parallella.org/

It is a very cheap but impressive computer ( 64-cores, they say it gives about 90 GFLOPS of computing power). The problem is the very limited amount of memory (1GB). See: http://www.parallella.org/board/ for specifications.

Now my question is: do you think Maple will run on this machine (acoording to the site it will run Ubuntu) and if so then does it make sense to try it given the small amount of memory it has? Or in another words: do there exist problems that could be solved by Maple on this powerful machine and that cannot be solved on a regular machine with let's say 4GB of RAM?

Hello everybody

I'm new at using Maple

so what I'm trying to do is " solve system of differential equations numerically " and plot the result 

I use the floweing code

PDEtools[declare]((u, v, w)(t), prime = t)

> params := z = 0;

Omega= 2.2758;

tau = 13.8;

T2 = 200; s = 1;

r = 0.7071;

\[CapitalDelta] = 1.7758;

s = 2.2758;

Eta= 1.05457173*10^-34;

omega = 0.5; k = 1666666.667;

> sys1 := {diff(u(t), t) = Omega*v(t)-u(t)/T2,

diff(v(t), t) = -Omega*u*{t}-2*s*exp(-r^2/omega0^2-t^2*1.177^2/tau^2)*cos(k*z-omega*t)*w(t)-v(t)/T2,

diff(w(t), t) = 2*s*exp(-r^2/omega0^2-t^2*1.177^2/tau^2)*cos(k*z-omega*t)*v(t)};

Cs1 := {u(-20) = 0, v(-20) = 0, w(-20) = -1}

> ans1 := dsolve*RealRange(Open({ICs1, sys1}), {u(t), v(t), w(t)});
%;
Error, (in RealRange) invalid arguments

plot([u(t),t=-20..20])
plot([v(t),t=-20..20])
plot([w(t),t=-20..20])

:::::::::

also I need to use the result of v(t) in another equation as,

x=2*v(t)*cos(k*z-omega*t)

How I can do that ?

I would like to define a graph in terms of its incidence matrix. It's very easy to go from the graph to the incidence matrix, but is it possible to go the other way? Here is the incidence matrix of interest:

M:=Matrix([[1,1,0,1,0,0,0],[0,1,1,0,1,0,0],[1,0,1,0,0,1,1],[0,0,0,1,1,1,0],[0,0,0,0,0,0,1]]):

Hello,

I was wondering if I can call Matlab R2012b with maple 14 on my macos 10.7.5.

When I try to do this:

> Matlab[setvar]("x", 3.14);

I get this:

Error, (in Matlab:-setvar) there was a problem finding or loading matlink.so. Refer to ?Matlab,setup for help configuring your system to work with the Matlab-link.

I read that I may have to change a script. Where are those scripts located?

Regards,

Three families grow vegetables in their backyards,and agree to form a small closed economy by sharing their produce. family andrews grow artichokes, family brown grows beans and family Cuthbert grows corn.Family Andrews receives 70% of its artichokes, 30% of the beans and 30% of the corn.Family Brown receives 20% of the artichokes,60% of the beans and 10% of thr corn.Family Cuthbert receives the remainder of the vegetables.

a) Write down the exchange matrix A.

b) If production is measured in doolars ,define the variables
                                    x[1]
,
                                    x[2]
,
                                    x[3]
 of the production vector X.

c) Find all solution of the Leontief closed with matrix equation (I-A)^-1¤X=0.

d) If family Cuthbert produces $100 worth of corn, how much will familes Andrews and Brown need to produce (in dollars) in order foe the economy to be in equilibrium?

Hey everyone,

I have seen similar questions a few times but I didn't really understand why this problem occurs in my system (as well as I am ney to maple). Therefore I would appreciate every help a lot!!!

I do have the following linear differential system:

As additional conditions I use:

And solve the system by doing:

additionally I define afterwards:

The result of for example Result[1] is:

Why do I get the unknown value of _U1? In my opinion I defined all parameters for the Linear differential system.

Thanks for every help!

Markus

hello 

please say me how can i make diff and implicit diff with Set Point answer?

for example how can i write these?

(x^2+y^2=6)' => 2*x+2*y*y'=0 => y'(-2,5)=2/5

or

(y=x^3+4*x)' => y'=3*x^2+4 => y'(3)=31

thanks very much

delete this post, not maple related

apologies

How would i label the axis and change its font and size.

Also, how would i move the horizontal axis down, so that the whole diagram can be seen. (Hope that makes sence)

Finally, there are some odd points between 0 and 1 on the vertical axis, how do you get rid of these...

Digits:=20: N:=10000: M:=100: x_max:=1: r_min:=0:
r_max:=4: for n from 0 to N do r:=r_min+n/N*(r_max-r_min):
x:=evalf(x_max*rand()/10^12):for m from 0 to M do x:=x*exp(r*(1 - x)): od:
X[n]:=x: od:
with(plots):
bifpoint:={seq([r_min+j/N*(r_max-r_min),X[j]],j=0..N)}:
pitchf:=pointplot(bifpoint,symbol=point):display(pitchf);

Thanks

I am trying to solve the rate equations for A<=>B<=>C with k1,k2,k3,k4 as the rate constants.  I put in the equations into Maple but can't seem to get it to solve them.  The initial conditions for A(0) = S, B(0)=C(0)=0

Here is what I plugged into the program:

Parameters(S, k1, k2, k3, k4);

sys_ode := d*A(t)/dt = -k1*A(t)+k2*B(t), d*B(t)/dt = k1*A(t)+(-k2-k3)*B(t)+k4*C(t), d*C(t)/dt = k3*B(t)-k4*C(t);

d A(t)
------ = -k1 A(t) + k2 B(t),
dt

d B(t)
------ = k1 A(t) + (-k2 - k3) B(t) + k4 C(t),
dt

d C(t)
------ = k3 B(t) - k4 C(t)
dt

ics := A(0) = S, B(0) = 0, C(0) = 0;
A(0) = S, B(0) = 0, C(0) = 0

dsolve({ics, sys_ode});
Error, (in dsolve) found functions with same name but depending on different arguments in the given DE system: . It is required an indication of the dependent variables

I have to questions. One is the main one, about simplifying fractions.

This is equations.

jj:=L*theta*gamma[2]^2*(-1+alpha)/(delta[2]*(-alpha*gamma[1]^(-2*sigma+2)*gamma[2]^(2*sigma)-gamma[2]^2+gamma[2]^2*alpha));

kk:=expand(denom(jj)/gamma[2]^2);

ll:=op(1,kk);

mm:=simplify(ll,power);

nn:=mm-delta[2]+delta[2]*alpha;

oo:=numer(jj)/gamma[2]^2/nn;

pp:=(numer(jj)/gamma[2]^2)/(denom(jj)/gamma[2]^2);

From now, I will explain. Though as I don't know how to put a maple screen on this board I think I will have you feel some inconvenience, please forgive it.

I want to make jj more simple, dividing both of numerator and denominator of jj, with gamma[2]^2.

But, what I could at best was this. First, taking denominator of jj and deviding it with gamma[2]^2 and expanding it like kk. Second, as the first term among that expression had terms of gamma[2] on both of the numer and denom, I picked out the first term as ll in the above maple commands, and simplified it with option: power, in mm. Thirdly, I rewrote kk as in nn.

and Forthly, I rewrote jj with modification of deviding numerator and denominator with gamma[2]^2, using the result which I could get as I wrote above. oo is that. And the result of oo is what I wanted to get.

But, on the other hand, when i divided denom and numer of jj directly by gamma[2]^2, the result I can get then is not 

good one: the output of pp is not what I wanted to get.

Is there more simple way to do what I did in the above? This is a first question.

And, another is how I can post the screen of maple on this question board?

I hope you will answer my question. I thank all of you who will answer these question in advance.

Hellp

for a task about the launch of a rocket we had to make a matrix OPL3 with in the first column the time T and in the second column the height of the rocket at that time.

I must Maple find the time where the rocket is at his highest point.

I can find how high the highest point is, but didn't succeed to find the corresponding time.
Can Maple tell me at wich time (= row # o/t matrix OPL3) the rocket is at his highest point.

I used the following to find the highest point.

ymin:=min(seq(abs(OPL3(i,2)),i=64..830));

Note: the diff equitation is numerically solved, so I can't use solve(....=0) (am I right here?)

Thanks!

How  can you create a loop for Monty Hall Problem when you have 3 door (1 opening) and then 4 Doors (with 2 openings and possible 2 switches)

I have this kind of problem. When I try to solve my command:

is there any way to avoid this? in order to get the answer directly...