What is the command for filling the space between two circular spacecurves? Something like a Saturnian ring?

Hi there,

I'm trying to get the values from the output of a dsolve command.

I have a system of differential equations:

de1 := diff(V(t), t) = V(t)-(1/3)*V(t)^3-W(t)+Ie;
de2 := diff(W(t), t) = 0.8e-1*(V(t)+.7-.8*W(t))

For a range of the independent variable t and for some given values of the parameter Ie, I would like to get the value of the dependent variable V, as well as its minimum/maximum values for each Ie.

Can anybody suggest a solution please?

This is the worksheet: MaplePrimes_dsolve_min-max.mw

Thanks,

jon

Hello all,

I'm experiencing an aggravating issue when substituting numerical values into a symbolic expression.  I'm using Maple 14.  I believe the issue may be related to floating point precision.  Either way, it's ruining my life.  If anyone has a solution or a work-around, I'd be most grateful.  The following code reproducibly produces the anomolay on my machine:

vx := Vector( 4, symbol=x ):

# A numerator and a denominator. They're equal.
num := (x[1]*x[3]+x[1]*x[4]+x[2]*x[3]+x[2]*x[4])^2:
den := ((x[1]+x[2])^2*(x[3]+x[4])^2):

# Prints true...
evalb(num=den);

# This equals zero
y := 0.1 - 0.1*num/den;

# Prints 0...
simplify(y);

# Prints 0...
simplify(subs(x[3]=1,x[4]=3,0.1-0.1*num/den));

# Prints 0...
simplify(subs(x[3]=1,x[4]=3,y));

# Prints 0.1 x 10^-10 !?!?
simplify(subs(x[3]=1,x[4]=2,y));

Hi,

I need you help to understand this problem.

Let alpha=0.1.;

When I do :  (-0.2)^alpha I get a complex number.

I must find a real number.

What's the problem

Thanks for your idea.

Hello,

I have a little problem. I have a system of ODEs, I can solve (with numeric solve) it but I can't plot it...
Any idea to help ?

Thank you !!

Here is my code:

with(DEtools);
beta := 1; lambda := 5; nu := 1; Delta := 1;


sys := {C(0) = 0, In(0) = .2, diff(C(u), u) = beta*In(u)*s(u)-C(u)/nu, diff(In(u), u) = C(u)/nu-In(u)/lambda+Delta*(diff(diff(In(u), u), u)), diff(s(u), u) = -beta*In(u)*s(u), s(0) = .8, (D(In))(0) = .2}

dsys := dsolve(sys, numeric, [In(u), s(u), C(u)])

plot(dsys[1], u = 0 .. .5, axes = boxed)

And I have this error:

Warning, expecting only range variable u in expression dsys[1] to be plotted but found name dsys[1]


PS: Maple solves the system without any problem:

 dsys(.5);

        [                                              
        [u = 0.5, In(u) = HFloat(0.33305297276372425),
        [                                              

           d                                       
          --- In(u) = HFloat(0.33950703755414946),
           du                                      

          s(u) = HFloat(0.7022184797811994),

                                           ]
          C(u) = HFloat(0.0781621551198558)]
                                           ]

Hi!

I use 'alias' to define variable dependencies.
But when I want to use the result of my calculation as an expression, I bite the dust; I cannot get rid of the dependencies.
Can anyone help me here?

In the simplified example below, I use the independent variables z[3] and z[4] to define z[1](z[3],z[4]) and z[2](z[3],z[4]).
I then introduce a function f[1](z[3],z[4]) and its derivatives df_dz.
I introduce the names dz[1]/dz[3]=Dp[1,3], dz[1]/dz[4]=Dp[1,4] etc...
When I have calculated df_dz[3], however, I cannot insert values of z[2] without affecting Dp[2,3].
Actually, now I have finished my calculation and want df_dz[2,3] to be a regular expression, so that I can insert values of
z without affecting the Dp:s.
How can I convert df_dz[3] to a regular expression?

alias( seq(z[i]=z[i](seq(z[j],j=3..4)),i=1..2)):
alias(f=f(seq(z[j],j=3..4))):
alias(seq(seq(Dp[m,n]=diff(z[m](seq(z[j],j=3..4)),z[n]),m=1..2),n=3..4)):

f:=z[1]-z[2]*z[3]:

for j from 3 to 4 do
df_dz[j]:=diff(f,z[j]);
od:

df_dz[3];
subs(z[2]=4711,df_dz[3]);
whattype(Dp[1,3]);
whattype(z[3]);




function

indexed

Hi there,

I'm relativly new to maple and right now I'm trying to plot a couple of points (as part of a TSP-Solution) by doing the following:

with(plots):

p := [[565,575],[685,595],[700,580],[770,610]];

pointplot(p, connect=true);

Unfortunately it doesn't number the points. What I'm trying to do is that the above four points have a number or an index so I can see which point is where. 

Is there a way to do this?

thanks

restart:

Eq1:=x=(phi*(theta[m]-1/x)/theta[m]+(1-phi)/2);

solve({Eq1},x);

