please guide me how to write a simple code for Newton method to solve the following system of nonlinear equations
D= diagonal of A
L= lower triangular
U= upper Triangular
 

Download NM.mw
 

Download NM.mw

1) I want to visualize the probabilities of normal law. I used a technique ( a little heavy) ( see appendix) and I wonder if there is a way to do better ?

2 ) How do I calculate the probabilities of p( a<X<b) and p( X>c / X<d ) ?

3) How do I calculate the inverse probability in Maple. I only find the inverf(p) command under Maple TA?

NormaleVisualizationCourbe.mw

partial differential equation question??

u_tt = 9u_xx

_{boundary conditions : u (0,t) = (2,t) = 0}

_{initial conditions : u (x,0) = f(x) = { x(2-x)   ,0<x<1     ,      x   ,1<x<2  = 0 , t>0 ,}

                         u_t (x,0) = 0 , 0<x<2

I have a physics question I need to program it by Maple , Can you help me to solve ? Its Problem 9.34 from Griffiths

Hi,

I'm a beginner in programing with Maple.

Maple has an algorithm "LinearSolve". 

I'm trying - just for fun - to generate an algorithm to solve a linear system of equations according to the Gauss-Jordan-method.

The linear System is given by:

M1*x=b

(in the worksheet: M1 is a matrix, point (1)

b is a vector, point (2) and x is a vector too and is the solution of system, point (3)).

Point (3) is to check - after the generated algorithm points out the solution - if the algorthm is working fine.

The augmented matrix M2 is given by <M1|b>, point (4).

First, I want to generate zeros under the main diagonal of matrix M1 in the augmented matrix M2.

In the worksheet there is my idea to do this, but I get a (perhaps simple) error message ("Error, Invalid if statement termination"). But I don't know what I'm doing wrong...

Perhaps someone can have a look on my idea (link to the file) and give me a hint?

Thanks a lot for your help!!

LGS_Gauss-Jordan.mw

I am working on an Isoperimetric problem.

I have two différential dL and dA. And I am seeking the function y for wich the two différentials are collinear.

So I want "dL/dA=constant"

I tried many expression, but I don't know how to express the "= constant"
I tried also (dL/dA)'=0 but no workable answer...

Thank's for your help !

Hi!

I have the following problem

> restart:with(plots): Digits := 60:
> A:=42403:
> B:=269.95:
> C:=269.95:
> J1:=0.55:
> J2:=1.02:
> a:=196200:
> b:=14342220:
> c:=1589220:
> k:=157.08/0.0875:
> sys := {A-B*diff(x(t),t)=J1*diff(x(t),t,t)+C*(x(t)-y(t)),k*(J2*diff(y(t),t,t)-C*(x(t)-y(t)))=-(a+b*(diff(y(t),t)^2)/k^2+c*diff(y(t),t,t)),  D(x)(0) = 157.08, D(y)(0)=0, D(y)(10)=157.08,  D(x)(10) = 154.98}:
> dsn1 := dsolve(sys,numeric):
Error, (in dsolve/numeric/BVPSolve) matrix is singular

I would like to extract the unique values of a few sets in a list (L1) and insert them in another list (L2) without the sets.

L1:=[{3,5},{5,4}]:

L2:=[3,5,4]:

Could anyone help me with a bit of code so this works for all values in sets in a list?

I use the command 'labels = [x, Im(f)]'

to change the labels of my plots, but this time I need to put some "different" characters. i. e., \^{c} (in LaTeX notation).

How can I do this on Maple?

I need a help from someone who knows the GRTensor commands.

How do I write the last two terms of this equation of motion?How do I define Lagrangian matter, too?

My master said that I will do these calculations in the classic worksheet  maple18.

K² G_mn= (1/2) T_mn + (C₀/2) R_abgl T^ag T^bl g_mn+ 3 C₀ R_brs(mT_n)^b T^sr – C₀ L_m R_ab T^ab g_mn + C₀ R_ab T^ab T_mn+ 2 C₀  L_m R_mrns T^rs + C₀ ∇_a∇_r(T_(m^r T_n)^s)- C₀ ∇_s∇_r(T^sr T_mn)

Lm =-(1/4) F_mn F^mn 

Hi,

I can’t export a csv file.
The 'export' command returns a numeric value ?QuestionDataExport2.mw

How to sketch 3D graph finite string problem of wave equation in partial differential equation using maple?

u_tt = 4u_xx

_{boundary conditions : u (0,t) = (5,t) = 0 , t>0 ,}

_{initial conditions : u (x,0) = f(x) = { 0, 0<x<4      , (5-x) , 4<x<5} 

                         u_t (x,0) = 0 , 0<x<5

how to sketch 3D graph for solution of the corresponding partial differential equations ?? 

I thought in Maple the standard was to use _C1, and _C2, etc... for constants in the solutions returned.

Sometimes Maple mixes _C1 and c[2] in the same result. Is this common, to be expected sometimes and is OK? I noticed this only recently. 

I was thinking may be some part of Maple code still was not updated to use _C1 notation? Here is an example

restart;
pde:=diff(u(x,t),t)+ diff( u(x,t),x )^3 + 6 * u(x,t)* diff(u(x,t),x) = 0;
sol:=pdsolve(pde,u(x,t));

which gives

sol := u(x, t) = -(3/2)*_C1^2+3*(t*_c[2]+x)*_C1-(3/2)*(t*_c[2]+x)^2-(1/6)*_c[2]

With latest Physics updates  268

Is there an option, like AllSolutions used with solve, so that pdsolve would return all solutions to a PDE when it is nonlinear?

I looked at pdsolve help and do not see a HINT that looks like might do this.

For example, this PDE, Maple returns one solution. But Mathematica returns 2 solutions

restart;
pde:= diff(u(x,t),t) = diff(u(x,t),x$5)+10*diff(u(x,t),x$3)*u(x,t)+25*diff(u(x,t),x$2)*diff(u(x,t),x)+
             20*u(x,t)^2*diff(u(x,t),x);
sol:=pdsolve(pde,u(x,t));

#sol := u(x, t) = -12*tanh(176*_C2^5*t+_C2*x+_C1)^2*_C2^2+8*_C2^2

But there is another solution

sol1:=u(x,t)=-(1/2)* _C1^2*(-2 + 3*tanh(x*_C1+ t*_C1^5 + _C2)^2)
pdetest(sol1,pde)
#0

Here is another example. Maple returns one solution and Mathematica 7 solutions

restart;
pde:= diff(u(x,t),t)= u(x,t)*(1-u(x,t))+ diff(u(x,t),x$2);
sol:=pdsolve(pde,u(x,t));

#sol := u(x, t) = (1/4)*tanh(-5*t*(1/12)+(1/12)*sqrt(6)*x+_C1)^2-
              (1/2)*tanh(-5*t*(1/12)+(1/12)*sqrt(6)*x+_C1)+1/4

But there are other solutions

pde = D[u[x, t], t] == u[x, t] (1 - u[x, t]) + D[u[x, t], {x, 2}];
DSolve[pde, u[x, t], {x, t}]

I've tested some (not all) of these 7 solutions in Maple using pdetest and Maple agrees they are solutions:

restart;
pde:= diff(u(x,t),t)= u(x,t)*(1-u(x,t))+ diff(u(x,t),x$2);
sol:=pdsolve(pde,u(x,t));
with(MmaTranslator);
sol2:=FromMma(`-(1/4) (-3 + Tanh[(5 t)/12 - (I x)/(2 Sqrt[6]) - C[3]]) (1 + 
   Tanh[(5 t)/12 - (I x)/(2 Sqrt[6]) - C[3]])`);
pdetest(u(x,t)=sol2,pde);
#0

I tried setting 

       _AllSolutions:=true

But it had no effect. Is there other options?