PatrickT

Dr. Patrick T

2168 Reputation

18 Badges

17 years, 174 days

MaplePrimes Activity


These are answers submitted by PatrickT

Do you mean this?

tvalues:= fsolve(subs(curveq, surfeq), t=-2..2);

which will return the same answers as before since they were in the interval specified. However, say:

tvalues:= fsolve(subs(curveq, surfeq), t=0..2);

                        tvalues := 1.291980321

 

plot understands the discontinuity without the option discont=true:

plot(tan(x), x = -5..5, y = -5..5, style=point, grid=[50,50]);


plot(tan(x), x = -5..5, y = -5..5, grid=[50,50]);

Like Doug I too have had difficulty uploading an animated gif on mapleprimes. I'd like to give it another try. As the gif is too large to be uploaded with mapleprimes, I uploaded it with imageshack:

http://img534.yfrog.com/img534/2471/phaseportrait3.gif

I saved the gif using the context-menu option "Export As" and select "Graphics Interchange Format (GIF)".

Let's see if it will display here and show the animation:

 

The animation was created with the code above, using
 

DEtools[DEplot](
  {ode1,ode2},
  {x(t),y(t)},
  t=0..20,
  [ic1,ic2,ic3,ic4,ic5,ic6],
  scene=[x(t),y(t)],
  x=-1.5..1.5, y=-1.5..1.5,
  arrows=smalltwo,
  dirfield=400,
  size=magnitude,
  color=magnitude, 
  linecolor=t,
  thickness=3,
  animatecurves=true,
  animatefield=false,
  numframes=20
);

you may want to start with using "shorthand" notation, like defining:

a := ((R4*l^2+sigma)/R4)^(1/2);
b := 7/6*varpi*(R1-1)*l;

etc., and using subs to substitute the as and bs into the messy expression.

Another thing you may want to do is replace the cosh and sinh by their exponential forms and see if that simplifies to something more compact.

your feedback is urgently expected.

not sure I understand your question ... do the following help?



with(plots):
plottools[reflect](p, [0,0],[1,1]);
plottools[rotate](p, Pi/2, [0,0]);

?pointplot will point you to the help page, does it do what you're looking for?

you will need to define Vectors (not vectors) to use Robert's code, thus:

X:= Vector([-1,0,1]);
Y:= Vector([1+I*2,3+I,4+3*I], datatype=complex);

I don't think you can control the position of the labels automatically. But you can remove them altogether and then redefine text labels that you could place where you wanted to. To remove the labels, you can do something like this:

plots:-implicitplot3d(....., labels=[``,``],.....)

or

plots:-display(theplot, labels=[``,``],.....)

From the context menu, you may be able to do something (untested):

http://www.maplesoft.com/support/help/AddOns/view.aspx?path=worksheet/plotinterface/axesstyle

VectorCalculus:-Jacobian(subs(x(t)=x,y(t)=y,z(t)=z,map(rhs,[rossler_sys])),[x,y,z]);

oh, but then it evaluates the Jacobian at (x,y,z) rather than (x(t),y(t),z(t))...

It all depend on how good you want the figures to look. Currently 2D plots look a lot better than 3D plots. There are also differences between standard and classic gui. I personally save a basic plot in postscript format and add fancy decorations with pstricks and pstricks-add.

Of possible interest:

http://www.mapleprimes.com/forum/changethicknessaxesandtickmarkssimpleplot

http://www.mapleprimes.com/forum/tickmarkshowsetspacingandthickness

http://www.mapleprimes.com/forum/exportingimplicitplot3dstructuresstandaloneuse

http://www.mapleprimes.com/forum/poorqualitygraphicsexportedeps

to get you started,

SteL:=-1: SteV:=-1: h:=1: R1:=1: #select parameters that make sense
is(SteL<0 and SteV<0 and h>0 and R1>0); # this should return true if your assumptions are satisfied
eq := sqrt(Pi)*beta=-1/h*SteL/(exp(h^2*beta^2)*erf(h*beta))
+R1*SteV/(exp((beta/R1)^2)*erfc(beta/R1)); # this is your equation
fsolve(eq,beta); # this is for numerical solutions, different  from solve(eq,beta)
is(%>0); # this should return true if the solution has beta>0
restart;

eq := theta0 = 10 * (2.718281828 - exp(S) ) / exp(S);
p := plots:-implicitplot(eq, theta0=-1..1, S=-1..1):
plots:-display(p);

these {} define a set (no order), while these [] define a list (with order), so what method have you tried?

the square brackets define a "list", replace them by round brackets.

[] -> ().

if problems arise, copy-paste your lines of code together with the error message (if any).

more code would make it easier to understand your needs, perhaps the following is of relevance:

http://www.mapleprimes.com/blog/scumath/howplotgeneralspacecurvemaple

http://wamp.mapleprimes.com/forum/intersectioncurve

First 9 10 11 12 13 14 15 Last Page 11 of 24