Hi,

I have a problem with maple solving diffusion equation. I solved the equation as follows

Es := 0.117108e12;
Ef := 0.78125e11;
l := 0.150e-6;
s := 0.500000e-3;
f := 0.5898334197e-6;
o := 0.9e-5;
d := 0.10e-17;
cb := 0.1e7/(19.9);
c := l*f/(d*cb);
PDE := diff(u(x, t), t)-(diff(u(x, t), x, x)) = 0;
IBC1 := {u(x, 0) = 0, (D[1](u))(0, t) = 0, (D[1](u))(1, t) = c};
S1 := pdsolve(PDE, IBC1, numeric, time = t);
S1:-plot(t = 0.6e-2);

i am plotting cincentration along the thickness of the material. i want to know why the concentration is negative. It should not be negative because constant flux 'c' is present at x=1. and flux at x=0 is 0, and also the initial concentration is 0.

Thanks

 

How can I plot the volume of revolution of the region between the curves y=ln(x) and y=-ln(x)on the interval [1,e]

around the y-axis.please specify the command. I used the command: VolumeOfRevolution(ln(x),-ln(x), x=1..e, scaling=constrained,axis=vertical,output=plot). But this command only plots the revolution of the curves not the region between them.

Best Regards

Yegan 

Hi,

I would like to plot E2(t) , but It gives errors. How can I avoid the singularity and solve this problem?

restart: with(plots):
with(Student[NumericalAnalysis]):
g1:=(x,t)->(-sqrt(t)/(2*sqrt(Pi*r^3)))*(sin((r*(x-1)^(2)/(4*t))+(Pi/4))-sin((r*(x+1)^(2)/(4*t))+(Pi/4))):

g2:=(x,t)->int(g1(x,t),r=1..infinity);

g3:=(x,t) -> (diff(g2(x,t),t)):

g4:=(x,t) -> (diff(g2(x,t),x,x)):

g5:=(x,t) -> ((1/2)*(g3(x,t)^2+g4(x,t)^2)):

E2:=t->(int(g5(x,t),x=0..100)):

evalf(E2(0));
Error, (in g1) numeric exception: division by zero
evalf(E2(1));
Error, (in g3) invalid input: diff received 1, which is not valid for its 2nd argument
plot(E2(t),t=0..20);

Best regards,

Hej guys,

