solvepde(model)

Здравствуйте!
Помогите пожалуйста решить проблему!
При решении системы дифференциальных уравнений в частных производных

model = createpde(3);

syms u1(x,y) u2(x,y) u3(x,y) ;
pdeeq = [2*(r_radius*r_radius-x*x)*diff(u1,x,y)-x*diff(u1,y)-x*sqrt((r_radius*r_radius-x*x))*diff(u2,x,y)+x*x*diff(u2,y,y)-sqrt((r_radius*r_radius-x*x))*diff(u2,y)+4*(r_radius*r_radius-x*x)*diff(u3,x,x)+2*r_radius*r_radius*(1+nu_coef)*diff(u3,y,y)-2*x*diff(u3,x) + 2*prEh_coef*(1+nu_coef)/(nu_coef*len); sqrt((r_radius*r_radius-x*x))*diff(u1,y,y)+sqrt((r_radius*r_radius-x*x))*diff(u2,x,y)+2*sqrt((r_radius*r_radius-x*x))*diff(u3,x,y); 2*(r_radius*r_radius-x*x)/r_radius*diff(u1,x)-2*x*sqrt((r_radius*r_radius-x*x))/r_radius*diff(u2,x)+r_radius*nu_coef*diff(u3,y) - (R_radius-y)/len*prEh_coef-(3-nu_coef*nu_coef)*prEh_coef];
symCoeffs = pdeCoefficients(pdeeq,[u1 u2 u3],'Symbolic',true);

coeffs = pdeCoefficients(pdeeq,[u1 u2 u3]);

specifyCoefficients(model,'m',coeffs.m,'d',coeffs.d, ...
'c',coeffs.c,'a',coeffs.a,'f',@fcoeffunction);

rmix = [0 0 0];
qmix = [0 0 0; 0 0 0; 0 0 0];
gmix = [0 0 -nu_coef*prEh_coef_2];
applyBoundaryCondition(model,'neumann','Edge',1,'g',[0 prEh_coef_2*(1-nu_coef)/nu_coef 0],'q',[0 0 0; 0 0 0; 0 0 0]);
applyBoundaryCondition(model,'dirichlet','Edge',2,'r',[prEh_coef 0 0],'h',[1 0 0; 0 1 0; 0 0 0]);
applyBoundaryCondition(model,'mixed','Edge',3,'h',@myufun,'r',rmix,'g',gmix,'q',qmix);
applyBoundaryCondition(model,'dirichlet','Edge',4,'r',[-prEh_coef 0 0],'h',[1 0 0; 0 1 0; 0 0 0]);

res = solvepde(model);

возникает следующая ошибка:
Error using formGlobalMatrix
Coefficient evaluation function "makeCoefficient/coefficientFunction" returned a complex matrix but previously returned a real matrix. This function must consistently return the same type of matrix.(Пользуюсь пробной версией R2021b)
Подскажите, пожалуйста, что может быть причиной данной ошибки. Заранее благодарю!