@Markiyan Hirnyk didn't you get the error?

@Markiyan Hirnyk didn't you get the error?

I'm always getting the error:

Error, (in gfun:-diffeqtorec) invalid input: getname expects its 1st argument, yofz, to be of type function(name), but received -(3/2)/x(n)+(3/2)/(x(n)-1)+(3/2)/(x(n)-z(n))

the DE is:

ode:=x^2*(x-1)*(x-z)*(diff(X(x), x, x, x)) = (2*x^2-x*lambda-z)*(diff(X(x), x))+(-(3/2)*x^3+(3/2)*x*z)*(diff(X(x), x, x))-X(x)*x+X(x)*lambda

I'm always getting the error:

Error, (in gfun:-diffeqtorec) invalid input: getname expects its 1st argument, yofz, to be of type function(name), but received -(3/2)/x(n)+(3/2)/(x(n)-1)+(3/2)/(x(n)-z(n))

the DE is:

ode:=x^2*(x-1)*(x-z)*(diff(X(x), x, x, x)) = (2*x^2-x*lambda-z)*(diff(X(x), x))+(-(3/2)*x^3+(3/2)*x*z)*(diff(X(x), x, x))-X(x)*x+X(x)*lambda

Thanks first. Yeah it indeed took very long before stating the error.

I used RootOf now which does give a result relatively quick.

0.3033380032 seems to be correct.

But regarding the other way Preben chose it's definitely something to consider in the future if similar problems occur.

Thanks first. Yeah it indeed took very long before stating the error.

I used RootOf now which does give a result relatively quick.

0.3033380032 seems to be correct.

But regarding the other way Preben chose it's definitely something to consider in the future if similar problems occur.

Hey, Thanks again.

Now for better understanding:

Do you perhaps have a typical example for when to use radnormal over evalc or vice versa?

and specifically why radnormal/evalc does not work in either case?

Hey, Thanks again.

Now for better understanding:

Do you perhaps have a typical example for when to use radnormal over evalc or vice versa?

and specifically why radnormal/evalc does not work in either case?

Thanks.

But why is evalc your friend here?

Other times the solution for simplification of an expression was radnormal.

I'm still not sure when to use what.

Thanks.

But why is evalc your friend here?

Other times the solution for simplification of an expression was radnormal.

I'm still not sure when to use what.

@Carl Love 

Yeah but in this case I parametrically plot y against x (for all t from the interval), whilst in each odeplot I plot x against t and y against t so using display I get 2 graphs in one plot but not y against x.

@Carl Love 

Yeah but in this case I parametrically plot y against x (for all t from the interval), whilst in each odeplot I plot x against t and y against t so using display I get 2 graphs in one plot but not y against x.

Thank you very much ;)

edit: well infact I meant:

Sol:= dsolve({ode3, ics3, ode, ics}, numeric):
plots:-odeplot(Sol, [x(t),y(t)], t= 0..3);

does this also work with display command?

Thank you very much ;)

edit: well infact I meant:

Sol:= dsolve({ode3, ics3, ode, ics}, numeric):
plots:-odeplot(Sol, [x(t),y(t)], t= 0..3);

does this also work with display command?

sorry for a third question:

But how can I do this:

Let xs be the numerical solutions as above:

xs:= dsolve({ics, ode}, numeric);
plots:-odeplot(xs, t= 0..3);

works without any problems

how can one plot:

plots:-odeplot(f(xs), t= 0..3);

with some function f

it says:

Error, (in plots/odeplot) input is not a valid dsolve/numeric solution