I want to plot curves with Japanese title and labels. I want to plot through an external java program. Is it possible?

When I tried, I am getting control character unexpected error.

Hi, Everyone.

I'm trying to figure out how to make a graph (actually a "line") change color based on certain criteria so I can create a Maple Application.

I've figured out the following, which indicates the sign of the graph by the color, but it's way more complicated than it ought to be.

It's also not robust as I need it to be, this happens to depend on the zeros occuring at integer values.  I could deal with that on a case by case basis, but not for something automatically.

I'm at a loss; maybe I'm not Googling the right things... maybe Maple can't do this (although I'd heard color tools were being added in Maple 16).

Does anyone have any suggestions?

Thanks In Advance!

Joseph

Am I right that using methods of the 1-forms (that implemented in liesymm or DESOLV), we can always generate determining equations for ODE, that solved for highest derrivative?

Hello, could someone please help me to understand this warning that i'm getting?

Warning, `N4ds` is implicitly declared local to procedure `SFunc`

Thanks

Maple is providing answers with the gamma function and LambertW function.  I am not a math major.  Although I did take advanced engineering mathematics in grad school, I am not familiar with transcendental functions.  

I did look at some of the help.  I tried to use the convert command to convert to elementary functions, however it did not help.  I looked at the definition for the Gamma function, GAMMA(z), and the upper incomplete Gamma function, GAMMA(a,z) that uses the generalized hypergeometric function 1F1.  Looking at the generalized hypergeometric function 1F1 defintion, hypergeom(n, d, z), it references the general pochhammer function, pochhammer(z, a), and product functions, product(f,k=m..n).  This all seems way too complicated.  

eqn1:= R = exp(t/MTBF) * (t/MTBF) * sum(1/(i!), i=0..n);

eqn2:= MTBF = solve(eqn1,MTBF);

Is there any easy way to convert to use basic functions like sum, int, etc.??

I am wanting to use the answer in Excel.  I have found some custom Excel VBA code that I could use to implement these non-standard function in Excel, but that seems like a lot of effort.

Lambert W function for Excel, work on real and complex number (VBForums)

Excel Add-Ins Incomplete Gamma Function (PC Review)

Update: Some more background on the problem is as follows. R is the expected reliability say 90%. MTBF is mean time between failures in hours. t is the time period of operation excluding allowable downtime. n is the allowed number of failures. I am looking to find the MTBF needed to achieve a certain required reliability, with the number of failures allowed, and the up-time operational period, being defined by the requirements as well.

To clarify the simplification does not necessarily have to be closed form, and an approximation would be acceptable. (in Excel I currently am using iteration over a number of rows, then using the Goal Seek tool to find the MTBF, however this is less than desirable to have to do everytime when looking at different trade studies) I am not that familiar with how to use assumptions in Maple, and do not know what would be the appropriate assumptions to make when simplifying. Thanks for your help.

Hi, I am learning GR and new to maple. I played with maple's DifferentialGeometry package.

While duplicated the maple GR example, in DGsetup([t, r, theta, phi], M1, verbose),

I got the error,

DGsetup([t, r, theta, phi], M1, verbose)

The following coordinates have been protected:
[t, r, theta, phi]

The following vector fields have been defined and protected:

Error, (in Typeset:-Tdisplay[true]) `alias` does not evaluate to a module

 I don't know how to fix it. Can anyone kindly help?

Thanks

I'm trying to build a Maple procedure that will generate vector fields on a metric with certain properties. Working with metric g over the coordinates {u,v,w}, call the field X = (a(u,v,w), b(u,v,w), c(u,v,w)). The field should satisfy <X, X> = 0 and have the directional covariant derivative of X in the direction of each coordinate vector field = 0 (with resepct to the Levi-Civita conenction).

Basically, these conditions yield a system of 3 PDEs and an algebraic expressionin terms of a,b,c. I've been trying to solve them using pdsolve, but I'm getting the error message:

>Error, (in pdsolve/sys) the input system cannot contain equations in the arbitrary parameters alone; found equation depending only on _F1(u,v,w): _F1(u,v,w)

I've attached my worksheet. Can anyone help me out?

Thanks! ppwaves.mw

Hello,

i am trying to get a latex output which equals the displayed maple equations, where i am using alias/surpress or declare to shorten the dependencies of my derivatives. Whereas in Maple this looks how i want it to look like, i can not get the latex() command to apply the aliases. Instead it replaces the short forms with the long terms before creating the latex code. Is there anyway to get a latex output directly from the displayed math WITH alias?

2 small examples:

Both yield  $ X \left( a,b,c \right)$ whereas i would like to get X. Simply substituing changes the partial differential symbol to "d" in the latex output and makes the equations unuseable, so this is no option i guess..

Thanks a lot for any suggestions!

I have the following integral equation to solve numerically:

v(x,t)=1 - h*\int_0^t JacobiTheta0(1/2x , \pi i s) v^4(1,t-s)ds

where h is a numerical parameter, and v(1,t) = 1-h*\int_0^t \theta_3(r)v^4(1,t-r)dr (theta3 is Jacobi theta3 function).

So I want to use an iteration method that will converge numerically to the solution, where v(1,0)=1.

How to use maple for this?

I want also to find the rate of convergence to the numerical solution.

 edit: I should note that v(x,0)=1, even though it's implied from v(x,t) above.

I have three cylinders as follows.

with(plots);
with(plottools);
p1 := cylinder([1, 1, 1], 1, 1);
p2 := cylinder([1, 1, 1], 2, 3);
p3 := cylinder([1, 1, 1], 1, 4);

I want to display in order where p1 is at the bottom, then p2, and then p3 on the top of p2. Thank you for your assistance.

