### desolve_system where boundary values are functions of the dependent variable

I am using `desolve_system()`

to solve a boundary value problem. The equations are very simple but I'm having a bit of difficulty trying to include the boundary conditions. The system of equations and boundary conditions are,

http://i.imgur.com/kpY3N.png

This is the code that I'm using in sage. How can I include the boundary values in this problem?

```
reset()
var('I_1 I_2 I0 R0 R1 x alpha D')
I_1 = function('I_1', x)
I_2 = function('I_2', x)
eq1 = diff(I_1,x) == -alpha * I_1
eq2 = -diff(I_2,x) == -alpha * I_2
sol = desolve_system( [eq1, eq2], [I_1, I_2], ivar=x)
view(sol)
# How can I include these?
#BV1: I_1(x=0)==A*I_2(x=0) + C
#BV2: I_2(x=D)==B*I_1(x=D)
```