Do you see that as being substantially different from having a single worksheet with all the code in one execution group?

@tuGUTS Yes, you can add initial and boundary conditions. These will, if they work at all, remove some or all of the arbitrary functions and constants from the solution. I also don't know what would be appropriate initial and boundary conditions for this PDE. But see ?pdsolve for numerous examples.

@tuGUTS Yes, you can add initial and boundary conditions. These will, if they work at all, remove some or all of the arbitrary functions and constants from the solution. I also don't know what would be appropriate initial and boundary conditions for this PDE. But see ?pdsolve for numerous examples.

You equations are effectively a system of ordinary differential equation. Try changing alpha(r,t) to alpha(r) and then using dsolve. (Although I wouldn't be surprised if pdsolve had already tried that.) It might help to make assumptions on c and kappa such as assume(c>0, kappa>0).

I would try these things myself, but I don't feel like typing in your equations. Please post them in plaintext or upload a worksheet.

Are you saying that you want to extract the point data from the plot?

What is Veq?

Can you post the equations?

Are you getting an error message from pdsolve, or no result at all?

A simplified version of your plotting command:

plot(
     [zip(`[]`,X1, Y1), zip(`[]`, X2, Y2)],
     color= [blue, red],
     style= point, symbol= solidcircle, symbolsize= 8,
     view= [0 .. 2, 0 .. 1], scaling= constrained
);

A simplified version of your plotting command:

plot(
     [zip(`[]`,X1, Y1), zip(`[]`, X2, Y2)],
     color= [blue, red],
     style= point, symbol= solidcircle, symbolsize= 8,
     view= [0 .. 2, 0 .. 1], scaling= constrained
);

It is not clear what you mean by "solve", "solving for", "equation", and "plotting". Do you mean

to solve the equations x/(x+1) = .5, x/(x+1) = .9, etc., for x?

Or do you mean

to evaluate the expression x/(x+1) at x = .5, .9, etc.?

And, in either of the above situations, it is not clear what "plotting" means. We can plot expressions (but not equations) with one variable. We can plot equations with two variables. We can plot a sequence (not a series) by using the index number into the sequence as one of the coordinates.

@andy.zhou.nuaa I recommend you download the free DirectSearch package from the Maple Applications Center. It contains a solver similar to fsolve, but with more flexibility in the precision.

@Markiyan Hirnyk He means the Answer by Alec Milhailovs to which you referred.

@Markiyan Hirnyk He means the Answer by Alec Milhailovs to which you referred.

@Sergio Parreiras Yes, I see what you mean. You can use simplify (non symbolic) after simplify(..., symbolic). I'm sure someone would take a closer look at how to simplify the expression if you post the code for the expression.

@Sergio Parreiras Yes, I see what you mean. You can use simplify (non symbolic) after simplify(..., symbolic). I'm sure someone would take a closer look at how to simplify the expression if you post the code for the expression.