This is a problem of Russian in 2005: Solve inequality 

(3-x-sqrt(5-x^2))/(cos((2*x-7)/4)-cos((x-5)/4))>=0.

The problem can solve by hand. We only solve

(3-x-sqrt(5-x^2))/((2*x-7)/4-(x-5)/4) >=0.

I tried

eq1:=(3-x-sqrt(5-x^2))/((2*x-7)/4-(x-5)/4);

solve(eq1>=0,x);

I got

RealRange(1,Open(2)), RealRange(Open(2),5^(1/2)).

If I solve directly by

eq:=(3-x-sqrt(5-x^2))/(cos((2*x-7)/4)-cos((x-5)/4));

solve(eq>=0,x);

I don't get the exact answer. 

P.S. with Mathematica, I only use

Reduce[(3 - x - Sqrt[5 - x^2])/(Cos[(2 x - 7)/4] - Cos[(x - 5)/4]) >=
0, x, Reals]

I have an expression 

f:= x1 * x2 * x3^4;

I would like to do something like subs( (x1,x3) = (x3,x1), f) and get 

x1^4 * x2 * x3

back.

However this does not work. Is there some _lightweight_ way of doing it?

Hello everyone,

Here is my question,

restart:with(plots):

Eq:= diff(f(y),y$3)+f(y)*diff(f(y),y$2)-diff(f(y),y$1)^2-S*(diff(f(y),y$1)+y/2*diff(f(y),y$2))=0;

bc:=D(f)(0)=1,f(0)=0,(D@@2)(f)(b)=0, f(b)=S*b/2;

How we can plot b vs S?

Thanks

Hello,

I am solving eq1 and then plot it using odeplot. Then I am extracting data from odeplot using op[1,1]. From that I only care about  one point so I am using pt[1] := dt[1]([1, 2]); then I ploted it using pointplot. I run a do loop and at the end I am ploting all the points.

My code works fine but I am trying to have a line contacting each points. I tried  style=line but it did not work.

How can I have a line contacting each points on my final plot (which is display([pl[1], seq(allpl[k], k = 1 .. 3)])) ?

This is my code:

> restart; with(DEtools); with(plots);
> with(DEtools); with(plots);
> A := 0.2e-1; B := 10^(-5); k := 0;
> eq1 := diff(X(t), t) = -(A+B*X(t))*X(t);

> ic[1] := X(365*k) = 1000;
> s[1] := dsolve({eq1, ic[1]}, X(t), range = 0 .. .365, numeric);
> p[1] := odeplot(s[1], [[t, X(t)]], t = 0 .. .365);
> dt[1] := op([1, 1], p[1]);
> pt[1] := dt[1]([1, 2]);
> pl[1] := pointplot(pt[1], axes = boxed);
> for k to 3 do
tk := 365*k;
A := rhs(s[k](tk)[2]);
ic[k+1] := X(tk) = 500.*A;
s[k+1] := dsolve({eq1, ic[k+1]}, X(t), range = tk .. 2*tk, numeric);
p[k+1] := odeplot(s[k+1], [[t, X(t)]], t = tk .. 2*tk);
dt[k+1] := op([1, 1], p[k+1]);
pt[k+1] := dt[k+1]([1, 2]);
pl[k+1] := pointplot(pt[k+1], axes = boxed);
allpl[k] := display([pl[k+1]])
end do;
> display([pl[1], seq(allpl[k], k = 1 .. 3)]);

Thank you

I installed the shared library from http://www.maplesoft.com/applications/view.aspx?SID=1691

and it seems to have loaded ok. Now I can do

with(share);
See ?share and ?share,contents for information about the share library
                               []

When browsing help, I noticed it has a function called IntSolve() for solving linear integral equations. Maple itself has function called intsolve() which is build-in.

Since Maple is case sensitive, I assumed when I type IntSolve() then it will find the shared version of IntSolve() automatically and not call Maple's own intsolve. And that is indeed what happed. I typed an example from the help on share as you can see

eq1 := p(x)-1/2*Int(x*y*p(y),y=0....1)=5/6*x;
IntSolve(eq1,p(x));
                             7775  
                             ---- x
                             7776  
intsolve(eq1,p(x));
                            p(x) = x

Ok, So one of them is a little wrong :).

But suppose I want to explicitly add the share library there, in the call, as prefix, so that when I look at the code later, I can see which function I am calling. How should one do this? I tried

share:-IntSolve(eq1,p(x));
    Error, `share` does not evaluate to a module
share[IntSolve](eq1,p(x));

The strange thing, is that this worked

IntSolve[IntSolve](eq1,p(x));

Basically, all what I want is to way to add the share name to all the calls somewhere, so I know, when looking at the code, which call is from which library or package.

Any idea how to do this? So to avoid any possible name clash since the share library seems large.

I am using Maple 17.02

I am using 17.02. I looked at this question, http://www.mapleprimes.com/questions/35783-Inverse-Cosine-And-Sine-Transform but I do not understand the answer.

I have a function F(w), and want to find f(t) using inverse cosine forurier transform.

Here is F(w)

