Given expr:=sin(n*x); and say I want to find its value at n=1,x=Pi/2, I found I can do

subs({n=1,x=Pi/2},expr); eval(%);

or combine into one

eval(subs({n=1,x=Pi/2},expr); ----(1)

But also one can write

eval(expr,{n=1,x=Pi/2}); ----- (2)

So if someone wants to just find value of expression as above, it seems (2) is much simpler than (1) since no need to call subs, it is done automatically. Is there anything else I am overlooking?

Only problem I see in the above, is that subs() uses different order: subs({},expr), while eval is:  eval(expr,{}); The other way around.

It would have been better if the order was the same in the API, to reduce confusion.

need to find steady state of this system ( finding this simultaneously) in maple 

How can you do it? 

I am trying to find the root of an equation that involves a procedure and a definite integral (solved numerically). Of course, I don't need the root to be found symbolically, but numerically would be fine. The problem is, I keep getting the error

"Error, (in fsolve) Can't handle expressions with typed procedures"

whenever I try to solve it. Anyone have any ideas? My worksheet is here: Table-1-duplication-mapleprimes.mw 

Hi,

I'm trying to fit some parameters to the data of multiple slightly different experiments.

I've written a function which returns the sums of the error for all experiments.

The function which has to be minimized does the following:

1) Set parameters for the model ODE based on the input parameters

2) calculate the difference (numeric solution of ODE <-> experimentel data)

3) repeat step 1+2 for all experiments

4) return the sum of all differences

The function works as aspected. But when I try to minimize it by calling:

I got an error:

"Error, (in dsolve/numeric/process_parameters) parameter values must evaluate to numeric, got A = A"

For some reason Maple isnt able to set the new parameters of the ODE.

Anyone got an idea how to fix this?

Thanks in advance.

I got a problem with solving a second order ODE. 

The ODE is :

-V(xi)+(1/2)*xi*(diff(V(xi), xi))+(1/4)*(diff(V(xi), xi, xi))=-(1/2)*k2*(diff(H(xi), xi))-k1*n*X/E+1+k2

where k1,k2,n,X,E  all are constant.

the condition is :

V(xi) tends to 2*xi^2 as xi tends to infinity.

I used 'dsolve' to solve the equation firstly, and got a solution with two constant C1 and C2, I want to use the condition to elimilate C2, so I used limit(sol,xi=infinity)=2*xi^2. But when I used the command 'limit', I can't get the answer.

Could any one help me? 

Many thanks!!!

I am using Maple 18 on a Mac. I also uploaded the tetrad package.  I setup my metric but when I call Christoffel I get the following error.

Download Equation_2.16_using_DiffernetialGeometry.mw

I appreciate the help.

Hi,

I'm a newer of maple .I want to prove x=z,y=z from equations of parametric z.But maple don't solve what I want,who can give me a help .thank you

sometimes when I download something from Maple apps centers, such as this package, http://www.maplesoft.com/applications/view.aspx?SID=33406  I find the example documents there are written in .mw and when I open them, there are in document mode style, the fancy word like style which I can't stand looking at.  It has all the math input in italic and maple commands look different from classic Maple text. It very confusing, since I see something as   (0<x<1) which is valid in this document mode, but in classical maple, this is not valid code.

It seems Maple has 2 different syntax. One that works in document mode and one that is classical text maple.

I like to use only input as Maple notation which these documents do not do.

Is there a way to convert such .mw file to become standard classic worksheet mode? I know I will lose the chapter/section heading and all that.  I also tried selecting all in the document, then did Tools->options->Display->Input display->Maple notation, but nothing happend to the open document. It remained document mode with 2D math input.

When I save it as .mws, and open it again, it remains on document mode.

How to convert such documents to classical Maple syntax?

Short question version:

