Trouble computing a 4th order piecewise ODE (3 parts) symbolically

Hello I have tried to run the following code, but after hours, the program has not found a solution and comes up with an error. what is wrong with my code

Thanks

%*************************************************************************

% Code

%*************************************************************************

%Variables

k=20e6; % N/m^3

Gp=10e6; % N/m

gamma=22e3; % N/m^3

We=6.54; % m

he=6.655; % m

phi=30; % deg

coh=150e3; % Pa K0=1; %

H=10; % m

q=0; % Pa ?

delta=0.4; % m

a=1.5; % m

E=200e9; % Pa

I=10e12; % Pa

b=0.4; % m

q0=1; % ?????????

D=0.102; % m

t=0.01; % m

Lf=30; % m

w0=0.014; % m

%Background Equations

B=We/2;

phirad=phi*3.14/180;

B1=B+he*tan(3.14/4+phirad/2);

d=he*tan(3.14/4+phirad/2);

bstar=b*(1+((Gp/k)^(1/2))/b);

%Geometry Equations

PointA=0;

PointB=1.5*a;

PointC=PointB+d;

PointD=Lf;

%Stress Equations

P=(gamma*B1-coh)/(K0*tan(phirad))*(1-exp(-((K0*H*tan(phirad))/B1)))+q*exp(-((K0*H*tan(phirad))/B1)); % Load acting on the element

%Symbols

syms w5(x) w6(x) w7(x) w(x)

dw5 = diff(w5);

d2w5 = diff(w5,2);

d3w5 = diff(w5,3); d4w5 = diff(w5,4);

dw6 = diff(w6);

d2w6 = diff(w6,2);

d3w6 = diff(w6,3);

d4w6 = diff(w6,4);

dw7 = diff(w7);

d2w7 = diff(w7,2);

d3w7 = diff(w7,3);

d4w7 = diff(w7,4);

[w5(x), w6(x), w7(x)] = ...

    dsolve(...

    diff(w5,4) == b*q0/(E*I), ...                                           %fourth order DE Section AB

    diff(w6,4) == Gp*bstar*d2w6/(E*I)-k*bstar*w6/(E*I)+b*q0/(E*I), ...      %fourth order DE Section BC

    diff(w7,4) == Gp*bstar*d2w7/(E*I)-k*bstar*w7/(E*I), ...                 %fourth order DE Section CD

    w5(PointA) == w0, ...                   %B1 Known displacement at start of element

    dw5(PointA) == 0, ...                   %B2 No rotation due to fixed reaction at start of element 

    w5(PointB)== w6(PointB), ...            %B3 CE between section AB-BC displacement

    dw5(PointB)== dw6(PointB), ...          %B4 CE between section AB-BC rotation

    d2w5(PointB)== d2w6(PointB), ...        %B5 CE between section AB-BC moment

    d3w5(PointB)== d3w6(PointB), ...        %B6 CE between section AB-BC shear

    w6(PointC)== w7(PointC), ...            %B7 CE between section BC-CD displacement

    dw6(PointC)== dw7(PointC), ...          %B8 CE between section BC-CD rotation

    d2w6(PointC)== d2w7(PointC), ...        %B9 CE between section BC-CD Moment

    d3w6(PointC)== d3w7(PointC), ...        %B10 CE between section BC-CD Shear

    d2w7(PointD)== 0, ...                   %B11 assumed no moment at end of element

    d3w7(PointD)== 0)                       %B12 assumed no shear at end of element

% B1-B12 Boundary conditions

% CE Continunity Equations

*************this this the error that comes from it***********

Error using mupadengine/feval (line 157) MuPAD error: Error: Wrong type of first argument. [output::tableForm]

Error in dsolve (line 206) dropped = char(eng.feval('output::tableForm',dropped,'String'));

Error in HelpNew (line 59) [w5(x), w6(x), w7(x)] = ...