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

how to make a table of newton raphson ? 
working is below 
printf("No           root                   x[i+1]                 abs(x[i+1]-x[i])\n"); for i to n do y[i] := evalf(x[i]+(-D1(x[i])+sqrt(D1(x[i])^2-2*f(x[i])*D2(x[i])))/D2(x[i])); x[i+1] := evalf(y[i]-f(y[i])/sqrt(D1(x[i])^2-2*f(x[i])*D2(x[i]))); printf("%d         %10.10f     %10.10f       %10.10e \n", i+1, x[i+1], f(x[i+1]), abs(x[i+1]-x[i])); if abs(x[i+1]-x[i]) < tol then print("approximate solution" = x[i+1]); print("No of iterations" = i+1); break end if end do
please help me 

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

I'm asked to solve the equations below simultaneously using fsolve.

f := x*exp((1 - x)*a);
f[2] := subs(x=f,f);
eq[1] := f[2] - x;
eq[2] := diff(f[2],x) + 1;

fsolve({eq[1],eq[2]},{x,a}):

fsolve doesn't seem to work here for some reason.

Any help help would be great!

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 all. I have a probably simple problem, that is drinving me totally crazy. I have a expression having:

 (sqrt(M+m)*sqrt(M-m).(sqrt(M+m)*sqrt(M-m))

and it refuses to simplify. I told maple assume(m>0, M>0), etc. But it doesn't want to simplify. Anyone has any ideas?

Regards,

Jelmew

Edit: uploaded the mw file, so it can be seen what the problem is.10.1.mw10.1.mw

[I added Physics to the tags.--Carl Love as moderator]

Hi, i'm trying to solve the simultaneous equations,

a[1]:=2*x^2 + 3*x^2 + x*y - x^2 + x;

a[2]:=3*y^2 + 4*x^2 - y;

eval(y,fsolve({a[1],a[2]},{x,y}));
0.

Even though y can be 0 it can also be 1/3 and two other complex numbers.

How do you get fsolve to show all four y solutions.

Secondly, how would i get maple to just to show the positive y solution ie. 1/3 only.

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)

Error, (in Statistics:-NonlinearFit) invalid input: no implementation of NonlinearFit matches the arguments in call

after tried many times, got this error

x11 := [0.208408965651696e-3, -0.157194487523421e-2, -0.294739401402979e-2, 0.788206708183853e-2, 0.499394753201753e-2, 0.191468321959759e-3, 0.504980449104750e-2, 0.222150494088535e-2, 0.132091821964287e-2, 0.161118434883258e-2, -0.281236534046873e-2, -0.398055875132037e-2, -0.111753680372819e-1, 0.588868146012489e-2, -0.354191562612469e-2, 0.984082837373291e-3, -0.116041186868374e-1, 0.603027845850267e-3, -0.448778128168742e-2, -0.127561485214862e-1, -0.412027655195339e-2, 0.379387381798949e-2, -0.602550446997765e-2, -0.605986284736216e-2, -0.751396992404410e-2, 0.633613424008655e-2, -0.677581832613623e-2]:
y11 := [ -21321.9719565717, 231.709204951251, 1527.92905167191, -32.8508507060675, 54.9408176234139, -99.4222178124229, -675.771433486265, 42.0838668074923, -12559.3183308951, 5.21412214166344*10^5, 1110.50031772203, 3.67149699000155, -108.543878970269, -8.48861069398811, -521.810552387313, 26.4792411876883, -8.32240296737599, -1085.40982521906, -44.1390030597906, -203.891397612798, -56.3746416571417, -218.205643256096, -178.991498697065, -42.2468018350386, .328546922634921, -1883.18308996621, 111.747881085748]:
z11 := [ 1549.88755331800, -329.861725802688, 8.54200301129155, -283.381775745327, -54.5469129127573, 1875.94875597129, -16.2230517860850, 6084.82381954832, 1146.15489803104, -456.460512914647, 104.533252701641, 16.3998365630734, 11.5710907832054, -175.370276462696, 33.8045539958636, 2029.50029336951, 1387.92643570857, 9.54717543291120, -1999.09590358328, 29.7628085078953, 2.58210333216737*10^6, 57.7969622731082, -6.42551196941394, -8549.23677077892, -49.0081775323244, -72.5156360537114, 183.539911458475]:
a1 := Diff(x1(t),t) = k1*x1(t)+ k2*y1(t)+ k3*z1(t);
b1 := Diff(y1(t),t) = k4*x1(t)+ k5*y1(t)+ k6*z1(t);
c1 := Diff(z1(t),t) = k7*x1(t)+ k8*y1(t)+ k9*z1(t);
ICS:=x1(1)=x11[1],y1(1)=y11[1],z1(1)=z11[1];
sol:=dsolve({a1,b1,c1,ICS}, numeric, method=rkf45, parameters=[k1,k2,k3,k4,k5,k6,k7,k8,k9]);
ans:=proc(p1,p2,p3) sol(parameters=[x1=p1,y1=p2,z1=p3]); end proc:
tim := [seq(n, n=1..27)];
ExperimentalData := <<x11>|<y11>|<z11>>;
FitParams:=Statistics:-NonlinearFit(ans, ExperimentalData,initialvalues=[k1=0,k2=0,k3=0,k4=0,k5=0,k6=0,k7=0,k8=0,k9=0], output=parametervalues);