print(plots[display](seq(seq...   ...)):  works fine!  Why??

Longer, more detailed question:
  print(plots[display](seq(seq...   ...)):   This Maple statement ends in a colon - so I would not have expected any output.  I'm using "pretty" print, as opposed to printf, because it automatically centers the output

    This question is actually an extension of one I posed  yesterday, under the title "How to clear Maple output" - but the question I should have been asking is: "Why is the previous output not automatically cleared?"  The main points of that question are stated below.  Any comments, suggestions gratefully received.

   David

"How to clear Maple output"

I'd like to know how the screen previous ouput can be cleared.   ...or what is the likely cause.

   At the start of my program I have:

restart;
interface(echo=0, verboseproc=0, warnlevel=0, prettyprint=1):

#I thought restart; would automatically have cleared any previous output.

The program reads necessary input (eg no of players, ..)  using:

x:=readline(terminal):

I would like to look at a 3d plot including an condition about the two variables of the plot

For this simple case, I only want to see the plot with the condition that x>0. Is that possible?

plot3d(x*y,[x,0,1],[y,0,1]);

Hi all,

which alternative options fo I have to prove the equality of two algebraic expressions if testeq fails?

The case im reffering to can be seen in the following document:

https://dl.dropboxusercontent.com/u/29147149/Exam%202006%20Question%20A.mw

Equation 14-18

I am pretty sure that the expressions are equal. evalb just returns false because it does not simplify expressions.

howdy.

evalf(map(log10,[25,5,1,10,4,20]))=

[1.39794000867203, .698970004336017, 0., 1., .602059991327960, 1.30102999566398]

How do I reverse this process , ie get back [25,5,1,10,4,20] using map command. Obviously it's 10 to the power but map(10^,[1.39794000867203, .698970004336017, 0., 1., .602059991327960, 1.30102999566398]) won't work

I'd like to plot the following inequalities:

sqrt(x)<=1/sqrt(2)

1/sqrt(2)<sqrt(x)<=1/sqrt(2)

Hey!

I have this MATLAB script, but as I work in Maple, I'll need to translate it to Maple. 
I know how to define symbols and functions, but I don't know which Maple commands to use, of if it needs to be done in another way, so the problem is translating from " dZ = [dx; dy; ax; ay; dm];" and down.

Any help would truly be greatly appreciated! Thank you!  

This is the script:

function dZ = meteor_step(~, Z)

  P = 1.2; % initial atmospheric pressure
  H = 1.39E-4; % scale height of atm pressure
  E = 8.11131859E6; % evaporation energy
  D = 1; % drag constant
  G = 9.814; % acceleration due to gravity
  PM = 3.3E3; % density of the meteor
  
  S = 3.986E14; % standard gravitational parameter of Earth (G*M)
  R = 6.3674447E6; % radius of the Earth (meters)

  x = Z(1);  
  y = Z(2);
  dx = Z(3);
  dy = Z(4);
  m = Z(5);
  
  atm = P*exp(-y*H);

  v = sqrt(dx^2+dy^2);

  area = pi * ( (3*m)/(4*PM) )^(2/3);

  dist = sqrt(x^2+y^2);
  Gv = -9.8;

  accel = -(D*atm*area)/m*v;
  ax = accel * dx;
  ay = accel * dy + Gv;
  

  dm = -(atm*v^3*area)/(2*E);
  
  dZ = [dx; dy; ax; ay; dm];
  
end

[t, R] = ode45(@meteor_step, [0 250], [0, 100000, 100, -300, 25]);

x = R(:,1);
y = R(:,2);

dx = R(:,3);
dy = R(:,4);

v = (dx.^2+dy.^2).^(1/2);

m = R(:,5);

figure(1);
plot(t, y);
  title('Meteor Kinematics: Height vs Time');
  xlabel('Time elapsed (s)');
  ylabel('Height (m)');

figure(2);
plot(x, y);
  title('Meteor Kinematics: Horizontal vs Vertical Position');
  xlabel('Horizontal (m)');
  ylabel('Vertical (m)');

figure(3);
plot(t, v);
  title('Meteor Kinematics: Speed vs Time');
  xlabel('Time elapsed (s)');
  ylabel('Absolute speed (m/s)');

figure(4);
plot(t, m);
  title('Meteor Kinematics: Mass vs Time');
  xlabel('Time elapsed (s)');
  ylabel('Mass (kg)');

  
figure(5);
plot(t, dy);
  title('Vertical Velocity vs Time');
  xlabel('Time elapsed (s)');
  ylabel('Vertical velocity (m/s)');

temp = abs(y - 52900);
[~, index] = min(temp);
  
t(index)
dx(index)/1000
dy(index)/1000

I've got a worksheet in which I have invested many hours of CPU execution time and if the computer goes down, or Maple fails for some reason, I'll lose it all.

If this happens I would like to be able to continue the calculation from where I left off.

By saving the worksheet periodically, I can save all the commands, but not the results, so if I have to restart, I'll have to wait many hours before the worksheet catches up to where it left off.

In ancient versions of Maple, you used to be able to save an executed worksheet including results (I sort of remember that you wrote a file with the extension ".M")  but the new help pages say that is now different from what it used to be. Obviously I can "save" individual symbols to a file, but for a complicated worksheet  that gets complicated.

I've read about "maplet" files, but that doesn't seem to fill the bill either.

So, is there any simple way to save a worksheet so you can continue seamlessly from where you left off after a crash, with all the previous results intact?