How to calculate the forces in a statically indeterminate beam.

%calculation of forces in a ststically indeterminate beam.

clear;

q=input('Enter the loads, n/m');

a=input('Enter the span lenght of the loads, m');

l=input('Enter the length of the beam, m');

b=input('Enter the width of the beam, m');

kk=input('Enter step for graphics');

k=a/l;

st=l/kk;

x=0:st:1;

for i=1:kk+1

Ma(i)=(q*x(i)^2*k^2*(3-(4*k)+(1.5*k^2)))/6;

Mb(i)=(q*x(i)^2*k^2*(k-(0.75*k^2)))/3;

A(i)=((Ma(i)-Mb(i))/x(i))+(q*x(i)*k*(1-(k/2)));

B(i)=((Mb(i)-Ma(i))/x(i))+((q*x(i)*k^2)/2)-(q*a)-A(i);

end

plot(x,Ma,'r',x,Mb,'bl')

Please can any one help I want to sove this problem. I want to know how to get the results as graph as well as put them in tables.

condition 0.1<=k<=1