Hi everyone

I have solved an ODE in Maple with the following command:

soln := dsolve({ode1, f(0) = 0, f(1) = 1, (D(f))(0) = 0, (D(f))(1) = 0}, f(eta), numeric)

where ode1 is:

ode1 := diff(f(eta), [`$`(eta, 4)])+R*(eta*(diff(f(eta), [`$`(eta, 3)]))+3*(diff(f(eta), [`$`(eta, 2)]))-f(eta)*(diff(f(eta), [`$`(eta, 3)]))+(diff(f(eta), eta))*(diff(f(eta), [`$`(eta, 2)]))) = 0

then I saved the soln plot in p1 as:

p1 := odeplot(soln, [eta, f(eta)...