7 years, 45 days

## MaplePrimes Activity

### These are replies submitted by mehdibgh

@dharr Thanks for your response. What about the eigenvectors? As you can see I expect to have 2000*2000 eigenvector matrix since the size of X is 2000*1, but this method only returns 1200!!!

How to discern my eigenvectors???How to plot my eigenvectors???

## Is not there other way to deal with thes...

@dharr Thanks for your knd answer. I am wondering if there are another ways to find the eigenvalues and eigenvectors of such eigenvalue problems.

It will be appreciated if somebody suggest another way.

There is also another using qr decomposition, which is also time consuming. Here I am attaching it.

 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >

## I provided relevat data for your attenti...

@mehdibaghaee I am adding the relevant matrices as well as worksheet as you requeste to present.

For specific problem the eigenvalues are as below:
1.5477,  2.2752,  2.9762,  3.3202 , ...

I have tried different ways in Maple to find them but unfortunately was not accurate enough and also time consuming.

The only way I could find them is by changing the variable (LV-Method) as below which gave me the nearest solutions.

What do you suggest to find the eigenvalue and eigen vectors of this problem fast and accurately in Maple (other than LV-Method)?

Change the txt to mat in the following files. This is because I could not upload mat format.

Ks.txt

Ms.txt

## @mmcdara i HOPE this problem will be sol...

@mmcdara i HOPE this problem will be solved in new versions of maple, where we can use 2D math safely.

## @Kitonum Why these problem occured and h...

@Kitonum Why these problem occured and how to solve it? I just typed in Maple

## @mmcdara Thanks. I am looking for the po...

@mmcdara Thanks. I am looking for the possiblity to plot with only one plot3d command in one line short command?

## @Rouben Rostamian excuse there were nota...

@Rouben Rostamian excuse there were notation errors in my equations which decoupled them. Here I corrected them as below:

And the example is:

Actually, I know these kinds of equations do not have symbolic solutions, I am looking for the best way to solve them by Maple numerically.

## @Rouben Rostamian  Actually I want ...

@Rouben Rostamian  Actually I want to solve the following differential equation set:

for example:

Do you know how to do this in Maple?

## @Carl Love your commands give error: pl...

@Carl Love your commands give error:

plot([seq(PPP[[1,k],..]^+, k= 2..rhs(rtable_dims(PPP)[1]))], color= black);
Error, invalid power

## @Kitonum Thanks, Actually I dont have fu...

@Kitonum Thanks, Actually I dont have function, I have many vectors, where I want to plot all in one figure. How to plot all of them in one fig with only one line of commands, as you did above for functions.

It is tedious to plot one by one and then gather them as below

with(plots); g := 2926;
P1 := plot(PPP[1, 1 .. g], PPP[2, 1 .. g], style = line, color = black);

P2 := plot(PPP[1, 1 .. g], PPP[3, 1 .. g], style = line, color = black);

P3 := plot(PPP[1, 1 .. g], PPP[4, 1 .. g], style = line, color = black);

P4 := plot(PPP[1, 1 .. g], PPP[5, 1 .. g], style = line, color = black);

P5 := plot(PPP[1, 1 .. g], PPP[6, 1 .. g], style = line, color = black);

P6 := plot(PPP[1, 1 .. g], PPP[7, 1 .. g], style = line, color = black);

P7 := plot(PPP[1, 1 .. g], PPP[8, 1 .. g], style = line, color = black);

P8 := plot(PPP[1, 1 .. g], PPP[9, 1 .. g], style = line, color = black);

P9 := plot(PPP[1, 1 .. g], PPP[10, 1 .. g], style = line, color = black);

P10 := plot(PPP[1, 1 .. g], PPP[11, 1 .. g], style = line, color = black);

P11 := plot(PPP[1, 1 .. g], PPP[12, 1 .. g], style = line, color = black);

P12 := plot(PPP[1, 1 .. g], PPP[13, 1 .. g], style = line, color = black);

P13 := plot(PPP[1, 1 .. g], PPP[14, 1 .. g], style = line, color = black);

P14 := plot(PPP[1, 1 .. g], PPP[15, 1 .. g], style = line, color = black);

P15 := plot(PPP[1, 1 .. g], PPP[16, 1 .. g], style = line, color = black);

P16 := plot(PPP[1, 1 .. g], PPP[17, 1 .. g], style = line, color = black);

P17 := plot(PPP[1, 1 .. g], PPP[18, 1 .. g], style = line, color = black);

P18 := plot(PPP[1, 1 .. g], PPP[19, 1 .. g], style = line, color = black);

P19 := plot(PPP[1, 1 .. g], PPP[20, 1 .. g], style = line, color = black);

P20 := plot(PPP[1, 1 .. g], PPP[21, 1 .. g], style = line, color = black);

P21 := plot(PPP[1, 1 .. g], PPP[22, 1 .. g], style = line, color = black);

P22 := plot(PPP[1, 1 .. g], PPP[23, 1 .. g], style = line, color = black);

P23 := plot(PPP[1, 1 .. g], PPP[24, 1 .. g], style = line, color = black);

P24 := plot(PPP[1, 1 .. g], PPP[25, 1 .. g], style = line, color = black);

plots[display]([P1, P2, P3, P4, P5, P6, P7, P8, P9, P10, P11, P12, P13, P14, P15, P16, P17, P18, P19, P20, P21, P22, P23, P24]);

## @Kitonum What about the font size of the...

@Kitonum What about the font size and type of the numbers in the axis? how adjust them?

## @tomleslie yes i am running pa...

@tomleslie yes i am running parallel computations, and monitoring cpu and ram usage during the run, min 44% of ram and 11% of Cpu are free. but error persists.

## @acer I modified my ws as below:  ...

@acer I changed GAMMA to PI and modified my ws as below:

 >
 >
 (1)
 >
 (2)
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 (3)

I created 2 UPs procs (UP1,UP2), the only difference between them is in using Rr as argument.

First of all I thick y1 is computed faster than y2. Why?

Second, why y1-y2 is not zero???

I am amazed.

## Thanks all of you....

@Carl Love Thanks acer, vv and Carl Love. It was great help.

 1 2 3 4 5 6 7 Page 3 of 9
﻿