I want to write a loop, which can give me the appropriate roots of a long equation.

I have tired to use a very simple loop, but the problem is that there are multiple solutions of the equation.

I want a loop which use newton methods to give me the satisfactory roots.

The funtion which I am dealing with is of the form sigma(l).

Plz help me out.

I have attached the maple file.

Hi everyone,

when I use:

c:=2e-6*Unit(A*s/V);

then maple answers with the unit [F]. but when I use:

epsilon[0]:=8.85e-12*Unit(A*s/V/m);

then in the middle of a large document, Maple (Version 12) ignores [m] and still gives me [F]. On a blank (new) document the same command results in correct units [m^-3 kg^-1 s^4 A^2]. The behaviour is reproducible.

anyone knows why?

cheers :-)

 

 

I want to graph a real-valued function e.g. plot y=x^2 from x=1 to 2.  There used to be a function called "plot2D". But it does not exist any more. When I call up the glossary or Help list, I get hundreds of plot-related functions - countourplot, loglogplot, etc - but not a plot for doing a real-valued function.  It is not under with(plots): 

In Advanced Engineering Mathematics there are a lot of math symbols within the text that don't execute. They're not simply an inert math statement, like Int or Diff, but rather the symbolic math used solely as text. How is this done?

Can I store data in a 3D Matrix in Maple ?

 

I want to solve the following PDE with boundary conditions, but unfortunately I do not know the meaning of _F1 and _F2 in this stuation. Would you please help me to find these arbitrary functions?

> PDE := diff(U(y, z), y, y)+diff(U(y, z), z, z) = 0;

> with(PDEtools);

Hi,all, I found a weird problem. When I use alias of RootOf in a procedure, it seem that it cannot work correctly. The code is as following: > restart; > cf := proc (n, ir) local P; alias(x = RootOf(ir)); P := evala(Expand(x^4)) mod 2; print(P); end; >cf(4,x^4+x+1); x^4 ---------------------> should be x+1 However, when I use it in seperate commands, it works fine. > restart; >alias(x = RootOf(x^4+x+1)); >P := evala(Expand(x^4)) mod 2; >print(P); x+1 I have no idea what is wrong here Thanks Gepo

The following program produces a plot in the 4th quadrant.

S1:=unapply(u*FresnelS(u)+((1/Pi)*cos((1/2)*Pi*u^2))-(1/Pi),u);

C1:=unapply(u*FresnelC(u)-((1/Pi)*sin((1/2)*Pi*u^2)),u);

R1:=1;

u1:=evalf(sqrt(0.4));

E1:=evalf(Pi*R1*S1(u1)+R1);

D1:=evalf(Pi*R1*C1(u1));

with(plottools);

with(plots);

c1 := circle([D1, 0], R1);

a:=evalf(Pi*R1*u1);

p1:=plots[display](plot([1+a*FresnelC(u),a*FresnelS(u)-E1,u=0..u1]),c1,scaling=constrained);

p3:=plot(0,x=-1..0);

display(p1,p3);

 

Download 16210_compressed form.mws
View file details

 

My friends, I am trying to compress a very long equation. Can any one help me to make it possible to compress the equation, means to write it in a very very simplified form. The real reason for this hardwork is that the equation is very long to write it in latex.

Your quick response is required.

Thanks to all

Hi everybody and thanks a lot for your help!

My problem is the following: I would like to draw a phase diagram for a system of 3 differential equations. So far this is not a problem but I would like to have the arrows of the vector field included in the diagram, like it is possible for systems of 2 diff. equations. Is that possible?

I use DEPlot3d, I suppose that is the correct tool?

Thanks and regards

Hans

I have the above program

S1:=unapply(u*FresnelS(u)+((1/Pi)*cos((1/2)*Pi*u^2))-(1/Pi),u);

C1:=unapply(u*FresnelC(u)-((1/Pi)*sin((1/2)*Pi*u^2)),u);

R1:=1;

u1:=evalf(sqrt(0.4));

E1:=evalf(Pi*R1*S1(u1)+R1);

D1:=evalf(Pi*R1*C1(u1));

with(plottools);

with(plots);

c1 := circle([D1, E1], R1);

a:=evalf(Pi*R1*u1);

plot([a*FresnelC(u),a*FresnelS(u),u=0..u1]);

p1:=plots[display](plot([a*FresnelC(u),a*FresnelS(u),u=0..u1]),c1,scaling=constrained);

p3:=plot(0,x=-1..0);

display(p1,p3);

I have found eigenvalues of a  2x2 matrix A, these are displayed in a vector. eg.

a + sqrt(b)

a - sqrt(b)

Is there any way I can write code which will find the difference between them? At the moment I'm doing it manually and it is taking a lot of time as I have many matrices.

Hello
Yesterday I was playing around with the Tolerances package in Maple 12. To familiarize myself I with it I tried to work out the equivalent resistance of two parallel resistors with tolerances.
R1:=80.0 +/- 4.0 [Ohm]
R2:=120.0 +/- 6.0 [Ohm]
Maple returned: Req:=R1*R2/(R1+R2) = 48.5 +/- 7.22 [Ohm]   (rounded values, there is some weirdness going on with this also.)
The right answer should be 48 +/- 2.4 [Ohm].

Dear Maple users

I read somewhere that plots with units is new in Maple 13. There is however something which is not working well:

> restart:
> T:=30 [[year]]:
> A0:=370000[[Bq]]:
> A(t):=A0*(1/2)^(t/T):
> plot(A(t), t=0..100[[year]])

I expected to see a plot of an exponential with units. A coordinate system was displayed but no function plot!! Can someone explain why? If I remove the units, the exponential appears!

Regards, Erik V. 

Download 16210_convection1.mws
View file details

I am try to write the solution in latex but its too long. Can any one help me to use some maple commads to write the solution in a compressed form.

Thanks