Hello.

I know that I can use Maple to fit a curve on  data points (Fit command).
However, I'd like to plot (1sig, 2sig, 3sig) confidence contours for the fit parameters.
Is there any way how to do this in Maple 12?

Thanks.

I just wonder what you guys thing about the direct search optimization package
written by Dr. Sergey Moiseev which is available at the Application center

www.maplesoft.com/applications/view.aspx

Can it be used as a substitute for the global optimization toolbox ? 
 

If I have a sequence of points occurring at specific times.  How can I plot them on a graph so that they appear at a relative time scale? 

Let's take 4 points as an example. 

a:=[ [3,2,11:15] , [4,5,11:17] , [6,3,11:25] , [7,9,11:49] ]

How would I get these points displayed insequence at their relative times?

I have a program (its content is not essential, but I inserted it. The values in outer loops were changed to have only one set of initial data).

printlevel := 5;

for a from .5 by .5 to .5 do

for m from -.999 by .4995 to -.999 do

for n from -5 by 5 to -5 do

ODE := diff(diff(x(t), t), t)+.2*(diff(x(t), t))-x(t) = a*sin(t); k[1] := 0;

ics[1] := x(k[1]) = m, (D(x))(k[1]) = n;

for i to 30 do

X[i] := dsolve({ODE, ics[i]});

u[i] := rhs(X[i]);

Seems that I don't know how to organize a nested loop.

Simple example:

for i by 2 to 19 do

for k from 2 to 26 do

A := i*k

end do

end do

Nothing is returned here (when I expected a list of different values for A).

___________________

Hi everybody,

I want to plot a function f(x,y). I know I can do this with plot3d(f(x,y),x=-1..1,y=-1..1) and this works. But now I want to do something else: x and y fulfill g(x,y)=0 and I want to create a new function h(x) that is basically the old function f(x,y) except that for each x he has to take the y so that g(x,y)=0. Then I can do plot(h(x),x=-1..1). 

OK, this is easy if g(x,y)=0 can be solved for y. But life ain't simple ;-) and in my example this is impossible. I hope anybody has an idea?

Regards,

Casaubon

I have a parameter r which is between a and R, a less then r, r less then R - real or complex numbers, whatever, I need to integrate following function:

1/4*Pi*int(dr1/|r-r1|*int(dr2/|r-r2|*int(dr3/|r-r3|,r3=a..R),r2=a..R),r1=a..R);

The problem is that Maple always says, that it's unable to determine if r is between a and R, how can I assign r into this interval?

[note: Axel Vogt, 20 May: formatted the posting, such that it is readable]

I create equations with so many terms.

And MAPLE print out me NOT all of terms of equation.

for example,              ...........  +`[...10000 terms...]`+..............

But, I want to See these terms!!!!

How can I SEE these terms??

Thank you for reading my Quastion!

Hello, in the attached graph I have scheduled some inequalities, but it shows error, you would know what is that?

Hello, anyone knows why the next page graphics not anything in the server, Maple works,

but seeing it from the link does not display anything,

what should be?

I'm missing a command to generate Bell Polynomials in Maple. Only Bell numbers can be computed. Or are they somehow tightly related so that a separate command would be overkill ?

Hi

I am trying to solve a pde equation with matlab , but i cant, please guide me:

eq_u_0 := diff(u[0](x,t),`$`(x,4))+a*u[0](x,t)+diff(u[0](x,t),`$`(t,2))=0;

bc_u0 := u(0,t) = 0, u(pi,t) = 0, diff(u[0](x,t),`$`(x,2))(0,t) = 0, diff(u[0](x,t),`$`(x,2))(pi,t) = 0, u[0](x,0) = a*sin(x), diff(u[0](x,t),t)(x,0) = 0;

pdsolve(eq_u_0,{bc_u0});

with convert and FormalPowerSeries one can find power series of some functions. WIth the function int one can sum / make integrals of expressions.

I am new to Maple and have tried finding the sum function of the product of two power series :

in both cases sum of n from 0 to infinity
 
how to program Maple to find the sum function of f(x), g(x)
and g(x)*(f(x) ?

a:=[ [ 6,3],[5,6],[6,6],[5,8],[4,7],[5.5,10],[6,12],[5,11],[5,13]]

how do I select the list sets, say that contain in the first position 6.  Something like

if a[i][1]=6 then put the pair in a list.

Hi everyone.

For example, I have

...

plo[i] := plot([u[i](t), v[i](t), t = k[i] .. k[i+1]])

...

Now I want to draw all these plots on the same picture:

with(plots):display([plo[1],plo[2],plo[3],..........................................]);

How can I place all plots in brackets? Number of iterations varies, and there are a lot of sets of iterations, so I need a universal way to do it.

Tried to use Array() but didn't suceed.

Thanks for help!