Hello, I am currently using Nelder Mead for an optimization problem. I used the following code taken out of the Maple Library (www.maplesoft.com/applications/app_center_view.aspx?CID=1&SCID=18&AID=1198) Concerning this code I have 2 questions: 1. The optimized values I get are sometimes negative. My solution should be limited on positive optimization values. I am not sure how to implement this. 2. The optimization is running quite slowly. Any suggestions how to improve the running time? Here are the most important parts of the code:

Hi to everybody, I'm computing these polinomials:

restart;
 pol:=(-1)^n/(2^n*n!)*(1 - x)^(-a) * (1 + x)^(-b) *diff((1 - x)^(a+n)* (1 + x)^(b+n),x$n);
 Lx:=unapply(pol,n,a,b);
 theta1(2) := Lx(2,0,0);

 

How can i expand the summation? i would like to have something like that:

evalf(theta1(2));
simplify(%);

How do I solve

int(int(exp(cos(theta)*sin(phi))*sin(phi), phi = 0 .. Pi), theta = 0 .. 2*Pi)

?

I think it is equal to

2*Pi*(int(sin(phi)*BesselI(0, sin(phi)), phi = 0 .. Pi))

but Maple also does not solve this integral?

 

Does anyone know how to solve these integrals (with or withouth Maple)?

Are there other programs that can solve these kinds of integrals?

(beta*p-beta-p)*alpha/(p*(beta-1)) reduces to alpha+beta/p assuming 0 < alpha, alpha < 1, 0 < beta, beta < 1, alpha+beta = 1. How do we show this with maple? Thanks.

u=[-2        v=[ 1           w=[3

      -5             7a               3

       5]             -2]              a]

Im newbie to Maple and trying to solve the following pde numerically:

PDE:=diff(u(t,z),t )-v*diff(u(t,z),z) = D*diff(u(t,z),z$2 z)z) -k*u(t,z);
ics := {D*`&PartialD;`(u(t, 0))/`&PartialD;`(z) = v*(C-u(t, 0)), D*`&PartialD;`(u)*(t, L)/`&PartialD;`(z) = 0, u(t, 0) = C};

sol := pdsolve(eval(PDE, {D = 1, k = 1, v = 1}), eval(ics, {C = 10, D = 1, L = 1, v = 1}), numeric, time = t, range = 0 .. 1);
Error, (in pdsolve/numeric/process_IBCs) improper op or subscript selector

I do not know how to correct this error. PLease help. Thank you.

I tried using 'maple tag' but the preview came out bad. So here's a clipboard copy of my maple code:

> f := proc (x, d) options operator, arrow; (1/4)*Pi*x^6-(3/5)*Pi*x^5+(1/60)*(30*x^3-72*x^2+80)*arctan(d/sqrt(x^2-d^2))*x^3+(2/3)*Pi*x^3+(1/60)*d*sqrt(x^2-d^2)*(16*d^4+x*(8*x-27)*d^2+2*x(3*x^3-9*x^2+20))-(1/60)*d(27*d^2-40)^3*log(x+sqrt(x^2-d^2)) end proc;
print(`output redirected...`); # input placeholder

I tried

taylor(f(x,d),x=0,d=0,4) but it didn't work. How do I do this?

I am having issues combining multiple graphs into one.  I have 2 sets of points that are defined as p1 and p2, and I have 2 functions defined as p3 and p4.  Here is the command & what it gives me when I try to graph them together...

> plots [display](p1, p2, p3, p4);

 

Hi everybody, probably a very basic question for most of you! I created a function with 2 lops n=1 to 10 and m=1 to 10. Each loop (n,m) I have a certain value for the local variable Xlocal. I want to know how I can put those value of Xlocal into a matrix. The place where I want to put it is (n,m) Matrix[ (1,1) (1,2)] [(2,1)(2,2)] Thanks a lot

graph the function and explain why the limit does not exist lim((x*y^3)/(x^2+y^6)) with (x,y)-(0,0) any helps

How can i find:

 

int(x*sinh(x)/sqrt(1-sinh(x)*sinh(x)), x = 0 .. log(sqrt(2)+1))

 

int(x*sinh(x)/sqrt(1-sinh(x)^2), x = 0 .. ln(sqrt(2)+1))

 

with maple

Hey, I'm using NonlinearFit to find parameterranges in theoretical equations for physical data. however some of the parameter values are being returned as negative values and I need to specify parameterranges but am having trouble doing this. 

For example in the equation below I would like to specify that alpha and ro should take on positive values, as otherwise with the curret data set I am using ro is assigned a negative value.

g1 := NonlinearFit(E(F, alpha, ro, .31, 165, e), X, Y, F)

investigate the family of functions f(x,y) = e^(c+x^2 + y^2). how does the shape of the graph depend on c? using maple

graph the surface, tangent line and normal line for:

xy + yz + xz = 3 (1,1,1). choose the domain so that avoid the extraneous vertical planes.