function [vq1] = neilijisuanvq1(x,x1,x2,l1,E,I)
vq1=zeros(size(x));
if x>=0&x<=x1
vq1=1/(6*E*I*l1)*((-0.5*x1^2+x1*x2-0.5*x2^2)*x^3+(0.5*x2^3-1.5*x1^2*x2+x1^3)*x^2);
elseif x>x1&x<=x2
vq1= 1/(6*E*I*l1)*(-0.05*x^5+0.25*x2*x^4-0.5*x2^2*x^3+0.5*x2^3*x^2+(0.75*x1^4-x1^3*x2)*x+0.25*x1^4*x2-0.2*x1^5);
else
vq1=1/(6*E*I*l1)*((0.25*x2^4-x1^3*x2+0.75*x1^4)*x-0.2*x1^5+0.25*x1^4*x2-0.05*x2^5);
end