Unanswered Questions

This page lists MaplePrimes questions that have not yet received an answer

Hi there,

I would like to have an operator (in this case, the natural logarithm) applied to a list/array of points defined as:

ydata := [0.572594976618e-1, 0.327865007249e-1, 0.280821589546e-1, 0.114365745192e-1, 0.578537931608e-2, 0.139154661062e-2, 0.641467839994e-3, 0.18013801847e-3];

How can I apply Maple's ln() operator to the whole array (i.e. avoid to apply it to ydata [1], ydata [2], etc.)?

Thank you,

jon

 

I need to change the color of all elements of a plot by command.

I know that it is possible in interactive mode (see online maple help) but I would like to select the elements of a plot by line command and change color.

I would at first produce the plot with some color and then change only color of elements without reproduce the whole plot.

How can I do it ?

Thank you.

Hello friends!

I 'm a student and I don't know a lot about Maple, so I would be really grateful if anyone could help me.

I want to solve a system of two equations and I have two unknowns, which are k and εα. However I don't know what I am doing wrong and I can't solve it.

I have attached my file.

Thank you very much in advance!agogos_2.mw

 


Dear people in Mapleprimes,

 

I have a question about the ordering of monomials in a polynomial.

I hope you will help me understand how Maple works about it.

I inputed the polynomial as is written in black below.

Then, the outcome was blue, which ordering I could understand well: total degree ordering where at first 

those who have the order of 6 are collected which are 14 x^3*y^3, 6x*y^5, and then the following was those which 

have the order of 5: 21*x^5, -35 x^4*y, 9*x^3*y^2,-15*x^2*y^3, ... and so on.

And, among those who have the same order, lexical ordering was done, that is among 14 x^3*y^3, 6x*y^5, one which 

came first was the one with the larger degree about x, and among 21*x^5, -35 x^4*y, 9*x^3*y^2,-15*x^2*y^3, 

the first was 21*x^5, the second was -35*x^4*y, and so one, which was the ordering following the exponent about x.

 

And, then, I calculated Factor(polynomial) mod 7, which meaning I know.

Then, the result was 2*(x*y+2)*(3*y^3+x^2+3x*y)y.

I can understand the ordering among x*y and 2 in x*y+2, and that among 3y^3, x^2 and 3x*y in 3y^3+x^2*3x*y.

But, I can't understand why (x*y+2) comes at the first term, with 3 y^3+x^2+3x*y following it, and with y coming last.

 

This might be a trivial question. But, I hope you will teach me about this.

 

Best wishes.

 

taro 

 

polynomial := 14*x^3*y^3+6*x*y^5+21*x^5-35*x^4*y+9*x^3*y^2-15*x^2*y^3+12*y^4+18*x^2*y-30*x*y^2

14*x^3*y^3+6*x*y^5+21*x^5-35*x^4*y+9*x^3*y^2-15*x^2*y^3+12*y^4+18*x^2*y-30*x*y^2

(1)

`mod`(Factor(polynomial), 7)

2*(x*y+2)*(3*y^3+x^2+3*x*y)*y

(2)

NULL

``


Download heck_mod.mw

Hi,

  Is there any benchmark of the performance of maple on windows 7 64 bit vs linux, especially in solving generalized eigenvalue problem?

 

Thank you very much

is it possible?

i want to prove the injective of new custom operator based on custom logic table according to previous post

and want to see why new custom operator for integers are unordered or sometimes ordered

I want to solve maximize of equation,but the maximize failed to solve it,who can help me.thanks.

c[1] := (1/8)*w*{(1/((x+y+z)^2+1))^(3/2)+(1/((x+y)^2+1))^(3/2)+(1/((x+z)^2+1))^(3/2)+(1/((y+z)^2+1))^(3/2)+(1/(x^2+1))^(3/2)+(1/(y^2+1))^(3/2)+(1/(z^2+1))^(3/2)+1}+(1/8)*{x/((x+y+z)^2+1)+x/((x+y)^2+1)+x/((x+z)^2+1)+x/(x^2+1)}:

c[2] := (1/8)*w*{(1/((x+y+z)^2+1))^(3/2)+(1/((x+y)^2+1))^(3/2)+(1/((x+z)^2+1))^(3/2)+(1/((y+z)^2+1))^(3/2)+(1/(x^2+1))^(3/2)+(1/(y^2+1))^(3/2)+[1/(z^2+1)]^(3/2)+1}+(1/8)*{y/((x+y+z)^2+1)+y/((x+y)^2+1)+y/((y+z)^2+1)+y/(y^2+1)}:

t[1] := diff(c[1], x);

