E30106 I500 A10 La1012 Lbsqrt826212 Lc1612 na0 1012La nb812

E=30*10^6;
I=500;
A=10;
La=10*12;
Lb=sqrt(8^2+6^2)*12;
Lc=16*12;

na=[0 10*12]./La;
nb=[8*12 16*12]./Lb;
nc=[8*12 0]./Lc;

pf=[10;15;15;0;0;0];

ka=frameElement2D(E,A,I,La,na);
kb=frameElement2D(E,A,I,Lb,nb);
kc=frameElement2D(E,A,I,Lc,nc);

K=zeros(12,12);
K(1:6,1:6)=K(1:6,1:6)+ka;
K(5:10,5:10)=K(5:10,5:10)+kb;
K(7:12,7:12)=K(7:12,7:12)+kc;

kff=K(6:11,6:11);

uf=kff^-1*pf;

  
when I run this code I get a blank command window, what is wrong?

Solve Problem by using a computer program For the rigid frame shown in Figure P5-19, determine (1) the nodal displacement components and (2) the support reactions. (3) Draw the shear force and bending moment diagrams. For all elements, let E 30 x 106 psi, I-200 in.4, and A 10 in 15 k-ft 15 kip 6 f 10 kip 16 ft 10 ft 8 ft Figure P5-19

Solution

firstly correct the value of I=200;

correct the values also:

na=[10*12]./La;
nb=[8*12 +16*12]./Lb;
nc=[8*12 ]./Lc;

I hope you can find the solution after correcting this.