F:= w-> cos(w*L)*cos(w*x)/(w^6+w^2*k);

All parameters are positive.

How can I use the answer shown above, or another function to find f(t)?

I can offcourse simply apply the definition myself, and I did:

int(F(w)*cos(w*x),w=0..infinity) assuming positive;

but the result is

So it does not converge. And I just wanted to check using a build-in function to verify. Strange than Maple does not have these in its inttrans package.

with(inttrans);

I am trying to use Maple to produce an image that will be part of a math paper.  The final image must be in .eps format.  Whenever I produce a plot in Maple17 and click on it, it looks like the image has the shape of a square (regardless of what's inside it), perhaps with a lot of empty space.  I want an image without a lot of empty space, and I'd like Maple to do this.  If I crop a square image using some other application, it is not going to improve the quality of the image and might hurt it.  

If the answer to question in the Title is "no", can anyone suggest a good way to produce a good .eps image from a Maple plot that is cropped properly to eliminate as much empty space as possible?  I have Photoshop, but I don't know if this is the best way to go.

I wrote a toy file that plots a simple curve.  I'll upload it and provide a link to it.  If you execute the code and click on the image, you'll see that it is a square with a lot of empty space.

Nov122013.mw

Download Nov122013.mw

Can any workarounds in Maple for debugging the code without procedure? 

With Maple, I tried

solve(sin(Pi*(x+1)/(4*x^2-4*x+2))=cos(Pi*(x-2)/(4*x^2-4*x+2)),x);

And I got

Warning, solutions may have been lost

RootOf(-sin(1/2*Pi*(_Z+1)/(2*_Z^2-2*_Z+1))+cos(1/2*Pi*(_Z-2)/(2*_Z^2-2*_Z+1)))

With Mathematica

Reduce[Sin[(x + 1) \[Pi]/(4 x^2 - 4 x + 2)] == Cos[(x - 2) \[Pi]/(4 x^2 - 4 x + 2)], x, Reals]

I got all solutions.

Hey everybody I am trying to get maple to do some partial differention for me. 

Where 

I have PDEtools on but, it doesnt return an expression correctly. It's I think the problem is because Ψ is a function of (x,y) and f is a function of η and η is a function of x and y . U and ν are constants. v_x and v_y are the x and y velocites. 

What I would like the program to return when I write 

is something along the lines of 

could someone give me a hand on this. I would really apperciate it. 

Best Regards, 

Kyle 

Hi Guys,

I am trying to do the following and here is the issue:

> System*Transfer*Functions;
Power Stage Transfer Function:
> He := s^2/`ωn`^2+s/(`ωn`*Qn)+1;
>
Compensation Netwrok Transfer Function:
> Av := gm*(1+s/`ωcz1`)(1+s/`ωcz2`)/((Cc+Rc)*s*(1+s/`ωcp`));

>
Loop Transfer Function:
> Loop := Vref*RL*(1+s/`ωesr`)*Av/(Vout*Ri*(1+s/`ωp1`)*He);

> sys3 := TransferFunction(Loop);

Warning, transfer-function(s) are not rational-polynomial(s) in s

I want to do the TF and Bode Plot for this complex system... 

I am having trouble getting a certain Maple worksheet to execute.  Using both Maple17 (I think it is 17.02) installed on two separate PC's at home (XP), and using Maple15 via Citrix from a server maintained by my employer, I get a long error message.  It starts:    "Kernel connection has been lost.  You should save this worksheet and restart Maple.  Executing commands in Maple requires a connection to the Maple kernel.  Firewalls have been known to cause problems with kernel connections in Maple.".  If anyone wants, I can add the entire error message to this question. I don't think I actually have a firewall at home, and this error message doesn't seem to explain why my local copy of Maple is having trouble.

I think I updated my home copy of Maple to Maple 17.02 after I started writing the worksheet and before this problem appeared, so there may be an issue with Maple 17.02.  But this would not seem to explain the trouble I have executing the worksheet using Maple15 via Citrix.

I have uploaded the file.  I removed the output to make it a reasonable length.  I created a "toy" file that plots x^2 from x = 1 to x = 4, and does nothing else, and it works just fine.  I apologize for the length of the file I uploaded, but I don't know exactly what is causing  the error.

NewPicture.mw

GS

I copied and pasted this from the example page and I have to get this same graph for an assignment

im using maple 17

when i hit enter the return was

                       (1)

Hi,

I'm having some difficulty in figureing our how to store the results of an fsolve computation in an indexed variable that I can access subsequently to perform plots and do other calculations.  I have a do loop that, for every increment, has fsolve calculate five variable  (x,y,z.. )from a set of equations that have been entered.  It does the calculations just fine.  What I can't seem to do is set up an array variable (or other indexed variable like a vector) so that with every increment of the loop the values that fsolve generates for x, y, z  ect. can be stored.  That just doesn't seem to work, what am I getting wrong? Any help would be most appreciated.

Al Stiegman

Hi,

I need help formulating a biased coin toss.

I want to toss a coin 0 being heads and 1 being tails with probabilty 0.495 of getting heads, how do I show whether heads or tails wins??

Thanks