(1/8)*w*{-(3/2)*(1/((x+y+z)^2+1))^(1/2)*(2*x+2*y+2*z)/((x+y+z)^2+1)^2-(3/2)*(1/((x+y)^2+1))^(1/2)*(2*x+2*y)/((x+y)^2+1)^2-(3/2)*(1/((x+z)^2+1))^(1/2)*(2*x+2*z)/((x+z)^2+1)^2-3*(1/(x^2+1))^(1/2)*x/(x^2+1)^2}+(1/8)*{1/((x+y+z)^2+1)-x*(2*x+2*y+2*z)/((x+y+z)^2+1)^2+1/((x+y)^2+1)-x*(2*x+2*y)/((x+y)^2+1)^2+1/((x+z)^2+1)-x*(2*x+2*z)/((x+z)^2+1)^2+1/(x^2+1)-2*x^2/(x^2+1)^2}

(1)

t[2] := diff(c[2], y);

(1/8)*w*{-(3/2)*(1/((x+y+z)^2+1))^(1/2)*(2*x+2*y+2*z)/((x+y+z)^2+1)^2-(3/2)*(1/((x+y)^2+1))^(1/2)*(2*x+2*y)/((x+y)^2+1)^2-(3/2)*(1/((y+z)^2+1))^(1/2)*(2*y+2*z)/((y+z)^2+1)^2-3*(1/(y^2+1))^(1/2)*y/(y^2+1)^2}+(1/8)*{1/((x+y+z)^2+1)-y*(2*x+2*y+2*z)/((x+y+z)^2+1)^2+1/((x+y)^2+1)-y*(2*x+2*y)/((x+y)^2+1)^2+1/((y+z)^2+1)-y*(2*y+2*z)/((y+z)^2+1)^2+1/(y^2+1)-2*y^2/(y^2+1)^2}

(2)

eliminate({t[1], t[2]}, w);

[{w = -{1/((x+y+z)^2+1)-y*(2*x+2*y+2*z)/((x+y+z)^2+1)^2+1/((x+y)^2+1)-y*(2*x+2*y)/((x+y)^2+1)^2+1/((y+z)^2+1)-y*(2*y+2*z)/((y+z)^2+1)^2+1/(y^2+1)-2*y^2/(y^2+1)^2}/{-(3/2)*(1/(x^2+2*x*y+2*x*z+y^2+2*y*z+z^2+1))^(1/2)*(2*x+2*y+2*z)/((x+y+z)^2+1)^2-(3/2)*(1/(x^2+2*x*y+y^2+1))^(1/2)*(2*x+2*y)/((x+y)^2+1)^2-(3/2)*(1/(y^2+2*y*z+z^2+1))^(1/2)*(2*y+2*z)/((y+z)^2+1)^2-3*(1/(y^2+1))^(1/2)*y/(y^2+1)^2}}, {{1/((x+y+z)^2+1)-x*(2*x+2*y+2*z)/((x+y+z)^2+1)^2+1/((x+y)^2+1)-x*(2*x+2*y)/((x+y)^2+1)^2+1/((x+z)^2+1)-x*(2*x+2*z)/((x+z)^2+1)^2+1/(x^2+1)-2*x^2/(x^2+1)^2}*{-(3/2)*(1/(x^2+2*x*y+2*x*z+y^2+2*y*z+z^2+1))^(1/2)*(2*x+2*y+2*z)/((x+y+z)^2+1)^2-(3/2)*(1/(x^2+2*x*y+y^2+1))^(1/2)*(2*x+2*y)/((x+y)^2+1)^2-(3/2)*(1/(y^2+2*y*z+z^2+1))^(1/2)*(2*y+2*z)/((y+z)^2+1)^2-3*(1/(y^2+1))^(1/2)*y/(y^2+1)^2}-{1/((x+y+z)^2+1)-y*(2*x+2*y+2*z)/((x+y+z)^2+1)^2+1/((x+y)^2+1)-y*(2*x+2*y)/((x+y)^2+1)^2+1/((y+z)^2+1)-y*(2*y+2*z)/((y+z)^2+1)^2+1/(y^2+1)-2*y^2/(y^2+1)^2}*{-(3/2)*(1/(x^2+2*x*y+2*x*z+y^2+2*y*z+z^2+1))^(1/2)*(2*x+2*y+2*z)/((x+y+z)^2+1)^2-(3/2)*(1/(x^2+2*x*y+y^2+1))^(1/2)*(2*x+2*y)/((x+y)^2+1)^2-(3/2)*(1/(x^2+2*x*z+z^2+1))^(1/2)*(2*x+2*z)/((x+z)^2+1)^2-3*(1/(x^2+1))^(1/2)*x/(x^2+1)^2}}]

