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,
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)