//
// Chris Hough
// Final SciLab Program for Mechanical System
// 

m1 = 1;
m2 = 1;
ks1 = 1;
kd1 = 1;
ks2 = 1;
kd2 = 1;
ks3 = 1;
kd3 = 1;
f = 1;

x1 = 0;         
v1 = 0;
x2 = 0;
v2 = 0;

h = 0.1;
steps = 100;


// Loop for integration



for i=1:(steps-1),

dx1 = v1;
dv1 = ((-kd1-kd2)/m1)*v1+((-ks1-ks2)/m1)*x1+(kd2/m1)*v2+(ks2/m1)*x2;
dx2 = v2;
dv2 = ((-kd2-kd3)/m2)*v2+((-ks2-ks3)/m2)*x2+(kd2/m2)*v1+(ks2/m2)*x1+(f/m2);

x1 = x1 + h * dx1;
v1 = v1 + h * dv1;
x2 = x2 + h * dx2;
v2 = v2 + h * dv2;
end

printf("%f, %f, %f, %f", x1, v1, x2, v2); // This will print out the values at t=10.