(3)

w = -(1/((x+y+z)^2+1)-y*(2*x+2*y+2*z)/((x+y+z)^2+1)^2+1/((x+y)^2+1)-y*(2*x+2*y)/((x+y)^2+1)^2+1/((y+z)^2+1)-y*(2*y+2*z)/((y+z)^2+1)^2+1/(y^2+1)-2*y^2/(y^2+1)^2)/(-(3/2)*sqrt(1/(x^2+2*x*y+2*x*z+y^2+2*y*z+z^2+1))*(2*x+2*y+2*z)/((x+y+z)^2+1)^2-(3/2)*sqrt(1/(x^2+2*x*y+y^2+1))*(2*x+2*y)/((x+y)^2+1)^2-(3/2)*sqrt(1/(y^2+2*y*z+z^2+1))*(2*y+2*z)/((y+z)^2+1)^2-3*sqrt(1/(y^2+1))*y/(y^2+1)^2);

w = -(1/((x+y+z)^2+1)-y*(2*x+2*y+2*z)/((x+y+z)^2+1)^2+1/((x+y)^2+1)-y*(2*x+2*y)/((x+y)^2+1)^2+1/((y+z)^2+1)-y*(2*y+2*z)/((y+z)^2+1)^2+1/(y^2+1)-2*y^2/(y^2+1)^2)/(-(3/2)*(1/(x^2+2*x*y+2*x*z+y^2+2*y*z+z^2+1))^(1/2)*(2*x+2*y+2*z)/((x+y+z)^2+1)^2-(3/2)*(1/(x^2+2*x*y+y^2+1))^(1/2)*(2*x+2*y)/((x+y)^2+1)^2-(3/2)*(1/(y^2+2*y*z+z^2+1))^(1/2)*(2*y+2*z)/((y+z)^2+1)^2-3*(1/(y^2+1))^(1/2)*y/(y^2+1)^2)

(4)

sub(w = -(1/((x+y+z)^2+1)-y*(2*x+2*y+2*z)/((x+y+z)^2+1)^2+1/((x+y)^2+1)-y*(2*x+2*y)/((x+y)^2+1)^2+1/((y+z)^2+1)-y*(2*y+2*z)/((y+z)^2+1)^2+1/(y^2+1)-2*y^2/(y^2+1)^2)/(-(3/2)*(1/(x^2+2*x*y+2*x*z+y^2+2*y*z+z^2+1))^(1/2)*(2*x+2*y+2*z)/((x+y+z)^2+1)^2-(3/2)*(1/(x^2+2*x*y+y^2+1))^(1/2)*(2*x+2*y)/((x+y)^2+1)^2-(3/2)*(1/(y^2+2*y*z+z^2+1))^(1/2)*(2*y+2*z)/((y+z)^2+1)^2-3*(1/(y^2+1))^(1/2)*y/(y^2+1)^2), c[1]);

sub(w = -(1/((x+y+z)^2+1)-y*(2*x+2*y+2*z)/((x+y+z)^2+1)^2+1/((x+y)^2+1)-y*(2*x+2*y)/((x+y)^2+1)^2+1/((y+z)^2+1)-y*(2*y+2*z)/((y+z)^2+1)^2+1/(y^2+1)-2*y^2/(y^2+1)^2)/(-(3/2)*(1/(x^2+2*x*y+2*x*z+y^2+2*y*z+z^2+1))^(1/2)*(2*x+2*y+2*z)/((x+y+z)^2+1)^2-(3/2)*(1/(x^2+2*x*y+y^2+1))^(1/2)*(2*x+2*y)/((x+y)^2+1)^2-(3/2)*(1/(y^2+2*y*z+z^2+1))^(1/2)*(2*y+2*z)/((y+z)^2+1)^2-3*(1/(y^2+1))^(1/2)*y/(y^2+1)^2), (1/8)*w*{(1/((x+y+z)^2+1))^(3/2)+(1/((x+y)^2+1))^(3/2)+(1/((x+z)^2+1))^(3/2)+(1/((y+z)^2+1))^(3/2)+(1/(x^2+1))^(3/2)+(1/(y^2+1))^(3/2)+(1/(z^2+1))^(3/2)+1}+(1/8)*{x/((x+y+z)^2+1)+x/((x+y)^2+1)+x/((x+z)^2+1)+x/(x^2+1)})

(5)

