集成和求解线性系统python中的竞争者

Question

我试图在此处集成2个功能（结果为绿色）

并根据初始值（a，b，c）求解常数C1 C2，以获得完整的方程式。到目前为止，这是我的脚本：

#Defining the symbols for sympy
x=sp.symbols('x')
C1=sp.symbols('C1')
C2=sp.symbols('C2')
#input
el = 7 #length
q=10*x #force
#Condition a) from the Picture M(0)=0, M(el)=0
cond1=(1,0) #condition one Q(0)!=0, M(0)=0
cond2=(1,0) #condition one Q(el)!=0, M(el)=0
#Integration
Q= sp.integrate(q,x)+C1
M= sp.integrate(Q,x)

现在我必须检查cond1和cond2中的元素是否为= 0。

if cond1[0]==0 :
    print("Q(0)=0")
    Q0=Q.sub({x:0}

if cond1[1]==0 :
    print("M(0)=0")
    M0.sub({x:0}

...

这里有更快的方法吗？ ...

现在我得到了Q0 M0 QL ML

求解我的线性系统A ^（-1）* b = [C1，C2]，其中b始终为0 b = np.array（[[0]，[0]]） 但是根据公式构建A矩阵会让我头疼。我怎样才能做到这一点？