I have the following PDE system to solve numerically and I am not sure how to use maple to solve it.

v_t = v_{xx} for 0<x<1 , t>0

v(x,0)=1

v_x(1,t)=-hv^4(1,t) (where h is some numerical number);

v_x(0,t)=0

To solve this pde numerically I need to use the following condition on v(1,t):

v(1,t) = 1-h*\int_{0}^t \theta_3(\tau)v(1,t-\tau)^4d\tau

this is the numerical boundary condition, where \theta_3 is Jacobi theta3 function.

I don't see how can I use maple for this numerical pde problem.

Here's my attempt at solution:

[code]

 PDE := diff(v(x, t), t) = diff(v(x, t), x, x);

JACOBIINTEGRAL := int(JacobiTheta3(0, exp(-Pi^2*s))*v(1, t-s)^4, s = 0 .. t);

IBC := {`&PartialD;`(v(0, t))/`&PartialD;`(x) = 0, `&PartialD;`(v(1, t))/`&PartialD;`(x) = -0.65e-4*v(1, t)^4, v(x, 0) = 1};

pds := pdsolve(PDE, IBC, numeric, time = t, range = 0 .. 1, spacestep = 0.1e-2, timestep = 0.1e-2, numericalbcs = {v(1, t) = 1-0.65e-4*JACOBIINTEGRAL}, method = ForwardTimeCenteredSpace)

[/code]

But I get the next error message:

Error, (in pdsolve/numeric/process_IBCs) improper op or subscript selector

How to fix this or suggest me a better way to solve this pde numerically?

multiple_parameters_-_sipa.mw 

@Carl Love please help me .. I tried solve a ode system of equation problem numerically with your code, using maple 15 but i get errors.

I dont know how to plot multiple solutions on one graph.

I want to vary one of the parameters....

e.g s1=0, s1=0.5, s1=1 where other parameters are kept constant

> restart;
> g := 1/7; alpha1 := .75; alpha2 = .5; `&alpha;v` := .3; xi := .4; `&mu;v` := 1/30; mu1 := 0.434931506849e-4; mu2 := 0.410958904109e-4; `&mu;h` := 1/(60*365); Nh := 10000; a := 1; Nv := 1000; s2 := .5; br := 4;
> A1 := diff(S0(t), t) = (1-s1)*mu1*Nh+(1-s2)*mu2*a-br*alpha1*S0(t)*V1(t)/Nh-`&mu;h`*S0(t);
> A2 := diff(Z0(t), t) = xi*Z(t)-br*alpha2*Z0(t)*V1(t)/Nh-`&mu;h`*Z0(t);
> A3 := diff(Iota(t), t) = br*alpha1*S0(t)*V1(t)/Nh-g*Iota(t)-`&mu;h`*Iota(t);
> A4 := diff(Y(t), t) = br*alpha2*Z0(t)*V1(t)/Nh-g*Y(t)-`&mu;h`*Y(t);
> A5 := diff(Z(t), t) = s1*mu1*Nh+s2*mu2*a+g*(Iota(t)+Y(t))-xi*Z(t)-`&mu;h`*Z(t);
> A6 := diff(V0(t), t) = `&mu;v`*Nv-br*`&alpha;v`*V0(t)*Iota(t)/Nh-`&mu;v`*V0(t);
> A7 := diff(V1(t), t) = br*`&alpha;v`*V0(t)*Iota(t)/Nh-`&mu;v`*V1(t);

> PD := unapply({A || (1 .. 7)}, s1);
> IBC := {S0(0) = 450, V0(0) = 100, V1(0) = 70, Y(0) = 50, Z(0) = 250, Z0(0) = 50, Iota(0) = 100};
> sList := [0, .6, 1]; Colours := table(`~`[`=`](sList, ["blue", "red", "green"])); for s1 in sList do pds := dsolve(PD(s1), IBC, numeric); Plots[s1] := pds:-plot[display](Iota(t), t = .5, linestyle = "solid", colour = Colours[s1], legend = sprintf("s1=%2.2f", s1), title = "Velocity Profile") end do;
Error, (in dsolve/numeric/process_input) invalid argument: {S0(0) = 450, V0(0) = 100, V1(0) = 70, Y(0) = 50, Z(0) = 250, Z0(0) = 50, Iota(0) = 100}

> plots:-display([seq(Plots[s1], `in`(s1, sList))]);
Error, (in plots:-display) expecting plot structures but received: [Plots[0], Plots[.6], Plots[1]]

Download multiple_parameters_-_sipa.mw

Hi EveryOne!

Let matrix M over GF(2ⁿ)/f(x). How to compute:

1. Rank of M.

2. Sum of the elements and Square of this Sum in any row of M (i.e: Sum=m_i0 + m_i1 + ... + m_in, and Sum² ).

3. New matrix M, where in row i  m_i0 = m_i0/Sum², m_i1 = m_i1/Sum², ..., m_in = m_in/Sum².

Please help me. Thank you very much.

Given a polynomial in several variables is it possible to split it so that all the coefficients of the monomials are +1 or -1.

Example:

p:=-z+2*x+4*y-3*x*y.

I would lie to obtain

f:=-z +x+x +y+y+y+y -x*y-x*y-x*y.

Given a polynomial expression I would like to obtain a list whose entries are the positive entries of the polynomial with multiplicity given by the coefficients:

Example:

Given the polynomial expression: p:=x^2*y-2*y*z+3*x^2+2*y-z

The positive terms are: x^2*y, 3*x^2, 2*y

Thus I would like to obtain the list L:=[x^2*y , x^2 , x^2 , x^2, y , y].

Notice x^2*y appears once since the coefficient is 1.

x^2 appears three times since the coefficient is 3.

y appears two times since the coefficient is 2.