MaplePrimes Questions

Hi there,
I have a simple task to do with arrays. I have an 8x3000 array, whose columns are sometimes all zeroes (if any element in the column is zero then the entire column is zero). I want to eliminate the zero vectors, so I have this loop (the array is A)

for i from 1 to 3000 do
    while A(1,i)=0 do
        A:=DeleteColumn(A,i);
    end do;
end do:

Setting "A:=DeleteColumn(A,i);" reduces the size of the array, so if all the zero columns have been deleted at "i=200", Maple will keep trying to check the next column, which doesn't exist anymore. This gives an "index out of bounds" error, but also records "A" as the new smaller array, which is fine.

If I run this loop inside a procedure, the error stops the procedure, but it doesn't record the new array, it stays as the older, larger one. Does anyone know of a way around this?

Any help would be much appreciated.

Below and above the iterators, there is a graphic with numbers rotated in all directions.  I am thinking it is based somewhat on what I have done here http://www.mapleprimes.com/posts/97195-Wordle-In-Maple

How did they get the rotations in Maple?

 

 

In functions such as fsolve, there is an optional parameter that allows one to specify the interval to perform the function on. Additionally, sometimes, if a solution is left out, one can specify an interval to search on to obtain the missing solution.

How does Maple determine the interval to search on if this is not specified?

 

(Additional Question - you don't have to answer this)

Ultimately, I am asking this question because I have a function for Newton's method; however, it requires an interval to run. I have read that fsolve uses Newton's method, so I am curious how to automatically select such an interval. Does anyone know how to implement such a thing?

Hi,

 

Let, fixed an integer i and 1<=j<=2^{i}-1, for each x and y in [0,1] let the following mapping

Then, with the above procedure we can obtained, for a fixed i, all the mappings for j=1,...,2^{i}-1

 

However, How can I to evalute the "components" of the above procedure? For instance, I can not to compute CreaF(2)[1](0.35,0.465) (i.e., the first function in the "vector" CreaF(2), in x=0.35, y=0.465). 

 

Thanks very much for your time.

 

hi everyone..i would like to ask,why do we need to type in (101-100λ) in our maple program..why do we need the continuation..thanks in advance for answering..

how to convert baby monster group into permutation group and iterate baby monster to display all permutation group?

if not convert to permutation group directly, 

how to convert to matrix group first and then convert to permutation group?

 

Hello,

I have a sequence of functions :  solution[i] , i = 1..n

I have a sequence of times:          Time[i], i = 1..n+1

I need help to plot in the same graph:

 plot(  Sol[1] , t = Time[1]..Time[2] )  

 plot(  Sol[2] , t = Time[2]..Time[3] )  

 plot(  Sol[3] , t = Time[3]..Time[4] )  

etc...

plot(  Sol[n] , t = Time[n]..Time[n+1] )  

Thank you

 

 

 

This question has been deleted by the author

 

Hello,

 

I got two problems creating a 3d-plot. First, as the title mentions, I want to cut out a distinct interval from my 3d-Plot.

For example when the graph is displayd from -y_max to y_max but I just want it to be shown from -y_max to -5 and from 5 to y_max again.

I also got a second question:

Another thing I want to do with this graph is redefining this axis with a parameter which is directly dependant from it.

Like for example I want to redefine the "y" axis with a" k" axis with k=y*2, so that in the end I have th

Hi 

Sometimes this happen to me.

why ?  and how can I fix it?

Please recover the content of the attached file if it is possible

Thanks

Check.mw

Hi dear experts

I'm wrote a complicated code for stress computation and define each component of stresses for n layers in a n dimension array.

for example (for srr component):  srr[1]:=f1(r,theta) , srr[2]=f2... , ..., srr[n]=fn(r,theta)

n is defined in the start of worksheet and is no. of layers!

so, i want to plot all of layer functions in a same coordinate plane.( beacause they are discontinuse and i want show that)

note: display command don't do that and plot them in a table form! :(

how can i do that?

 

thanks for participate.

hi.i have a problem with gain Eigenvectors

please help me

thanks

lorenz_1.mw

restart; with(LinearAlgebra[Modular]); A := Matrix(3, 3, [[-a, a, 0], [0, 0, -sqrt(l.b.c.(h+k))/(h+k)], [(2.*h.sqrt(l.b.c.(h+k)))/(l.(h+k)), (2.*k.sqrt(l.b.c.(h+k)))/(l.(h+k)), -c]])

A := Matrix(3, 3, {(1, 1) = -a, (1, 2) = a, (1, 3) = 0, (2, 1) = 0, (2, 2) = 0, (2, 3) = -(`.`(l, b, c, h+k))^(1/2)/(h+k), (3, 1) = 2.*(h.((`.`(l, b, c, h+k))^(1/2)))/(l.(h+k)), (3, 2) = 2.*(k.((`.`(l, b, c, h+k))^(1/2)))/(l.(h+k)), (3, 3) = -c})

(1)

Sol := LinearAlgebra[Eigenvectors](A)

Error, (in evala/preproc3) floats not handled yet

 

``


Download lorenz_1.mw

My question is related to changing the plot size and modifying the number of ticks.

For example, I am plotting y vs. x, and x is from -1.0*Pi/2..1.0*Pi/2. I want to modify the plot size with size=[h,w]. If I have three ticks in x-axis if I change to small width, then the number of ticks is decreased. How can I force it to have the certain number of ticks? Also, why if I set view=[-1.0*Pi/2..1.0*Pi/2,0..1] then the -1.0*Pi/2 disappear if the width size is small. 

Thanks,

Baharm31

Hello people in the mapleprimes,

I have a question, so I hope someone give me answers to it.

I calculated for the solution of the follwing differential equation.

restart
b:=diff(y(x),x)+a*y(x)=f(x);#where a and f(x) is not specified.
dsolve(b,y(x));

subs({f(x)=exp(x),a=2},%);where f(x) and a are specified.

c:=value(%);

The solution of the above was

y(x) = (1/3)*exp(x)+_C1/(exp(x))^2,  (A)

where please note that the second term takes

the form of fraction _C/(exp(x))^2.

 

On the other hand, next I calculated the following differential equation where f(x) and a are specified from the start.

restart
d:=diff(y(x),x)+2*y(x)=exp(x);

dsolve(b,y(x));

Then,

y(x) = (1/3)*exp(x)+exp(-2*x)*_C1  (B)

was the obtained solution.

 

Each (A) and (B) are the same substantially mathematically. But, for Maple, the variable powered to minus brabra

is not the same as one over variable powered to brabra, so that (A) and (B) takes different forms, and maple will see them 

different with each other.

 

  Surely, with algsubs, algsubs(_C1/(exp(x))^2=exp(-2*x)*_C1,c) transforms (A) to (B).

But, I want to know whether there are some other ways than that  to modify (A) to (B).

If there are any good ways for it, I will be happy if you teach them to me.

Thanks in advance.

 

taro

First 1111 1112 1113 1114 1115 1116 1117 Last Page 1113 of 2434