How do I solve for eigenvalues symbolically?

    adoty = [Na      0      0    -1    ];
    adotp = [0       Na     1     0    ];
    qdot  = [Mpa*wr  Ma     Mq   -I*wr ];
    rdot  = [-Ma    -Mpa*wr I*wr -Mq   ];

    A    = [adoty; adotp; qdot; rdot];
    EigA = eig(A)

So, for this matrix A, I can find the eigenvalues using this function. However, I would like to be able to get the eigenvalues as a function of wr. I can make a for loop to post a new eigenvalue per unit wr, but I don't have the ability to find the correlation between all those outputs since they're outputted as 4 different values. I tried to use eig(sym(A)), but the results end up being about 5 pages per answer, so that is useless to me. How do I find the correlation between my eigenvalues as a function of that wr value?

Attached is the program with all the numerical values if that is of any help.

Answer 1

[madhan ravi]

syms Na Ma Mq I wr Mpa
    adoty = [Na      0      0    -1    ];
    adotp = [0       Na     1     0    ];
    qdot  = [Mpa*wr  Ma     Mq   -I*wr ];
    rdot  = [-Ma    -Mpa*wr I*wr -Mq   ];

      A    = sym([adoty; adotp; qdot; rdot]);
      Eig(wr)=eig(A)

like this you mean?

[Albert Garcia]

Absolutely! Thank you so much mate! I was trying to turn those variables into strings but that didn't work out. Didn't know I could just make stuff symbols to begin with. Great help

[madhan ravi]

you're welcome :)