As far as I know, ‘bc2’, Maxima’s built-in function for solving two-point boundary value problem only handles the type:
For example, solving
But ‘bc2’ can not be applied to
since it is not the type of (*). However, you can roll your own:
An error occurs on the line where the second boundary condition is specified. It attempts to differentiate the solution with respect to 
diff(rhs(sol), 0);
which is absurd.
The correct way is to express the boundary conditions using ‘at’ instead of ‘ev’:
The following works too as the derivative is obtained before using ‘ev’:
Nonetheless, I still think using ’at’ is a better way:
