MaplePrimes Questions

Hi All,

I would like to plot3d a function say z=f(x,y). x=0..1, y=0..1.

But I only want to display a surface of the plots, say z=0.1. Is there any simple function to do that?

I am looking for soemthing like: plot3d(f(x,y)=0.1, x=0..1,y=0..1)

thanks.

Is it possible to split an execution group containing 2D input?
(without conversion to 1D which destroys the format).
Using F3 or the menu seems to work only for lines with a prompt.
But usually an execution group has a single prompt; lines with prompts appear e.g. when two execution groups are joined (with F4).
Is copy&paste the only solution?

Hi!

I am trying to export the plots I generated into black and white images. Currently, I am trying to use a proc as:

ExportFunctionPlot := proc(p::evaln, pname)
    local name, place, opts:
    name := cat(pname, ".eps"):
    opts := `landscape,width=768,height=768,noborder,axes=boxed,color="Black"`:
    plotsetup('eps', 'plotoutput'=name, 'plotoptions'=opts):
    print( plots:-display( eval(p), 'axesfont' = [ TIMES, 30 ],
                        'labelfont' = [ TIMES, ROMAN, 30] ) ):
    plotsetup(default):
end proc:

However, the eps files I got still have colored lines. Is there any way I can export to just black and white eps images?

Thanks!

William

Hi All,

I am thinking is there any simple function that could extract part of the points that have been ploted by plot() or plot3d(). And more ideally, we can specify some constraints such that the points being extracted must satisfy.

I guess Maple must stored the points somewhere in the memory in some kinds of data structure....

thanks.

I have a system with unknown functions: R(x), mu(x), Y(x). I need to find them AND also mu'(x), Y'(x), Y''(x). How do I write it in the code?

I tried something like this:

F := dsolve({cond, sys}, [R(x), mu(x), Y(x)], numeric, output = listprocedure);
Y1 := eval (Y(x), F);
Y2 := x -> diff(Y1(x), x);

It works, but does it work correctly? Y1 is a massive of numbers, is it correct to write "a derivative from a massive of numbers"?

 

I also tried this:

F := dsolve({cond, sys}, [R(x), mu(x), Y(x), diff(mu(x),x), diff(Y(x),x), diff(Y(x),x$2)], numeric, output = listprocedure);
Y1 := eval (Y(x), F);
Y2 := eval (diff(Y(x),x), F);

But it doesn't work at all and gives me an error on F: 

Error, (in dsolve/numeric/process_input) dependent variables must be functions of a single unknown, the independent variable. Got [diff(Y(x), x), diff(diff(Y(x), x), x), diff(mu(x), x)]


How should I really do it?
3.1.mw
3.2.mw

I am seeking the limit of a series representation of a waveform to within so many significant figures.  Unfortunately, as the series contains more & more terms the demands on computation memory become a problem.  Perhaps there is a method within MAPLE I am not familiar with that can possibly circumvent this problem?

Below is the link to my worksheet.  The presumed theoretical limit is supposed to be 0.0894945 (this is based on what I have read)  (perhaps it may be wrong).  So far I have 3 of the sig figs agree, but would like to take it further.  Is this possible with different MAPLE commands than what I am currently employing?  Currently, MAPLE crashes if I venture beyond 100000 terms in the series.

seeking_the_limit.mw

Hi All,

I would like to generate a sequence (say 10) of points <x_i,y_i>. (i=1, 2,..,10.)

There is a predefined constraint like f(x_i,y_i)=0.01. And I don't know the close form of f() (becasue there is a intermediate parameter in f() that was caculated numerically by optimization.).

Is there any simple function, e.g. a variant of seq(), that can do this? I am looking for something looks like:

seq(f(x,y)=0.01, x=0..1, y=0..1, 0.1)

Thanks.

Ok for some reason Get-Property doesnt work here but it does in another worksheet but anyway that problem isnt really what im posting for help for, i want general advice on the asthetic lay out of education tools.

 

The example i am doing tonight is pretty standard, i want to have a window where the student can play around with the recurrence relation making little alterations as he or she is curious about,  then hit a button to produce the first N values of the sequence, but see the problem i always get is that i dont care how messy or easy to use anything is because i am the one that wrote the code, so the features of such an analysis tool are redundant in the fact that i know how to use maple to a reasonable extent, but i want to make education  tools that are focused purely on the mathematics, removing the programming and all of that nonsense from the students focus here, so, i need advice for this given this stated objective.

 

Also forgot to mention, i already know how to get all the embedded components working together, this part isnt the problem i am referring to. I mean to say how can i best present the interface of these tools given what i can make in maple.

 


 

restart; with(DocumentTools); with(numtheory); with(StringTools); with(DocumentTools)

S[0]:   

s := proc (n) options operator, arrow; S[n-1]^2-2 end proc:

INIT0 := 4:

SOLVEfirstNterms(4)

[14, 194, 37634, 1416317954]

(1)

``


 

Download MAPLE_HELP_RECURRENCE_SEQUENCE_ANALYSIS_TOOL.mw

I was recently cleaning up a worksheet to make things more succinct.  In that process I modified how I expressed the series coefficients, Ck.  What seemed to be an innocous change apparently upset MAPLE to the point it cannot process my results.  In the abbreviated worksheet link below I process the results in the previous manner, Ck1, as opposed to my current modification, Ck2.

In the previous manner the results are generated under 3 minutes.  After modifying the expression for Ck, MAPLE cannot seem to process the results at all.  As far as I can tell the Ck1 & Ck2 concur.  So I am perplexed.  Can anyone see what is wrong?  The only thing I can think of is that sin(2*pi*k/T*x) in the denominator might cause the problem.  However, the sin term is cancelled out by the same sin term in S4.

Ck_modification.mw

Probably I am overthinking this but in any case I have hit a road block. If anyone can point me in the right direction I woould greatly appreciate it! Peter_Ursa.mw
 

restart

with(Student[MultivariateCalculus])

with(VectorCalculus) 

with(plots)

2. The surface of a mountain is modeled by the equation "h(x,y)=400-2 x^(2)-y^(2)."  A mountain climber is at the point 10, 10.  

"h(x,y):=400-2 x^(2)-y^(2)"

proc (x, y) options operator, arrow; VectorCalculus:-`+`(VectorCalculus:-`+`(400, VectorCalculus:-`-`(VectorCalculus:-`*`(2, x^2))), VectorCalculus:-`-`(y^2)) end proc

(1)

p := `<,>`(10, 10)

Vector[column](%id = 18446746427608139406)

(2)

a) In what direction should the climber move in order to ascend at the greatest rate?

 

 

H := Gradient(h(x, y), [x, y])

Vector[column](%id = 18446746427608141086)

(1.1)

eval(H, [x = 10, y = 10])

Vector[column](%id = 18446746427608133630)

(1.2)

``NULL

 

 

b) Determine a rectangular equation for the path of the mountain climber.

 

 

restart

  `<,>`(10, 10)+`<,>`(-40*t, -20*t)

Vector[column](%id = 18446746427574846638)

(2.1)

 

 

sys := {10-40*t = x, 10-20*t = y}

{10-40*t = x, 10-20*t = y}

(2.2)

solve(sys, {x, y})

{x = 10-40*t, y = 10-20*t}

(2.3)

t := solve({x = 10-40*t, y = 10-20*t}[1], t)

1/4-(1/40)*x

(2.4)

10-20*t

5+(1/2)*x

(2.5)

````

``

NULL

c) Plot the path of the climber along a 2-D contour diagram of the surface of the mountain.

``

 

 

restart

with(plots):

p1 := contourplot(-2*x^2-y^2+400, x = -15 .. 15, y = -15 .. 15)

p2 := plot(5+(1/2)*x, x = 0 .. 10, y = 5 .. .10)

display(p1, p2)

 

``

d) Show a 3-D view of the path of the climber along the surface of the mountain.

p3 := plot3d(-2*x^2-y^2+400, x = -15 .. 15, y = -15 .. 15)

p4 := plot3d(5+(1/2)*x)

display(p3, p4)

 

e) Animate the path along the surface of the mountain.NULL

``


 

Download Peter_Ursa.mw

 

Hello everyone,

  How do I go around this PDE. I think I have input all the necessary required but its given the error code "Error, invalid input: subs received [beta = .5, v = .5, v__l = 1.0, eta = 1.5, false], which is not valid for its 1st argument". 

Below is attached the code.

pde.mw

Please, anyone with any information should help. Thanking you in anticipation.

how i can pdsolve this equation

 

Doc186.pdf

I am working with the following differential equation:

$\frac{d^2z}{dx^2}+z=\frac{\cos 2x}{1+\epsilon z},\:\:\:z(-\pi/4)=z(\pi/4)=0$

where modulus of $\epsilon$ is much less than $1$.  The task is then to use perturbation theory (with Maple, if necessary) to show that the second-order approximation to the solution to this DE is:

$z=-\frac{1}{3}\cos 2x +\epsilon\bigg(\frac{1}{6}-\frac{8\sqrt{2}}{45}\cos x - \frac{1}{90}\cos 4x \bigg) + \epsilon^2 \bigg(\frac{2\sqrt{2}x}{45}\sin x - \frac{\sqrt{2}}{90}(\pi + 1)\cos x + \frac{7}{720} \cos 2x - \frac{\sqrt{2}}{90}\cos 3x - \frac{1}{1050}\cos 6x \bigg).$

I will then likely have to use Maple to determine the third-order term $\delta^{3}z_{3}(x)$ and evaluate $z_{3}(x)$ at $x=0$ and $x-\pi/8$.

My starting point is to use the theory for a regular perturbation (since the modulus of $\epsilon$ is much less than $1$).  For the unperturbed equation, I could set $\epsilon=0$ as that would give a simple differential equation which should be solvable.  I can then see that $1/{1+\epsilon z}$ can be expanded to second-order in $\epsilon$ as $1 - \epsilon z + \epsilon^2 z^2 + O(\epsilon^3), which looks promising.  Could someone advise how I put this together?  Do I then have to multiply the unperturbed solution by the expansion in $\epsilon$? 

Dear Users!

Hoped everyone is good. I am facing problem to write the following sigma notation for any m.

Please help me to fix this problem. Thanks

Probably a stupid question, but is there a free version of Maple for students or is this something we have to purchase?  

 

First 838 839 840 841 842 843 844 Last Page 840 of 2427