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 ?

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?

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

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

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...

I encounter "insufficient  initial/boundary value" error message,  do know how to proceed from there, search with "insufficient initial value" gets no result. Any help will be appreciated.

> restart; alias(r = r(t), f = f(t)); with(plots);
r, f
> DE := diff(r, t) = 2*r+alpha*r*f, diff(f, t) = -f+alpha*r*f;
d d
--- r = 2 r + alpha r f, --- f = -f + alpha r f
dt dt
> NULL;
> params := alpha = .3;
alpha = 0.3
> initv := r(0) = 101, f(0) = 2;
> NULL;
>
> dvars := [r, t];
> chaodisplay := proc (chartname) EQ := [op(subs(params, [DE])), initv]; EQ1 := dsolve(EQ, numeric); odeplot(EQ1, dvars, t = 0 .. 300, axes = frame, numpoints = 50000, color = green, orientation = [-30, 100], title = chartname) end proc;
Warning, `EQ` is implicitly declared local to procedure `chaodisplay`
Warning, `EQ1` is implicitly declared local to procedure `chaodisplay`
>
> chaodisplay("Rabbit and Fox");
Error, (in dsolve/numeric/type_check) insufficient initial/boundary value information for procedure defined problem
>

Dear people being here,

I have a question.

When I wrote as

aa:=(delta[1]*n[1])/(n[f1]*delta[1]);

maple returned as 

aa := n[1]/n[f1]

How can I have maple express the original one, aa:=(delta[1]*n[1])/(n[f1]*delta[1]), not what returned to me?

I hope you will answer to my question.

And, thanks in advance.

taro

Please help me, here my problem:

Parsing a string works correctely like this:

ABC := 3443;

print("ABC = ", eval(parse(cat("A", "B", "C"))));

But using PROCEDURE it doesn't work:

myParse := proc ()

     local PEA;
   

      PEA := 4334;

      print("PEA = ", eval(parse(cat("P", "E", "A"))))

end proc;

myParse();

                        "PEA = ", PEA

I suppose to have "PEA = ", 4334

Hi, I need to show some how by using Maple, that n^3>(n+1)^2 where n>=3. Iam new at Maple and not sure how I can graph inequality as n goes to infinity. Any help would be appriciated. Thanks!!!

Hi Everybody,

I'm having some problems simulating an algebraic system in Maple.  I want to do this for a larger system of equation (~20 equations/20 variables, with 3 inputs), but I'd like to get an idea of the form this should be in.

So, using an example out of Maple help:

with(DynamicSystems);

sys1:=AlgEquation([y=sin(x)],[x],[y],discrete=true,sampletime=0.001);

inp1:=Vector[column]([1,2,3,4,5,6,7]) * 0.1;

res:=Simulate(sys1,[inp1]);

Error, (in DynamicSystems:-Simulate) should never be here

 

I don't understand what this error means.  Reading the example in the help files, says the second argument should be a vector of input values.  What I would assume this would do is plugin x = 0.1, then solve y=sin(0.1), and output y.

Am I missing some syntax, or is my inputs wrong?  There are examples of how to set up an algebraic system in the help files, but no examples of how to simulate them.

Thanks in advance for any help!

Hello,
my question may be simple but I don't find the answer in any help guide.
when I define a function I cannot use a linearalgebra expression such as Trace.
Here is an example of what I would like to do:

If anyone can help me...
Thank you

I need to plot the band structure of Graphene using Maple, Ive written out the eqn as given in the homework, but when i hit enter it says "error. unable to match delimiters"

i have no idea what that means, here is the eqn i have typed into Maple:

http://imageshack.us/photo/my-images/39/j7ri.jpg/

http://img39.imageshack.us/img39/733/j7ri.jpg

i have been given values for gamma and acc, do i just type them in after, E is a function of k (kx,ky,kz) do i put that in brackets after the E?

E(kx,ky,kz)