I'm trying to solve a system of four pdes and I know that the Newton method won't converge.

Are there other numerical methods that I can use?

Any help would be greatly appreciated!

Thanks,

Eve

Hi,

I'm an industrial engineering student who's running into problems with solving a simple system of non-linear equations.

Why do I get this error when issuing the solve command? I'm pretty sure I'm only passing two lists into the function?

Thanks a lot in advance,

Joshua

Hi!

Is there any way to remove the empty space that comes under images when printing the project, while using document mode?

Regards

Nicolai

Please fix the MapleCloud service. Is no longer working. E uploaded jobs and are not shown. Fix it please.

I have a matrix A for which the basis of the left null space using NullSpace() is the empty set {}  while the column space is {e1,e2,e3}. By definition, we need every vector in col space . every vector in basis of left null space =0 but how would I show that in this case? Can I determine another basis for the left null space?

Determine using determinants the range of values of a (if any) such that
f(x,y,z)=4x^2+y^2+2z^2+2axy-4xz+2yz
has a minimum at (0,0,0).

From the theory, I understand that if the matrix corresponding to the coefficients of the function is positive definite, the function has a local min at the point. But, how do I get the range of values of a such that f is a min? Is this equivalent to finding a such that det(A) > 0?

2.

Now modify the function to also involve a parameter b: g(x,y,z)=bx^2+2axy+by^2+4xz-2a^2yz+2bz^2. We determine conditions on a and b such that g has a minimum at (0,0,0).
By plotting each determinant (using implicitplot perhaps, we can identify the region in the (a,b) plane where g has a local minimum.

Which region corresponds to a local minimum?

Now determine region(s) in the (a,b) plane where g has a local maximum.

I don't understand this part at all..

I've got the following matrix :

A:=[<a,a-1,-b>|<a-1,a,-b>,<b,b,2a-1>] where <> are the column elements of A, a is  a real number defined on [0,1] and b^2=2a(1-a) 

a) to show A is an orthogonal matrix, I understand that I need A.Transpose(A)=Identity(3*3) but is there a way in which I can let a take a random real numbered value between 0 and 1? The rand() only returns an integer within a range. Directly multiplying A and Transpose(A) will return an expression in a, so what's the right approach?

b) from a) we can infer that A is a matrix that describes a rotation in e1,e2,e3 where these are the standard bases vectors in R3. How can I determine the rotation axis? The hint I've been given says I need to consider the Eigenvalues and eigen vectors but I don't quite understand how.

I have the function:   f(x,y) = 1/(sqrt(2*Pi)) * e^{-1/2(x^2+y^2)}

I need to take the total derivative of this function in maple, but I don't know the syntax.

So i got a procedure test, she is kind of numeric, i whant to optimaze test([.5, .5, .5], 1, 3, 100, 100, true, [x, 0, 0, 0, 0])=0 by x. But optinization substitutes x like a symbol, i tryed all methods but they all do the same.

f := proc (x) options operator, arrow; abs(test([.5, .5, .5], 1, 3, 100, 100, true, [x, 0, 0, 0, 0])-.4) end proc; Minimize(f(x));
%;
Error, (in test) cannot determine if this expression is true or false: 0 < -43.0+100*x

Can i some how use optinization on such procedure?

file link  - >    primset.mw

A created an overly complicated worksheet thaI can no longer open.  It's about 17 Mb, has several interactive plots and equations in it.  Is there a way to force Maple to open a worksheet without executing it?

(Lesson learned:  split this up into separate worksheets next time.)

I wish to use closed Newton-Cotes with n=2, also known as Simpson's Rule to numerically integrate an improper integral.

If it matters the integrand is (cos(2x))/(x^1/3), integrating between x=0..1

I've tried a few different (but similar) code but to no avail. Here is some stuff I've tried:

1.

with(Student[NumericalAnalysis]):

with(Student[Calculus 1]):

Simp1 := ApproximateInt(cos(2*x)/x^(1/3), x = 0 .. 1, method = newtoncotes[2]);

This gives me an output message that says "Float(infinity)".

2.

with(Student[NumericalAnalysis]):

with(Student[Calculus 1]):

Simp2 := int(exp(-x)/sqrt(1-x), x = 0 .. 1);

This doesn't have Simpson's rule as an option.

I think I'm on the right track with my first try, since I guess it wasn't tecnically an error message, but I'm not sure how to alter the code accordingly to get a numerical value instead. Thanks for any help.

Thanks.

Download dsolve_question.mw

Hi,

I got the Real and Imaginary of an expression J1 

assume(d,real):

Gamma:=0.04:tau:=10*Pi:j:=0:

J1:=(exp((1-I*d)*Gamma*tau)-1)/((1-I*d));

J1mod:=simplify((Re(J1))^2+(Im(J1))^2): (I works here this amont is real)

################

but when I change the expression  for J1 to be

J1:=((2*e^(-2^(-j-1)*(1-I*d))-e^(-2^(-j)*(1-I*d))-1)*exp((1-I*d)*Gamma*tau)-1)/((1-I*d)):

J1mod:=simplify((Re(J1))^2+(Im(J1))^2): 

J1mod here is complex(I dont know why? it doesnt separate the real and the im )

Any comments will help

Thanks

Hello everybody,

can somebody tell me how to use subscripts in textmode? 

Thanks for your help!

Picard

OK, what am I doing wrong here?

This polynomial can be factored by almost anyone who knows algebra:

factor(x^2+y^2+2*xy);

but Maple refuses:

 x²  + y²  + 2 xy

but if I pose the question this way:

ans := expand((x+y)^2);

factor(ans);

(x + y)²

I get the expected result.  What is wrong with how I am using the factor?