Hello,

I would like to determine the position jacobian matrix from a set of constraint equations.

Here my constraint equations :

eq1:=l1*cos(theta(t))+l2*sin(beta(t))-x(t)=0
eq2:=l1*sin(theta(t))-l2*cos(beta(t))=0

The jacobian matrix that I would like to determine is :

Can you help me to make a general procedure to calculate a jacobian position matrix from a set of constraint equations ?

Thank you for your help

What is the easiest way to ask roots of a polynomial on a finite field. For example asking roots of x^2+xy+y on GF(8)? I was thinking to run a two for on members of GF(8) and ask to check it but I couldn't do it using Galois package or maybe I couldn't use that package. Thanks for any help.

I am wondering why Maple does this.

> f:=x^2
> f(3)
output: x(3)^2

I understand the difference between an expression and a function. If f is an expression, shouldn't it ouput  x^2(3). Why is the output x(3)^2?

When I enter f*3 or f*(3) then I get the correct expression.

Here is a screenshot. http://prntscr.com/a7u9hm

Here is image inserted with a slightly different function

Also while I am here, what exactly does g(x):= x^2 do? when i enter g(3) I get g(3) back.

g(x) is neither a function nor an expression.

screenshot http://prntscr.com/a7ua75

Dear All

I have third party Maple package saved along path E:/Maple work/General Maple Workout/TWS.mpl, but after using march command for other package, Maple is reporting error like "unable to read; E:/Maple work/General Maple Workout/TWS.mpl". What could be possible reason for this?

Moreover when I type "currentdir()" it shows me "C:\WINDOWS\system32" which is right path from where Maple is working

Regards

Hello,

I  put a caption under a plot and I want to use subscript. Since this is not the math mode, [] or __ does not work. Anyone knows, how to use subcript under the plot caption?

Thanks,

Hi all, 

I was wondering if it is possible to add a colour gradient to a point plot to represent another perameter in the data.

For example, each point corresponds to a value of 'f' ranging from 0.1 to 1, and I was wondering how to display this by means of a colour gradient.

Thanks,

 I am haunted by an equation which I could simplify it with hand and it equals to −i (complexe number). But, I don't know how to use maplesoft to simplify it. The equation takes the form of:

-Omega*a*sqrt(2)*sqrt(-Omega^2*a^2-2*k*m+sqrt(Omega^2*a^2*(Omega^2*a^2+4*k*m)))/(-Omega^2*a^2+sqrt(Omega^2*a^2*(Omega^2*a^2+4*k*m)))



I could simplify this equation with hand calculation and it equals to -i (the complexe number). I'am sorry for not clarifying that the a, k, Ω and m are positive variables.

Thank you in advance for taking a look.

This might be a mis-understanding on my part, so I figured I would ask a question first.  Narrowing down my code to something minimalistic, suppose I have 2 functions, each of which take a 'context' (abbreviated ctx) as a keyword parameter.  Now, suppose that the first one calls the second, like so:

foo := proc({ctx :: list := []}) bar(ctx = ctx) end proc:
bar := proc({ctx :: list := []}) nops(ctx) end proc:

and then a call "foo(ctx = [a,b,c])" returns the (completely unexpected) value 0.  See if you can puzzle out why.  If I change my code to use different names, like

foo1 := proc({ctx :: list := []}) bar1(_ctx = ctx) end proc:
bar1 := proc({_ctx :: list := []}) nops(_ctx) end proc:

Then the call "foo1(ctx = [a,b,c])" returns the expected 3.

I have tried a number of variants, like changing the call to bar('ctx' = ctx) in foo, but that does not work.  The completely un-intuitive :-ctx does work.

Is this documented anywhere?  Is this really the intended design?  Not being to re-use a name for a keyword parameter without going through some contortions seems, a little, shall I say, odd?

Just starting some work on numerical engineering solutions by going through a workbook, which includes a Maple Document on CD (with one index.mw and several .mws files).

As I don't have acces to a Maple system, I downloaded Maple Player. However - Maple Player will only accept .mw and .mwz files(at least that shows up on 'FILE>OPEN'); so trying to open an .mws file results in - NOTHING. In your forum I found some information that .mwz is a nonpublic Maple format !?...... Now I don't see a way to work with PLAYER. Isn't PLAYER meant for exactly that scenario - sharing MAPLE worksheets ?