in some of my Maple worksheets I'm using Maples fit-function to obtain a rational function fit to some data, as I'm interested in the zeros of enumerator and denominator, I chose my fit function accordingly. Afterwards I need to extract the results (mainly the zeros) and do some more calculation steps with them. When evaluating rational functions of high degree, I noticed that my way of extracting the results (I'm using patmatch) becomes increasingly slow up to the point where it takes a couple of minutes to extract 10 numbers.
E.g. the following minimal example (download link at the bottom) takes more than 20 minutes on my machine (Maple 18.01 Build ID 935137; Red Hat Version 6.6 running on system with an i5-3570, 8 GB of RAM). Additionally, the calculation requires an unreasonable amount of RAM.
I can work around the problem to some degree by extracting the leading coefficient beforehand, resulting in a substantial speed up, still I think that there might be some error in the code or my usage of patmatch.

Is there another way to extract these numbers in a similarly simple fashion?

Cheers,

Sören

Download patmatch-breakdown.mw

I'm trying to plot the region bounded by the coordinate planes, the plane z=2-y and the cylinder x=4-y^2

This is what I have so far.  I want to know how to get rid of the portions of the planes and cylinder that are not enclosing the region.

restart;

with(plots);

xmin := 0;
xmax := 4;
ymin := 0;
ymax := 2;
zmin := 0;
zmax := 2;

a := plot3d([2-y], x = xmin .. xmax, y = ymin .. ymax, orientation = [40, 70], transparency = 0, color = green, filled = false, title = "");
%;

b := implicitplot3d([z = 0, y = 0, x = 0], x = xmin .. xmax, y = ymin .. ymax, z = zmin .. zmax, transparency = .5, axes = normal);
%;
c := implicitplot3d([x = -y^2+4], x = xmin .. xmax, y = ymin .. ymax, z = zmin .. zmax, transparency = .5, filled = true, axes = normal, color = blue, orientation = [45, 35]);

display(a, b, c);

A duck, pursued by a fox, escapes to the center of a perfectly circular pond. The fox cannot swim, and the duck cannot take flight from the water. The fox is four times faster than the duck. Assuming the fox and duck pursue optimum strategies, is it possible for the duck to reach the edge of the pond and fly away without being eaten? If so, how?

http://www.crazyforcode.com/fox-duck-puzzle/

there is an animation here

https://www.youtube.com/watch?v=Zw9cHEnhzWo

wonder if the equations of motion can be derived usingg maple and an animaton...?

Hallo. I'd like to make a loop for two variables. Let say:

for i from 1..10, j from 1..10

How to do it?

Thanks.

Hi everyone:

I have two functions in terms of theta variable, how can I expressed function Y(X)?

X :=(theta)->cos(theta)+0.8e-1*cos(3.*theta)

Y :=(theta)->-sin(theta)+0.8e-1*sin(3.*theta)

I will earn Y(X) infact.

Thanks alot...

Mehran.

How to calculate the integral

(symbolically or/and numerically) with Maple?

How can I stop this text from appearing in the worksheet output?

What is the purpose of this text as it's a GIANT DISTRACTION from my perspective?

SAMPLE EXECUTED WORKSHEET OUTPUT:

> with(plots);

> with(geom3d);

outputs lots of blue text...

[Archimedean, AreCollinear, AreConcurrent, AreConjugate,

  AreCoplanar, AreDistinct, AreParallel, ArePerpendicular,

  AreSkewLines, DefinedAs, DirectionRatios, Equation, FindAngle,

  FixedPoint, GlideReflect, GlideReflection, GreatDodecahedron,

  GreatIcosahedron, GreatRhombicuboctahedron,

  GreatRhombiicosidodecahedron, GreatStellatedDodecahedron,

  HarmonicConjugate, HexakisIcosahedron, HexakisOctahedron,

  InRadius, IsArchimedean, IsEquilateral, IsFacetted, IsOnObject,

  IsQuasi, IsRegular, IsRightTriangle, IsStellated, IsTangent,

  MidRadius, NormalVector, OnSegment, ParallelVector,

  PentagonalHexacontahedron, PentagonalIcositetrahedron,

  PentakisDodecahedron, QuasiRegularPolyhedron, RadicalCenter,

  RadicalLine, RadicalPlane, RegularPolyhedron,

  RhombicDodecahedron, RhombicTriacontahedron, RotatoryReflect,

  RotatoryReflection, ScrewDisplace, ScrewDisplacement,

  SmallRhombicuboctahedron, SmallRhombiicosidodecahedron,

  SmallStellatedDodecahedron, SnubCube, SnubDodecahedron,

  StereographicProjection, StretchRotate, TangentPlane,

  TetrakisHexahedron, TrapezoidalHexecontahedron,

  TrapezoidalIcositetrahedron, TriakisIcosahedron,

  TriakisOctahedron, TriakisTetrahedron, TruncatedCuboctahedron,

  TruncatedDodecahedron, TruncatedHexahedron,

  TruncatedIcosahedron, TruncatedIcosidodecahedron,

  TruncatedOctahedron, TruncatedTetrahedron, altitude, area,

  center, centroid, coordinates, cube, cuboctahedron, detail,

  dilate, distance, dodecahedron, draw, dsegment, duality, faces,

  facet, form, gtetrahedron, hexahedron, homology, homothety,

  icosahedron, icosidodecahedron, incident, intersection,

  inverse, inversion, line, midpoint, octahedron, parallel,

  parallelepiped, plane, point, polar, pole, powerps, projection,

  radius, randpoint, reflect, reflection, rotate, rotation,

  schlafli, segment, sides, sphere, stellate, tetrahedron, tname,

  transform, translate, translation, transprod, triangle,

  vertices, volume, xcoord, xname, ycoord, yname, zcoord, zname] [animate, animate3d, animatecurve, arrow, changecoords,

  complexplot, complexplot3d, conformal, conformal3d,

  contourplot, contourplot3d, coordplot, coordplot3d,

  densityplot, display, dualaxisplot, fieldplot, fieldplot3d,

  gradplot, gradplot3d, implicitplot, implicitplot3d, inequal,

  interactive, interactiveparams, intersectplot, listcontplot,

  listcontplot3d, listdensityplot, listplot, listplot3d,

  loglogplot, logplot, matrixplot, multiple, odeplot, pareto,

  plotcompare, pointplot, pointplot3d, polarplot, polygonplot,

  polygonplot3d, polyhedra_supported, polyhedraplot, rootlocus,

  semilogplot, setcolors, setoptions, setoptions3d, spacecurve,

  sparsematrixplot, surfdata, textplot, textplot3d, tubeplot]

Hello those who are in this maple primes.

I downloaded maple2015 for Mac, the newest version of the software, and tried to use it to draw a simple figure with plot3d.

But, in the end, I used gnuplot which I had happened to install though it was my first time of using it.

The problem of using my Maple was that it was too heavy and with my macbook, I found it impossible

to have it work properly. So, I am writing my question now.

I used terminal, but as I think it  might be easy to deal with to use editor not terminal itself

though I use the editor from the terminal, I want to know how I can use maple from the editor and what things I 

can do when I do so. In my image, I want to use it like using TeX from the editor.

As for coloring the code, is the best tool one for Maple V yet? And, how can I do after I wrote the code of Maple?

What can I do with emacs and Maple? Can I use emacs like IDE of TeXshop or TeXworks in LateX?

I am glad to hear from you the general concept or procedures of using it if I use emacs and Maple.

Thanks in advance.

Dear All,

Colud you please tell me how to obtain the answer f(b)-f(a), by manipulating the following expression in maple:

Thank you very much. Best regards, Jand

Hello, I have a similary issue with the following post :

http://www.mapleprimes.com/questions/200835-Algsubs-With-Subscripts-And-Powers

In a constraint equation (in mechanics), I need to make these algebraic substitutions :

--> sqrt(xA1² + yA1²) =rF1 
--> arctan(abs(yA1)/abs(xA1)) = aF1
--> sqrt(xC1^2+yC1^2) = rBTP
--> arctan(abs(yC1)/abs(xC1))=a1

I try to do it with this code :

ListAlgsubs:={sqrt(xA1^2+yA1^2)=rF1,sqrt(xA2^2+yA2^2)=rF2,sqrt(xA3^2+yA3^2)=rF3, sqrt(xA4^2+yA4^2)=rF4,arctan(abs(yA1)/abs(xA1))=aF1,-arctan(abs(yA2)/abs(xA2))=aF2,pi+arctan(abs(yA3)/abs(xA3))=aF3,pi-arctan(abs(yA4)/abs(xA4))=aF4,sqrt(xC1^2+yC1^2)=rBTP,arctan(abs(yC1)/abs(xC1))=a1};
for i from 1 to nops(ListAlgsubs)
do Branch1Cons:=applyrule(op(i,ListAlgsubs),Branch1Cons)
end do:
Branch1Cons:= simplify(Branch1Cons);

I have tried with algsubs or applyrule but it wasn't successfull in both cases.

Do you have an idea how can I do to simplify my expressions ? and namely conduct a good use of applyrule function

I join an extract of my code in order to facilitate the troubleshooting.

example.mw

Thanks a lot for your help.

Dears, When I run calculation in Maple I found an error in matrices. See the file

Hi, I have 10 equation system and 10 unknown variables. I however, want to reduce the equations to 2 with two unknowns. I'm wondering how this could be done in maple. The variables are Y, q, y^d, y^*, y^x, H,  p^d, w, P and Pv. I intend to solve the equations for Y and w. 

Thanks in advance for your help. The maple file has also been attached.

Regards

Small_Open_Economy.mw