Pls, I tried solving the system of PDE numerically.... When I did for just 1 plot, the graph was plotted easilyrunning.mw but When I varied some parameters its not coming out.... Also, its not bringing any error so I can't trace out my mistakes....Not_runnin.mw PLS HELP ME OUT with the multiple plots... Attached are my source codes links....

When using Maximize, is it possible to now what parameter settings that resulted in complex numbers?

KR Fredrik

Maximize(Eq1_new, {R_K*rho*`&sigma;_K`+R_S*`&sigma;_S` >= .38780613588076011954, r_K*rho*`&sigma;_K`+r_S*`&sigma;_S` >= .21290471165685214846, R_K*rho*`&sigma;_K`+R_S*`&sigma;_S` <= .38780613588076011954, r_K*rho*`&sigma;_K`+r_S*`&sigma;_S` <= .21290471165685214846});

Error, (in Optimization:-NLPSolve) complex value encountered

intsolve(int(f1(u)*(x-y), u = a .. b) = (12)*cos(x)+5*x^2-12, f1(x));I try to resolve the integral equation but i received all the time the followinf message error 

Error, (in intsolve) invalid input: rhs received {f1(x) = ?}, which is not valid for its 1st argument, expr
I try to use different form but i do not how to resolve the help center do not help to much.

I preciate any help thanks

Hi.

please, I need some code about stepsize adaptative Runge Kutta method ( any order ) or other methods.

Thanks

i cant find the error the program is saying i have 

> L := -1.576674; MU := 0; DE13 := {D(y)*t = -x(t)*(1-2*x(t)*x(t))+MU*y(t)*(x(t)*x(t)-3*y(t)*y(t)-L), (D(x))(t) = y(t)*(1-y(t)*y(t))+MU*y(t)*(x(t)*x(t)-3*y(t)*y(t)-L)}; DEplot(DE13, [x(t), y(t)], t = -20 .. 20, [[x(0) = 0.1e-1, y(0) = .99], [x(0) = 0.5e-1, y(0) = .95], [x(0) = .1, y(0) = .9], [x(0) = 0.4e-1, y(0) = .96]], stepsize = 0.1e-2, scene = [x(t), y(t)], title = "phaseplane plot", linecolor = black, thickness = 1, number = 1000);
 
Error, (in DEtools/DEplot/CheckInitial) too few initial conditions: [x(0) = 0.1e-1, y(0) = .99]

it might be hard to read but if someone could help me it would be very appreciated 

Hi

I need a temperature distribution inside a barrier during a heating process.
I will be appreciated for any help.

wz

Download heating.mw

restart:

Eq1:=1/r*diff(r*diff(w(r),r),r)-(cos(z1))^2*(w(r)+1)-DP;

dsolve({Eq1});

bcs:=D(w)(0)=0,w(r2)=A1;

dsolve({Eq1,bcs},w(r)):

where z1, Dp, r2, A1 are constants.

Hi, i'm wondering if theres anything i can do printf and print and have the outputs on the same line. I want the string wihtout any quotations so i think i can only do printf to get this, but then i'm putting an entry out of an array after this. This is a small example out of my code, i would like the output to be all on one line:

The synchronising vector is (S[j-1]) The synchronising word is (X[j-1]) The length of the synchronising word is ...

if (sync=true) then 
printf("The Synchronising vector is"); print(S[j-1]);
printf("The sychronising word is");
print(X[j-1]);
printf("The length of the synchronising word is");
print(length(X[j-1]));
else printf("There is no synchronising word"); 
end if;

 At the moment its outputing each print on a new line.

Any help would be much appreciated, thanks.

hi all.
i have a system of ODE's including 9 set of coupled OED's . 

for example

a*b v a^-1 = b

i guess Disj or Conj are Max or Min respectively

however i do not know how to max(a,b) where a and b are permutation group

reference from L group in 

http://books.google.com.hk/books/about/Lattice_Ordered_Groups.html?id=_qAFX38-pKkC&redir_esc=y

if can not calculate this, do not know how to determine whether equal in a*b v a^-1 = b

I came to know that Chrome no longer supports MathML. Will this be an issue with Maple T.A.9/9.5? Should I require my students to switch to Firefox instread? 

One annoying feature in using document mode that I found is this. As I hit the return key to execute one command and automatically jump to the command below it (a feature that I like), but after few returns, I find myself at the bottom edge of the window, with the cursor pointing at the last command in the window.

Now when hitting return, the current command is execuated, and the next command below it appears, and it is also at the bottom edge of the window.  The cursor is from now on, always sitting at the bottom edge of the window, and I am not able to fully see the rest of the command and what is below it any more.

What would be better, is to try to keep the window centered, so that as I hit a return, the actual cursor remains in the middle of the window, and the window itself scrolls up. This way the command I am about to evaluate, and anything below it, remain in full view.

What I end up doing now, is when I get to the bottom of the window, is to manually scroll the window down using the mouse, so that now the cursor is up higher, and I can see more the commands which I am about to evaluate, then I repeate this process.

Is there  a way to configure Maple to try to keep the window centered? it is really very annoying to have to keep adjusting the window every few returns, and when one has 100 commands to evaluate, one by one (since I want to see the result before I go to the next command), this process gets tiring.

(does Maple do any usability studies and evaluation on its user interface by any independent outside group?)

Hi!

I have the following function:

and I want to calculate it for a certain set of r, theta and t, but when I use subs, the theta variable is also subsituted int the integral's parameters and I get the folllowing result:

I really don't understand why it does that given that there is no issue for the r variable...

Thanks in advance for your help!

P-S: please forgive my poor English : I'm a French student...

Florent

Dear all

Please I need your help to simplify by the coefficient a in this Matrix

I have The matrix A defined by  A:=Matrix(2,2, [[a,a],[3*a,4*a]]);

I want with maple transform A to  A:=a*Matrix(2,2, [[1,1],[3,4]]);

Thanks for your Help.

Dear all;

Please see only the last lines of this procedure,

1)I ask if this procedure give an output epsilon(x0,h). Really, I need your help. Thanks.

2) Can we plot epsilon(x0,h), versus h

> heun := proc (x0, h)

local x, y, i, N, k, f, ode, k1, k2, x1, x2, y1, y2;

f :=(x,y)-> 1/(1+cos(y)) ;

ode := diff(y(x), x) = f(x, y);

N := round((1/2)*x0/h);

y2 := Array(0 .. N);

x2 := Array(0 .. N);

y1 := Array(0 .. 2*N);

x1 := Array(0 .. 2*N);

x2[0] := 0; y2[0] := (1/4)*Pi;

for i from 0 to N-1 do

x1[2*i+2] := (2*i+2)*h;  k1[1] := f(x1[2*i], y1[2*i]); 

k1[2] := f(x1[2*i]+h, y1[2*i]+h*k1[1]);

y1[2*i+2] := y1[2*i]+(1/2)*h*(k1[1]+k1[2]);

x2[i+1] := (2*i+2)*h;

k2[1] := f(x2[i], y2[i]);

k2[2] := f(x2[i]+2*h, y2[i]+2*h*k2[1]);

y2[i+1] := y2[i]+h*(k2[1]+k2[2])

end do;

return firstresulat = evalf([seq([x2[i], y2[i]], i = 0 .. N)]);

return secondresulat =evalf([seq([x1[2*i], y1[2*i]], i = 0 .. N)]);

epsilon:=(x0,h)->add((firsttest[i][2]-secondtest[i][2])^2 , i=1..round(x0/(2*h))+1)

end proc;