So - what to do ?

Thanks for some hint,

GW.G

How do I write a procedure to find a root of f(x)=0 in the vicinity of a given value x0. The procedure should initially use the rearrangement method to produce a linearly convergent sequence of values, and should, when appropriate, switch to Aitkin's Method. The input for the procedure should be the re-arranged function and the velue for x0. The output should be the root and the number of iteration taken. The procedure should check that re-arrangement will converge. This program should do in Maple V Release5.

Thank you for your help.

Hey,

I wonder if there is an operator to skip inputs and outputs of functions. Consider this example:

The output would be:

However my system is much bigger and I only need the second output of the command and I dont want to waste memory on the first.

In Matlab I would write

~,ans := GenerateMatrix(...)Is there a similar shortcut availible in Maple?

Thanks in advance!

Hello,

I would like to simplify this following trigonometric expression :

I would like to make groups like : cos(a)cos(b) - sin(a)sin(b)=cos(a+b)  but keepind the maximum of expression products

On the following example (2 equations below), the function combine(expr,trig) works well 

But, I would like maple do only the first simplifications in order to the maximum of expression products. The function combine(expr,trig) goes too far in the first equation and I obtain only expression sums. 

Do you have ideas to simplify the first trigonometric equations
- with groups like : cos(a)cos(b) - sin(a)sin(b)=cos(a+b)
- and keeping products of expressions ?

Thank you for your help

f=x^2.  It's easy to compute the fouier transform. F(y)=fourier(f,x,y).

However, I want to do more. replace the variable y with a new formula g(z)=z^2.

I tried  "subs(y=g,F)". But failed. Need your help.

Download code.mw

I have:

r:=t-><4cos2t+4sint,5cos3t+4cost,(sin3t+2sint)^2>

with(plots):

then I tried:

spacecurve(r(t), t = 0 .. 2*Pi, thickness = 2, color = black, axes = normal, labels = [x, y, z], numpoints = 150)

But I keep getting the error: Warning, unable to evalute the function to numeric values in the region. 

But I thought <4cos2t+4sint, 5cos3t+4cost, (sin3t+2sint)^2> was defined on the region t=0 to 2Pi??? 

I'm using Maple 2015 if that helps. 

Dear all;

I open this good discussion, and hope can get a nice and strong idea in this domain of approximation of Hankel funciton and order truncation of infinite series. Thanks for all idea, can improve the discussion. 

Using the asymptotics of Hankel function for large argument and large orders ( both together) and
find   an order of truncation N of the obove series so that we can ensure an error bound  of epsilon( epsilon very small given).  abs(sum ( c[m]*HankelH1(m,x)*exp(I*m*theta), m=-infinity..infinity)-sum ( c[m]*HankelH1(m,x)*exp(I*m*theta), m=-N..N))<epsilon.

A first idea come in mind: the series converge, so that the general terms of this series converge to zero, and in particularity,  abs(c[m]* HankelH1(m,x))<1:
then abs(c[m])<1/ abs(HankelH1(m,x)).
So we can ensure an error bound of epsilon on the coefficient c[m] by imposing  (HankelH1(m,x))<epsilon  this lead to abs(c[m])<epsilon.
I consider the case where m and x are very large, I can suppose for example m=x*(1+zeta), with 0<zeta<1. zeta parameter. So that our truncation N depend on zeta.
And then how can I find and approximation of the inverse of Hankel function for large argument and large order. using m=x*(1+zeta). I think this give us N the truncation order.
I hope get a good discussion in this subjet.
Of course maybe there are other strong idea to compute the truncation series.
I get the following error in the code:

Error, (in MultiSeries:-multiseries) unable to expand with respect to parameter

###### Code### and error
restart:

with(MultiSeries):
assume(0 <= x);

assume(0 <= zeta<1);

HankelH1(v,x):

sum(c[m]*HankelH1(m, x)*exp(I*m*theta),m=-infinity..infinity);

1/HankelH1(x*(1+zeta),x);

MultiSeries:-asympt(%,x, 4);

eval(%, O=0);

convert(%,exp);

simplify(%);

truncated_series.mw