x := (1/4)*(phi*theta[m]+theta[m]+sqrt(phi^2*theta[m]^2+2*phi*theta[m]^2-16*phi*theta[m]+theta[m]^2))/theta[m];

Now plotting the first root

p1:=plot(subs(phi=0,x),theta[m]=0..14,color=red):
p2:=plot(subs(phi=0.1,x),theta[m]=0..14,color=green,linestyle=2):
p3:=plot(subs(phi=0.5,x),theta[m]=0..14,color=blue,linestyle=3):
p4:=plot(subs(phi=1,x),theta[m]=0..14,color=black):
display({p1,p2,p3,p4},axes=boxed,view=[0..14,0..1]);

Note: I'm unable to reproduce the following plot. Any idea?

 SA.mw

The following is taken from p 328 of the Maple TA manual.

\documentclass[12pt]{article}
\usepackage{ed}
\begin{document}

\begin{question}{Formula}
% This question uses randomized variables.
\qutext{Find the formula for the inverse of the function
$f(x)=\var{a}x+\var{b}$.}
\answer{(x-\var{b})/\var{a}}
\code{
$a=int(rint(9)+2);
$b=int(rint(9)+2);
}
% The above code chooses random integer coefficients a and b between
% 2 and 10.
\end{question}

\end{document}

Running it though the online compiler generates the following when Service Pack 1 checkbox is UNCHECKED

Note: Successful conversion from LaTeX to QU does not guarantee that your testbank will load into Maple T.A. correctly. Please examine the contents of the box below for any additional errors that may be generated when your testbank is loaded into Maple T.A.

l2e_pkIwct.qu: Error at line 3.
Line ended while reading variable name (missing "=" at end of name?)

(A response of "[OK]" indicates that there should be no problems uploading your converted testbank file to Maple T.A.)

and the following when it is CHECKED

...

[30] [31] [32] [33] [34] [35] [36] [37] [38] [39] [40] [41] [42] [43] [44]
[45] [46] [47] [48] [49] [50] [51] [52])
(/usr/local/latex2edu/newt4ht/texmf/tex/generic/tex4ht/html-mml.4ht)
(/usr/local/latex2edu/texmf/tex/local/ed/ed.4ht [53] [54]))
(/usr/local/latex2edu/newt4ht/texmf/tex/generic/tex4ht/html4.4ht)
(/usr/local/latex2edu/newt4ht/texmf/tex/generic/tex4ht/unicode.4ht)
(/usr/local/latex2edu/newt4ht/texmf/tex/generic/tex4ht/html4-uni.4ht)
(/usr/local/latex2edu/newt4ht/texmf/tex/generic/tex4ht/mathml.4ht)
(/usr/local/latex2edu/newt4ht/texmf/tex/generic/tex4ht/html-mml.4ht)
(/usr/local/latex2edu/texmf/tex/local/ed/ed.4ht))
No file l2e_yQsgNQ.aux.
)
*
! Emergency stop.
<*> \def\production{}\input l2e_yQsgNQ.tex

Output written on l2e_yQsgNQ.dvi (54 pages, 25316 bytes).
Transcript written on l2e_yQsgNQ.log.

What's going on?

I've got two matrices A and B that I've verified are similar using IsSimilar()

The description of the function says the third argument 'P' stores the transformation matrix. Does this P have to be initialized and if so, then how can it be done?

I'm a student. I want to purchase oline Maple Student, but I can't.

Please send me a feedback as soon as possible. Thanks.

I am attempting to plot an initial value problem in Maple 18.  I have my equation defined, as well as a general solution and two particular solutions at y(0)=3/4 and y(0)=1/2.  To graph, I entered the command

DEplot(de,y(x),x=-3..3,{[0,1/2],[0,3/4]},dirgrid=[12,12],color=black,linecolor=blue,thickness=2);

but instead of returning a graph, the software gave me the error message

Error, (in DEtools/DEplot/CheckDE) extra unknowns found: sinx

The Maple support site lists this as an unknown error, and as a new user, I'm not sure what to do.  What does this mean?

hello

I have a system of equation. to solve i using DirectSearch pakage. i think in answer DirectSearch use rounding number but in answer it's second and third digit after the decimal point is important. please help me. bbw.mw

if b and c increase 0.01 and -0.01 it's ok. another main question is why by increasing the intensity answer don't change.

Hi everyone,

I've been having a problem trying to correct this procedure so that it will accept empty tables as well. This is the code I've been messing up with:

6.10.mw

If someone could help me with this, it would be really appreciated.

Hi,

So I have a simple while loop and keep receiving this error and have no idea how to fix it. I don't get why it works for the first few, and then the error starts. Here is my code and output:

Input:

restart;
n:=0:
x:=Pi/3;
m:=0:
total:=0:
while abs(1/2-total) > 10^(-8) do
k:=((-1)^(n/(2))*x^(n))/(n!):
total:= total +k;
evalf(total, 11);
print(%);
n:=+2;
end do:

OUTPUT:
1
- Pi
3
1.
0.45168864437
Error, cannot determine if this expression is true or false: 0 < -50000001/100000000+(1/18)*Pi^2

Why does it think i'm asking if something is greater than 0? Any help??