function [intlifttwist,intliftbend,intmassandother]=integralbend(x,j,v)
W=567000;%v=245.565;
s=35.35/2;
lamb=35*(pi/180);
c=3.75;
seff=s/cos(lamb);
%M=0.8;altitude =8230
cmac=-0.01;
ei=2.23*10^8;
gj=1.431*10^8;
m=575;
e1=0.25;
e2=0.1;
a=5.42468;
pho=0.5112;
y1=(0:0.001:x);
y2=(x:0.001:1);
lifttwistload1=(1/ei)*((y1.^2)/2).*(((-j.*exp(-j)*y1)-(1.*(exp(-j*y1)))+1))*a*(1/c)*cos(lamb)*cos(lamb);
liftbendload1=-(1/ei)*((y1.^2)/2).*((-j*exp(-j)*y1)-(exp(-j*y1))+(0.5*j*j*exp(-j)*((y1-1).^2))-(j*j*0.5*exp(-j))+1)*a*(1/c)*cos(lamb)*sin(lamb);
massload1=-(1/ei)*((y1.^2)/2)*(4*m*9.81*(1/c)*cos(lamb));
lifttwistload2=(1/ei)*(((x*y2)-(0.5*x*x)).*(((-j*exp(-j)*y2)-(exp(-j*y2))+1)))*a*(1/c)*cos(lamb)*cos(lamb);
liftbendload2=-(1/ei)*(((x*y2)-(0.5*x*x)).*((-j*exp(-j)*y2)-(exp(-j*y2))+(0.5*j*j*exp(-j)*((y2-1).^2))-(j*j*0.5*exp(-j))+1))*a*(1/c)*cos(lamb)*sin(lamb);
massload2=-(1/ei)*((x*y2)-(0.5*x*x)*(4*m*9.81*(1/c)*cos(lamb)));
lifttwistmom2=(1/ei)*((sqrt(y2.*y2)-x).*((-j*exp(-j)*y2)-(exp(-j*y2))+1))*e1*a*cos(lamb)*cos(lamb)*sin(lamb);
liftbendmom2=-(1/ei)*((sqrt(y2.*y2)-x).*((-j*exp(-j)*y2)-(exp(-j*y2))+(0.5*j*j*exp(-j)*((y2-1).^2))-(j*j*0.5*exp(-j))+1))*e1*a*cos(lamb)*sin(lamb)*sin(lamb);
massandothermom2=(1/ei)*(sqrt(y2.*y2)-x).*(((cmac*cos(lamb)*sin(lamb)))-((1/(0.5*pho*v*v*c*c*seff))*4*m*9.81*c*e2*seff*cos(lamb)*sin(lamb)));
intlifttwist=trapz(y1,lifttwistload1)+trapz(y2,lifttwistload2)+trapz(y2,lifttwistmom2);
intliftbend=trapz(y1,liftbendload1)+trapz(y2,liftbendload2)+trapz(y2,liftbendmom2);
intmassandother=trapz(y1,massload1)+trapz(y2,massload2)+(((0.5*pho*v*v*c*c*seff))*trapz(y2,massandothermom2));
end