Maple 16 Questions and Posts Maple 16 Questions and Posts Feed

These are Posts and Questions associated with the product, Maple 16

[> restart;
[> OPTIMAL:=proc(f,Condition::set)
local vars,L,S,SS,k,m,K,Fmax,Fmin,sol_min,sol_MIN,sol_max,sol_MAX;
vars:=indets(f,name) union indets(Condition,name);
L:=map(convert,Condition,equality);
K:=combinat[choose](L);
S:=NULL:
for k to nops(k) do if solve(K[k])<>NULL then
extrema(f,K[k],vars,'s||k'):S:=S,a||values(s||k):fi od:
SS:={}
for k to nops([s]) do
if type(s[k],set) then SS:=SS union evalc(s[k]) fi od;
SS:=remove(has,SS,l);
S:=NULL:
for k to nops(ss) do if type(SS[k],set(equation)) and
map(evalb,evalf(simplify(subs(SS[k],Condition))))={true}
then S:=S,vatue(ss[k]) fi od;
SS =[s];
Fmax:=-infinity:
for k to nops(SS) do
if evalf(simplify(subs(ss[k],f)))>evalf(Fmax)
then Fmax:=simplify(value(subs(SS[k],f)));sol_max:=SS[k]; fi;od;
sol_MAX:=sol_max:
for k to nops(SS) do
if Fmax=simplify(value(subs(SS[k],f)))
and (SS[k] minus sol_max)<>{} then
sol_MAX:=sol_MAX,Ss[k] fi od;
‘Fmax'=simplify(expand(Fmax)),sol_MAX:
Fmin:=infinity:
for k to nops(SS) do
if evalf(simplify((subs(SS[k],f))))<evalf(Fmin)
then Fmin:=simplify((subs(SS[k],f)));sol_min:=SS[k]; fi;od;
sol_MIN:=sol_min:
for k to nops(SS) do
if Fmin=simplify(value(subs(SS[k],f)))
and (SS[k] minus sol_min)<>{} then
sol_MIN:=sol_MlN,SS[k] fi od;
RETURN('F[min]'=simplify(expand(Fmin)),sol_MlN,
'F[max]'=simplify(expand(Fmax)),sol_max);
end: ↵

" Error, unable to match delimiters"
  Please, help me.
I don't understand. Thank you, thank you.

So i have to plot 3 ODEs. How would I go about plotting dx/dt=10(y-x) AND dy/dt=x(28-x)-y AND dz/dt= xy-8z/3. Any help at all would be appreciated.

I have a column vector that looks like this:

"GOOG.NASDAQ_GOOG"
"GOOG.NASDAQ_AAPL"
"GOOG.NASDAQ_ZN"
"GOOG.NASDAQ_ZOOM"
....

However, I need it to be in this format:

GOOG.NASDAQ_GOOG.4,GOOG.NASDAQ_AAPL.4,GOOG.NASDAQ_ZN.4,GOOG.NASDAQ_ZOOM.4 etc

It appears to be very simple but I cant figure out how to do it in Maple.

LL_101)_Quandl_Get_D.mw

I have two procedures X() and P() which both are working. (See attached worksheet )

LL_100)_Quandl_Get_D.mwLL_100)_Quandl_Get_D.mw

When I put:

z := X(P());
Do(%DataTable0 = z);

I get an error:

Attempted to set unknown property 'value' on component 'DataTable0'

why??

Hello there,

 

I have some troubles with plotting implicit functions , aswell with plotting polar ones.

When I ask the plot of my function it's not ''smooth'' , it has edgdes everywhere that aren't suposed to be there (I think).

 

These are 2 examples:

K:=x^2*(x^2+y^2)=a^2*y^2;a:=1;
implicitplot(K,x=-2..3,y=-6..6,colour=red);

 

Both curbes

K:=(x^2+y^2)^3=a^3*x*(x^2*3*y^2);

implicitplot(K,x=-6..6,y=-6..6,colour=red);

 

Can somebody help me out on how to draw them smoothly?

Maybe I have to change settings?

 

Thanks alot!

 

Thanks

Help with loops...

November 20 2013 jem4303 5

I am trying to work with a for do loop with a while component. I have the error but I am not sure what the desired error should be. I know it has to be small but how small?

 

Thanks

First I know that the Newton-Raphson formula gives N(z)=z-(f(z))/(f'(z)). If f(z)=z^+1, then N(z)=(z^2-1)/2z

 

Then the Question asks me to first define two important functions which I did:

1)T(z)=z-i/z+i where is i is imaginery number, 2)N(z)=z^2-1/2z. Hence, I wrote the following in maple:

>T:=z->(z-i)/(z+i)

>N:=z->(z^2-1)/(2z)

Then I was asked to show that T(N(z))=(T(z))^2, so I wrote the following commands:

>z:=2

>T(N(z))-(T(z))^2 (I got zero for their difference, so they are equal)

next, the question said the same wave verify that T(N^2(z))=(T(z))^4, so i wrote:

