syms x x1 x2 x3 x4 x5 x6 F3w F4w F5w F3p F4p F5p
assume([x x1 x2 x3 x4 x5 x6 F3w F4w F5w F3p F4p F5p],'real')
L2=60000;
L1=40000;
%SECOND TRUSS
C1=[1 x3 x3.^2; 1 x4 x4.^2; 1 x5 x5.^2];
C1_=inv(C1);
N1=[1 x x^2]*C1_;
u1=(simplify(N1));
%pretty(u1);
B1=diff(u1,x);
%pretty(B1);
E=30E3;
A=pi*(3500^2-3000^2);
K1=B1'*A*E*B1*L2;
%pretty(K1);
x3=40000;
x4=70000;
x5=100000;
K1p=eval(K1);
%fun = @(x) K1p;
fun = matlabFunction(K1p);
K1pp=integral(fun,40000,100000,'ArrayValued',1)