MaplePrimes - answers and comments on Question, Plotting Real and Imaginary Parts?

sorry I meant to include:

ndattani — Fri, 22 Aug 2008 05:11:25 Z

sorry I meant to include: deltaE = 0.355

try

acer — Fri, 22 Aug 2008 06:48:40 Z

plot( Re@subs(deltaE=0.355,M=sqrt(5.4e5),eval(Eig1)), 1..10);

plot( Re(subs(deltaE=0.355,M=sqrt(5.4e5),Eig1(R))), R=1..10);

plot( R->evalf[100](Re(subs(deltaE=0.355,M=sqrt(5.4e5),Eig1(R)))), 1000..1010);

And for fun,

> limit( Re(subs(deltaE=0.355,M=sqrt(5.4e5),eval(Eig1))(R)), R=infinity);
                                 0.3550000000

However,

> limit( Re(eval(Eig1)(R)), R=infinity);
                                                          2
       (1/3)           3 (1/3)           (2/3)      deltaE
-1/9 27      Re((deltaE )     ) - 1/27 27      Re(------------)
                                                         3 1/3
                                                  (deltaE )
 
     + 2/3 Re(deltaE)
 
> simplify(%) assuming deltaE>0;
                                       0

acer

Thanks a lot! I used the

ndattani — Mon, 25 Aug 2008 02:09:12 Z

Thanks a lot!

I used the second command,

plot( Re(subs(deltaE=0.355,M=sqrt(5.4e5),Eig1(R))), R=1..10);

and it worked!

In your last line, you evaluated that limit to be 0 as R->infinity, although the plot shows that it goes to 0.355 (which is deltaE). And this is what I expected from theory as well.

I'm curious now as to how I could plot the same function against R^3 instead of vs R, without making a new variable equal to R^3 and redefining the function in terms of the new variable then plotting. Any ideas???

I mean, I can do plot(

ndattani — Mon, 25 Aug 2008 02:22:23 Z

I mean, I can do

plot( Re(subs(deltaE=0.355,M=sqrt(5.4e5),Eig1(R^3))), R=1..10);

but ideally what I would like to be able to do is

plot( Re(subs(deltaE=0.355,M=sqrt(5.4e5),Eig1(R^3))), R^3=1^(1/3)..10^(1/3));

but Maple won't let me do that... is there an alternative?

That doesn't work, When I

ndattani — Mon, 25 Aug 2008 20:26:44 Z

That doesn't work,

When I do

> plot( Re(subs(deltaE=0.355,M=sqrt(5.4e5),Eig1(R))),
>       R=1..10^(1/3)

Maple returns the correct plot,

When I change " Eig1(R) to Eig1(t) " , without changing "R =1..10^(1/3)", I get the same result.

But when I change "R =1..10^(1/3)" to "t=1..10^(1/3)" I get an error message.


Also, I'm not sure if that would give me the output I want,
because if I plot x^2 vs x I should get a parabola ,
and if I plot x^2 vs x^2 I should get a straight line,
so the shape actually changes as well, and it's not just a matter of labeling.

Is there a way I can plot a function of x vs x^2 instead of vs x?

Actually

ndattani — Mon, 25 Aug 2008 20:32:41 Z

Actually, it would be even better if I could plot a function of x vs 9x^2