subs(w = -(1/((x+y+z)^2+1)-y*(2*x+2*y+2*z)/((x+y+z)^2+1)^2+1/((x+y)^2+1)-y*(2*x+2*y)/((x+y)^2+1)^2+1/((y+z)^2+1)-y*(2*y+2*z)/((y+z)^2+1)^2+1/(y^2+1)-2*y^2/(y^2+1)^2)/(-(3/2)*(1/(x^2+2*x*y+2*x*z+y^2+2*y*z+z^2+1))^(1/2)*(2*x+2*y+2*z)/((x+y+z)^2+1)^2-(3/2)*(1/(x^2+2*x*y+y^2+1))^(1/2)*(2*x+2*y)/((x+y)^2+1)^2-(3/2)*(1/(y^2+2*y*z+z^2+1))^(1/2)*(2*y+2*z)/((y+z)^2+1)^2-3*(1/(y^2+1))^(1/2)*y/(y^2+1)^2), c[2]);

-(1/8)*(1/((x+y+z)^2+1)-y*(2*x+2*y+2*z)/((x+y+z)^2+1)^2+1/((x+y)^2+1)-y*(2*x+2*y)/((x+y)^2+1)^2+1/((y+z)^2+1)-y*(2*y+2*z)/((y+z)^2+1)^2+1/(y^2+1)-2*y^2/(y^2+1)^2)*{(1/((x+y+z)^2+1))^(3/2)+(1/((x+y)^2+1))^(3/2)+(1/((x+z)^2+1))^(3/2)+(1/((y+z)^2+1))^(3/2)+(1/(x^2+1))^(3/2)+(1/(y^2+1))^(3/2)+[1/(z^2+1)]^(3/2)+1}/(-(3/2)*(1/(x^2+2*x*y+2*x*z+y^2+2*y*z+z^2+1))^(1/2)*(2*x+2*y+2*z)/((x+y+z)^2+1)^2-(3/2)*(1/(x^2+2*x*y+y^2+1))^(1/2)*(2*x+2*y)/((x+y)^2+1)^2-(3/2)*(1/(y^2+2*y*z+z^2+1))^(1/2)*(2*y+2*z)/((y+z)^2+1)^2-3*(1/(y^2+1))^(1/2)*y/(y^2+1)^2)+(1/8)*{y/((x+y+z)^2+1)+y/((x+y)^2+1)+y/((y+z)^2+1)+y/(y^2+1)}

(6)

"#"Iwant to maximize the equation (5)and (6),under the conditon of x,y,z are negative or positive at the same time.

 

NULL

 

Download maximize.mw

 

I recently downloaded a Maple reader software program from

http://www.crystaloffice.com/

I'mpuzzled as I'm not sure if this has anything to do with the Maple mathematics software,  On cursory observation it just looks like a text editor.

  I was curious to know if there was any software available which would allow people to read & execute a Maple program on the Internet. 

Cheers

   David

 

Dear All,

I'm trying to solve the following in Maple.

minimize(int(0.1e-3+.5*t+0.2e-2*t^2-b*t-a, t = 0 .. 300), location = true)

But Maple told me that the answer is

Float(-infinity), {[{a = Float(infinity), b = Float(infinity)}, Float(-infinity)]}.

I really need to get a kind of numerical answer. Would it be possible? Please Help me!!

Hi everybody,

is it possible to define an homokinetic joint in MapleSim for multibody modeling? How can I do that?

 

thanks.

In brief the problem can be stated as follows:

 

Given dependent variables Qi i=1,...,N and independent variables xi, yi, and zi i=1,...,N

which are related via the following system of N linear equations with parameters P1, P2 and P3 :

Qi = P1xi+P2yi+P3zi   i=1,...,N

How to find the optimal values of  P1, P2 and P3 which satisfy the above system of linear equations subject to the following constraints:

Pi>=0   i=1,2,3

and  P1>=P2P3

 Without the requirement of P1>=P2P3, the problem can be solved with the Non-negative Least Squares Method of Lawson and Hanson.  But with this additional constraint, I am stuck.  

 

Your suggestions are welcome.

 

 

 

 

I would like to implement a rule for typeseting delayDotProduct in latex exports (in which this expressions appears in plain form). It corresponds to dot products of vectors/matrices with the "." operator. I tried to define 

`latex/Typesetting:-delayDotProduct`:=proc(f1,f2,f3)
sprintf("{ {%s}{%s} }",cat(`latex/print`(f1)),cat(`latex/print`(f2)))
end proc:

which does not do the job. 

Edit:

A related question: How can I make

D[1,1](x)(t)

produce x''(t) in latex output?

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?

First 242 243 244 245 246 247 248 Last Page 244 of 363