Hi,

I need to solve the heat diffusion equation in a system with two enclosed walls.  The first wall (made of a material called fr4) goes from x=0 o x=d1, and the second wall (made of copper, that is cu)from x=d1 to x=d2. 

u(x,t) is my unknown (change of temperature with time and space) for the first wall, and v(x,t) the unknown for the second wall.

So the system of PDEs to solve is:

>PDE:=[ diff(u(x,t),t) = afr4*diff(diff(u(x,t),t) ,

I am using Maple 12 on Windows XP.  When I export a plot to postscript as follows

plotsetup(ps,`filename.eps`,plotoptions=`width=2in,height=3in`);
plot(...);

the width and height options seem to have no effect.   The graphic is always 19.579 cm by 19.579 cm.  Other options such as portrait and noborders do have an effect.  Any suggestions for why this is not working?

Can anyone tell me how to work out correlation coefficients!!! I need it in very simple terms as I just dont get it ie if 2 variables has a correlation coefficient of .690 what does it mean?  Can someone explain it to me?  With thanks

I want a continous anti-derivate (cos(2*u*m)-alpha*u*sin(2*u*m))/(4*u^2+1)
w.r.t. u, for alpha and m Reals, 0 < u.

What I get with Maple in terms of Ci and Si jumps in m=0 (and for m=0 it
is not defined, due to Ci).

The best I got was

  Si(-conjugate(v))-Si(v); eval(%, v=m*(2*u+I));
  psi:= unapply(%, u,m);
  chi:= (u,m) -> piecewise(0 <= m, -Ci((2*u-I)*m)+Ci(m*(2*u+I)), 

Hi everybody:

Concerning exporting a worksheet as RTF (Maple12, win XP) I note two things: a) if you use inline output things such as equal sign, arrow with "at 5 digits" text, etc. appear slightly higher than the rest of the expression. also some part of the output is chopped off. b) the labels don't export.

Are these things consistent with what you could currently expect from maple or am I missing something?

thanks

Dimitrios

Suppose that the Hamiltonian is invariant under time reversal:

[H,T] = 0.

Show that, nevertheless, an eigenvalue of T is not a conserved quantity.

v/r,

1. How do I take a 2xN matrix and plot it on a 2D graph by its columns?

2.How do I take a number n and a 2xN matrix and return the product of A^n with the 2xN matrix, where A is a 2x2 matrix?

3.How do I, using your two functions above, write a procedure which plots on the same graph the result of multiplying the 2xN matrix by A^n for n = 0,1,2..10. How do i cause these sequence of plots to be displayed as an animation?

Hi,

I'm working on a computational physics project and prefer to use Fortran codes (compiled by MVS-9) with Maple power together.

In recent days I've changed some compiler options in MVS and unfortunatly (I don't remember the options) nothing works well as before.

For example, as I try to use a simple code like this:

This forum is the place for you to post all of your MapleSim questions.

Hi there,

So my problem is to find the solution(s) of a 3-D system. I can get Maple to spit out a solution (at times), but I am looking for reassurance that this is indeed the solution (if unique, and if not what the other solutions may be). Increasing Digits makes Maple "lose" the solution and return a blank-- normal behavior? Looking at the 3-D plot of the system offers something to puzzle over.

I set up the system in  (x,c,q):

xd := x^.5-.1*x-c;
cd := (.5/x^.5-.1-q)*c;
qd := -(q-.2e-1)*(q-.3e-1)+(q-.25e-1)*(.5/x^.5-.1-q)*c/(x^.5-.1*x-c);
 

I look for a solution with fsolve:

Digits:=10: interface(displayprecision=10):
ss:= fsolve({xd,cd/c,qd}, {x,c,q}, avoid={x=0, c=0, q=0});
xs:=subs(ss,x): cs:=subs(ss,c): qs:=subs(ss,q):

Maple gives:

ss := {c = 2.399289639, q = .2511089584e-1, x = 15.97164840}

If I raise digits to 15, Maple can't solve it anymore. Perhaps this is normal behavior.

I'd like some reassurance that this solution is indeed a solution. Simply feeding the fsolve output back into the system will not work because of a division by zero:

eval(xd,{x=xs,c=cs,q=qs}); eval(cd,{x=xs,c=cs,q=qs}); eval(qd,{x=xs,c=cs,q=qs});

yields:

0.
-.9597158556e-10
Float(-infinity)

I decided to look at the 3-D plot. The plot of the system suggests that there may be several solutions.

The plot:

plotopts:=style=patchcontour, shading=none,lightmodel=light3,axes=boxed:
px:=implicitplot3d({xd}, x=0..25,c=0..5,q=0.02..0.05, numpoints=10000, plotopts, colour=brown): pc:=implicitplot3d({cd/c}, x=0..25,c=0..5,q=0.02..0.05, numpoints=10000, plotopts, colour=green):
pq:=implicitplot3d({qd}, x=0..25,c=0..5,q=0.02..0.05, numpoints=10000, plotopts, colour=blue):
display3d({pq}, orientation=[-25,70]);
display3d({px,pc,pq}, orientation=[95,90], projection=1);

plot of third equation
plot of 3-D system

The plot is not very clean, and no matter how hard I've tried, by turning it around and zooming near the Maple solution, it is not clear that there are no other solutions. (well perhaps it is clear from the equations themselves, but I couldn't work it out and resorted to the numerical approach)

Does somebody have a sheet covering that algorithm for continous
anti-derivatives of rational functions (Bronstein's book §2.8
"Integration of transcendental functions")?

Or can show me, how to extract it from Maple's library in case
it is there (and I guess Maple covers that ...) for explicite
use?