Hi all!

I would like to isolate a variable from an equation but after trying both solve and isolate I noticed they result in two different solutions. Which one is correct or in other words are there any issues with my syntax resulting in these differences?

The Maple ?Compiler can compile a limited subset of Maple commands to native code. The result is substantially faster than running interpreted Maple code. This article shows how you can save a compiled procedure for reuse.

When a Maple procedure is compiled, a shared object library (dynamic-linked library for Windows) is created on the file system, as a temporary file.  When the Maple...

To my mind I did all the necessary assumptions, to avoid that the summand will become singular when calculating the expected value.

What did I do wrong?

Here´s the code:

with(Statistics):
interface(showassumed=0):
assume(t::posint):
assume(x::posint):
assume(0 < q, q < 1):
X1 := RandomVariable(Geometric(q)):
DARA := t->log(t):
ExpectedValue(DARA(X1));

Hi everyone;

I would like to solve these ODE's.

1) d/dX(X*dΘ/dX)+ρ^2*X*Θ=0;Θ=Θ(X),ρ=constant

2)d/dr(r*d/dr(φ)-m^2(r*φ-n/m^2)=0;φ=φ(r),m,n=constant

Boundry Conditions are:

1)d/dX(Θ(L))=0,  2)Θ(0)=φ(R), 3)d/dX (Θ(0))+d/dr(φ(R))=0, 4)d/dr(φ(0))=0,   R,L=constant

Two of boundry conditions are Joint

Thanks

I'm using Maple 14 on 4-way and 6-way AMD 64 hardware running Windows Vista to numerically estimate inverses of a cumulative distribution function at hundreds of points, and am using the Task model to assign each point's inverse computation to a separate task. The cdf is a weighted sum of thousands of simpler functions, and so each task needs the weights. Initially I had set these up as a global variable, but after reading Darin Ohashi's MaplePrime's post ...

Hi all!

solve({q[R1] >= 0, q[S] >= 0, `d&Pi;`[R1] = 0}, q[R1])

This command gives me a piecewise solution, showing somewhat like:

[] q[S] < 0.13...

[{q[R1] = ...(q[S])   }] q[S] => 0.13...

Is there some way to define that part of this piecewise function for which q[S]=>0.13 as a seperate function? I need to do some more work on q[R1...

Hello there,

so i have a maple algorithms which produces about a thousand set of results. The algorithms produces about 7 values.. and this is iterated about 100 times.

Is there any way I could transfer all these values into a spreadsheet?

basically I want to end up with a speadsheet with the values from each iteration forming a row.

Thanks in advance!

How to flatten a list without ListTools? Thanks.

The Linear Algebra package seems only able to handle Vector and Matrix algebra when the dimensions are given specifically.

It would be very helpful if the package were to be extended to cope with the situation where some or all of the dimentions were given symbollically. Of course, assume could be used to indicate when two dimensions were equal.

Are there any plans in this direction for future releases of Maple?

In the following code Maple simplifies the last expression different, based on whether assume(mu::real, sigma>0) is set or not.

The point is, that the assumptions got -to my mind- nothing to do with the displayed results.

Why does Maple do such weird things?

Here is the code. The last output is different based on whether the '#' before assume(mu::real,sigma>0): is set or not:

restart;
with(Statistics):
Y9 := RandomVariable(LogNormal(b, c)):

I discovered, that Maple did substitute every single parameter but c!

Anyway, when I isolate the problem and do not assume anything for mu and sigma, the problem does not occur.

Therefore I print here the complete relevant code.

Is there a way to "convince" Maple to do those trivial substitutions of c also, without seperating each problematical expression?

Here is the code:

restart;
with(Statistics):

#Definitions and assumptions

My problem is the following:

deq1:=2*[(y(x)*(diff(y(x), x, x))-1)/sqrt(1+(diff(y(x), x))^2)+(diff(y(x), x, x))/(1+(diff(y(x), x))^2)] = cos(x);
I do this way (error anyway):
deq := (y(x)*(diff(y(x), x, x))-1)*(1+(diff(y(x), x))^2)+(diff(y(x), x, x))*sqrt(1+(diff(y(x), x))^2)-cos(x)*sqrt(1+(diff(y(x), x))^2)*(1+(diff(y(x), x))^2) = 0
sol := dsolve([deq, y(0) = 0, (D(y))(0) = 0], y(x), numeric, range = 0 .. 3)

So, my question is, how can I determine the integral of the left side of deq1? ...

Hello,

I have this differential equation:

deq := (diff(y(x), x, x))*y(x)-1+x*y(x)^2 = cos(x);

sol := dsolve({deq, y(0) = 1, (D(y))(0) = 1}, numeric);

After that I have a numerical solution for y, and for diff(y(x),x). But I would like to plot diff(y(x),x,x). (and then diff(y(x),x,x)+diff(y(x),x)*cos(x)). How can I do that?

Any help will be appreciated,

L

Hello there,  How can I compile this? I used   boundarycompiled := Compiler:-Compile(findboundary)    but this didnt work. Thanks in advance

> findboundary:=proc(a,b,c,d,e,f,g,h)  ;
> global A,B,lowerbound,upperbound;
> A, B := solve(a = b, x), solve(c = b, x);
> x1a, x1b := evalf(solve((x-e)^2+(a-f)^2-g^2)); if x1a < x1b then minx1, maxx1 := x1a, x1b else minx1, maxx1 := x1b, x1a end if;
> x2a, x2b := evalf(solve((x-e...

Hi. I have to solve the following differential equation numerically

f := [(diff(y(x), x))*(diff(y(x), x, x))/sqrt(y(x)^2+diff(y(x), x))+diff(y(x), x, x)] = x^6

dsolve({f, y(0) = 2, (D(y))(0) = 1.5, (D(D(y)))(0) = 1}, y(x), numeric, range = 0 .. 10)

Error, (in DEtools/convertsys) unable to convert to an explicit first-order system

I don't know what to do, I've read help and work on this for 3 days now.