Can anyone explain the rationale behind plot command, which generates the same graph differently for an identical function?

> restart;
> with(plots); with(plottools);
m := 4;
k := 20;
omega[n] := sqrt(k/m); 
T[n] := evalf(2*Pi/omega[n]); 
td := 0.5*T[n]; 
ratio := evalf(td/T[n]);
p3 := 4*sin(Pi*t/td);
dur := evalf(2*Pi);

enc := u(0) = 0, D(u)(0) = 0; 
eq := m*(diff(u(t), t, t))+k*u(t) = p3; 
sol := unapply(rhs(dsolve([eq, enc])), t); 

I have a function 

f:=x->(x^2 + 2*m*x + m + 1)/(x -1);

Now I want to define a function as simplify(diff(f(x),x)), i did 

g:=x->simplify(diff(f(x),x)). Then, i could not calculate g(3). Please help me. Thank you.

http://aw.twi.tudelft.nl/~koekoek/askey/ch3/par28/par28.html

Which function represent phi in 3.28.1

Hi,

I would like to ask, what is the way to use the Direct Search (or other similar package) to solve not a static optimization problem for, say f(x_1,x_2), but the dynamic one, say for f(x_1(t),x_2(t),t)?

I think there must be the way of adapting the algorithm for that.

Any ideas welcome.

Best

Want create my own geometric distribution with rather only natural values i.e. without zero:

c:=Distribution(PDF=(t->PDF(Geometric(p),t-1)));
But such form is luck to calculate even Mean value:

simplify(Mean(c)) fails. How to apply to get even Mean(F(c)) for complicated enough F?

In which context runs subj? Why all time when i start debugging with 'stoperror' it stops at 'interface' call i.e. in debugger

interface('prettyprint')

doesn't return 3 as worksheet returns. It return interface('prettyprint') i.e. symbol.

I have a matrix with values of interest from .04<4 and I want to round all the terms that are ~10^-10 and lower to zero.  Is there an easy way to do this?  fnormal only works for lists!

The following example shows some typical computations with Householder- or reflextion matrices. Why are the second and third variants only slightly better than the first one? Could we get a real speedup without rtable/NAG/BLAS/etc. tricks?

$ maple15
    |\^/|     Maple 15 (X86 64 LINUX)
._|\|   |/|_. Copyright (c) Maplesoft, a division of Waterloo Maple Inc. 2011
 \  MAPLE  /  All...

I'm writing a maplet and I would like to put in my MathMLViewer two equations. Each equation in different line. How can I do this using MathML[Export]() while setting the value of MathMLViewer ?

As an intermediate result in a longer calculation I got the following equation:

Isolating p, Maple delivers

Calculated manually,...

I'm working on a project that involves solving a large system of polynomial equations using 'solve', with some inequalities of the form "a<>b" or "a<>0". Initially, I was using Maple 9.5 with great success; the machine I was using was upgraded to Maple 11, and now 'solve' struggles with problems that were previously easy (<30 seconds).

Were there any major changes in 'solve' between these two versions? (No, according to the documentation; Yes, according to the source....

So I have something like this:

> with(student);
> with(plots);
> with(plottools);
> with(DEtools);

with the information in the page

http://aw.twi.tudelft.nl/~koekoek/askey/ch3/par28/par28.html

what should do?

i guess to use Hn directly and use sum(Hn*z^m/m!)

is it possible in Maple?

Solve the equation (x-2)^{(25-x^2)}=1