Hi everybody,

i have two problems by using Sage to solve a technical difficulty and hope you can help me with it.

I want to calculate a bending line of a girder with two regions with different second moments of areas.

The first Problem is About the handling of sage. I get a solution, but i can't extract the single variables of the solution list.

The second Problem is that the solution isn't correct. I calculated the bending line at Excel and the area of the greater diameter has a wrong position and gradient. I don't know where the error is.

I hope you can help me and please excuse my bad english talent!

Greatings

Forceman

```
reset
var ('a_s,a,b,c,M,F,Delta,x_I,x_II')
var ('E,I_R,I_W')
var ('C_1_I,C_1_II,C_1_III,C_1_IV,C_2_I,C_2_II,C_2_III,C_2_IV')
MI=function('MI')(x_I)
MII=function('MII')(x_II)
MI = F*(a_s-x_I)
MII = M-M/c*x_II
w_sI = MI.integral(x_I)+C_1_I
w_I = (w_sI.integral(x_I)+C_2_I)/(E*I_W)
w_sI = w_sI/(E*I_W)
w_sII = MII.integral(x_II)+C_1_II
w_II = (w_sII.integral(x_II)+C_2_II)/(E*I_W)
w_sII = w_sII/(E*I_W)
w_sIII = MII.integral(x_II)+C_1_III
w_III = (w_sIII.integral(x_II)+C_2_III)/(E*I_R)
w_sIII = w_sIII/(E*I_R)
w_sIV = MII.integral(x_II)+C_1_IV
w_IV = (w_sIV.integral(x_II)+C_2_IV)/(E*I_W)
w_sIV = w_sIV/(E*I_W)
sols = solve([
w_I(x_I=0)==0,
w_sI(x_I=0)==w_sII(x_II=0),
w_II(x_II=0)==0,
w_II(x_II=a)==w_III(x_II=a),
w_sII(x_II=a)==w_sIII(x_II=a),
w_III(x_II=b)==w_IV(x_II=b),
w_sIII(x_II=b)==w_sIV(x_II=b),
w_IV(x_II=c)==0],
C_2_IV,C_2_III,C_2_II,C_2_I,C_1_IV,C_1_III,C_1_II,C_1_I,solution_dict=true)
C_1_I = C_1_I.subs(sols[0])
C_1_II = C_1_II.subs(sols[0])
C_1_III = C_1_III.subs(sols[0])
C_1_IV = C_1_IV.subs(sols[0])
C_2_I = C_2_I.subs(sols[0])
C_2_II = C_2_II.subs(sols[0])
C_2_III = C_2_III.subs(sols[0])
C_2_IV = C_2_IV.subs(sols[0])
```