I'm having a few problems with differentiating in Maple. I have a potential function U given by:

U[c] := (1/2)*r^2+M[1]/r[1]+M[2]/r[2];

U[r] := r^2*(M[1]*M[2]-3)/(2*c^2)+((x(tau)+diff(y(tau), tau))^2+(y(tau)-(diff(x(tau), tau)))^2)^2/(8*c^2)+3*(M[1]/r[1]+M[2]/r[2])*((x(tau)+diff(y(tau), tau))^2+(y(tau)-(diff(x(tau), tau)))^2)/(2*c^2)-(M[1]/r[1]+M[2]/r[2])^2/(2*c^2)-M[1]*M[2]*(1/r[1]+(1/r[1]-1/r[2])*(1-3*mu-7*x(tau)-8*(diff(y(tau), tau)))+y(tau)^2*(M[2]/r[1]^3+M[1]/r[2]^3))/(2*c^2);

Hi All,

Is using the NAG package the right way to go about creating the delaunay triangulation of a set of points in 3 dimensions. My points all lie on a surface in R^3.

Thanks. 

Hi All,

Recently I upgraded to Maple 14 from Maple 13, and I've encountered a strange behavior in Maple 14 that doesn't appear in Maple 13. By the definition of Spherical harmonic functions in Maple,  

SphericalY(2*l, -2*l, theta, phi) = SphericalY(2*l, 2*l, theta, phi).

And this is true in Maple 13, however not Maple 14. Can you reproduce this behavior in Maple 14? Or have I done something else that might be causing this?

Thank you for any help you can provide.

Hi there,

I have defined an own type MyType wich is internally just a list ['MyType', A, B] with let's say two lists A,B and the symbol 'MyType' as first element (to identify the type).

Now, if I define

X := MyType('MyType', A, B);

is there a possibility to extract A out of X by something like

X[2]

or even better X[1] (thus skipping the first one)? The only thing I know is op(2, X).

Hello,

I need to create a new bernoulli distribution such that P(I=1) = q = 1-P(I=-1).

I did U:=Distribution(PDF=(t->piecewise(t=-1,1-q,t=1,q,0))); but Maple gives me back:

module() option Distribution, Continuous; export...   instead of:

module() option Distribution, Discrete; export...  since the distribution is discrete!!

Do you know a way to force U to be a discrete distribution, cause otherwise I get all moments=0.

hi all

is there any (possibly simple) way to distribute a piece of maple software not allowing external user to access the source code ?

Translating to C, compiling and then including as an external function is very awkward procedure.

thanks in advance

Hello!

I'm trying to solve numerically an ODE system with piecewise. And this piecewise is very important for this task.

This system describes behavior of a pulley with friction. There are some constants: m, c, g, mu and J. Values of this constants are not important.

> sys := m*a(t) = piecewise(a(t) < a0, F0*time, a(t) >= a0, 0), v(t) = diff(x(t), t), a(t) = diff(v(t), t);
> m := 5; F0 := 10; a0 := 5;
> initialconditions := x(0) = 0, v(0) = 0;

Dear folks,

 In mapple programming while assigning a fraction to a variable...it directly simplifies it....

eg: x:=4/6 and it directly gives 2/3

can we make mapple not to simlify this fraction....

It is very important for me to find out....

please help me

thanks in advance

regards

uday

I have looked at the maple help for the ErrorPlot function, but i can't make it work for my plot. 

I have the following procedure:

K := 1.; Xrng := [seq(dt*i, i = 0 .. ntime)];

value := array(1 .. ntime+1);

ave := array(1 .. ntime+1);

for j to ntime+1 do:

discountfactor := exp(-r*j*dt);

Stot := 0.;

Stot1 := 0.;

for i to npath do:

Spath := max(Spot[i][j]-K, 0);

Stot := Stot+Spath;

What does it mean when I get the message, "Warning, model is not of full rank"?

Thank you.
Janice

Hi all,

Is there some way i can find the width of a plot?

I have the equation below:

Prob:=(w,T)->(B^2)*(sin(w*(T/2))^2)/((w/2)^2);

plot(Prob(w,6),wf=-3..3);

Which produces a gaussian looking plot, the plot reaches zero at some...

mem2.mw

Although i have defined K5 as a constant before solving the equations , maple returns an error message saying that

"Error, (in fsolve) k5 is in the equation, and is not solved for"

And i am alos trying t find the roots of the equations in a certain range as specified above.

Dear all,

Find here a part of my precedent post "2D finite element method" which represent a worksheet for 2D triangular Mesh.A part of the mesh can rotate with regard to the other by modifying theta.

restart: with(geometry): with(plots):

ms := 8: n := 6: mr := 5: theta := (1/180)*(23*Pi):

wrs := (1/1000)*<30.8, 33, 36, 39, 40, 43, 46, 49, 53.25>:

wrr := (1/1000)*<20,23,26,28,30,30.8>:

wt := (1/180)*Pi*<0,10,20,30,40,50,60>:

Dear All!

Could you explain how the coulditbe() function works?

> assume(q, integer, w, integer); about(q, w); coulditbe(q/w < 1); coulditbe(q < w);
Originally q, renamed q~:
  is assumed to be: integer

Originally w, renamed w~:
  is assumed to be: integer

                              FAIL
                              true

BR,

Zoltán Faigl