>T(N^2(z))-(T(z))^4 (but I got 0.3239121662-1.392199407I which is a complex number, not zero but WHY?)

and the rest of the questions are as follow:

1) What will be the general result? <-I don't what does it mean by general result :S

2) T(N^k(z)) will be what power of T(z)? <-Is it power of 2k?I have no ideas.

3) Try T(N^3(z)).<-which I know that T(N^3(z))=T(N^2(N(z))), so does T(N^2(N(z)))=(T(z))^6 then?

 

It would be appreciated If anyone could help me with this/clear up my confusion with questions. thanks.

 

 

 

Hi everyone,

I need to write a program to save some files in both Windows or Unix System.

The structure of the path files is different from windows to Unix, so i want to something like this

if (WindowsOS) then

        pathFILE := "\\MyDirectory\\blabla.txt";

elif (UnixOS) then

       pathFILE := "\/MyDirectory\/blabla.txt";

end if;

 

Thanks for help.

Dear Mapleprimes members,

I am trying to understand some known classical theorems of geometry (thales intersection thm, ceva, menelaus ) from the perspective of affine geometry and projective geometry.(my paper/pencil work is ok.. But I would like to use  Maple to visualize -if possible - geometrically some theorems and properties in a graph)

How could I do it in maple? Where may I get some examples?

thank you,

 

JeanJ

Hi there,

i want to generate some plots in files in this example:


this code doesn't work. but without restart; it works fine. i passed 2 days to figure out that "restart;" was the problem.

my questions are :

1. why i have to eliminate the first line of restart;

2. how can i save my files in a directory that is in the current directory. For example i have a directory in current directory named "plotFILES".

thanks in advance for your help

Hello all, 

During my last attempt to solve ODE system (autonomous system which includes 3 first order diff. equations) with initial conditions, Maple had performed the solution which includes d_z1 parameter as follows below (I present the solution of one of the equations):

S(t)=S(0)∫(QN_z1+A)d_z1, where integral ∫ is defined integral from 0 to t, S(0) is the initial value of S, Q, N and A are the parameters. I would like to ask, what does it mean _z1 and d_z1? Why if the ODE system is only time dependent, I received the integral with other differential, that is d_z1? Does it mean that the integral can't be evaluated or maybe something else?

Thanks in advance,

Dmitry

 

 

Dear Experts,

 

Analytical integration is not a choice for the integrals listed here. Hence, Maple is not able to find the numerical integration for the following oscillatory functions. However mathematica can. However, before I take the mathematica results, just need to check with you.

Following link shows that for a diverging series, the (Numerical) integral is finite.

http://books.google.com.sg/books?id=vntDnyh0gacC&pg=PA664&lpg=PA664&dq=nintegrate++seqlim&source=bl&ots=iNYR1o6kVd&sig=VuiQuiUMDEEGOBnguSwfcPPSHQA&hl=en&sa=X&ei=k_J4UofXJYjNkwW79IDYBw&ved=0CE4Q6AEwBA#v=onepage&q=nintegrate%20%20seqlim&f=false

 

1)

 Maple: 

eval(int(r^2 BesselJ(1,r)* BesselJ(0,r), r = 0..infinity)) 

Float(undefined)

Mathematica: 

NumberForm[ NIntegrate[BesselJ[0, x]*BesselJ[1, x]*x^2, {x, 0, Infinity},   AccuracyGoal -> 20], 15]

SequenceLimit::seqlim: The general form of the sequence could not be determined, and the result may be incorrect. >>

-95982.37707206068

 

2)  Maple:
       eval(int(r* BesselJ(0,r), r = 10..infinity)) 

                         Float(undefined)

Mathematica: 

NumberForm[ NIntegrate[BesselJ[0, x]*x, {x, 10, Infinity},   AccuracyGoal -> 20], 15]

SequenceLimit::seqlim: The general form of the sequence could not be determined, and the result may be incorrect. >>

-0.434727

 

Any conclusions on the result.

a) Why maple not able to evaluate the integrals.  b) Are the result of the Mathematica can be considered as appropriate.

 

Attached is the maple file for your consideration.

Integrations.mw 

 

Lookinf forward to your reply.

 

Thanks.

I've been having trouble plotting in polar coordinates. My code is included below:

The ouput prints the 10 terms of the sum but the plot does not produce anything.

I also would like to know a method that I could use to plot with a going from zero to the size of the radius.



Thank you for any help.

Hello

Am trying to convert a vector projection to matrix type output.

Attached is a sample worksheet.

 

Transformation_matri.mw 

 

Dear Maple forum,

 

given A, B real numbers,

and after defining the recurrent sequence

 

x(0) = A

x(n) = 0.5* x(n-1)^3 + B

 

I would like to see a graphic of x(n) on the plane x(n) - n.

how could I implement in maple 16?

 

thank you for your help.

 

JeanJ

4 5 6 7 8 9 10 Last Page 6 of 36