//Jason Heath
//Scilab program to integrate a 4 state variable mass-spring-damper system

// Define Constants in System
F = 1;
M1 = 1;
M2 = 1;
Kd1 = 1;
Kd2 = 1;
Kd3 = 1;
Ks1 = 1;
Ks2 = 1;
Ks3 = 1;
x1 = 0;
x2 = 0;
v1 = 0;
v2 = 0;


// Define time step

h = .1;
steps = 100;

for i = 0: (steps-1),

dx1 = v1;
dv1 = v1 * ((-Kd1 - Kd2)/M1) + x1 * ((-Ks1 - Ks2)/M1) + v2 * (Kd2/M1) + x2 * (Ks2/M1);

dx2 = v2;
dv2 = v2 * ((-Kd2 - Kd3)/M2) + x2 * ((-Ks2 - Ks3)/M2) + v1 * (Kd2/M2) + x1 * (Ks2/M2) + (F/M2);

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

end

printf("The values of the differential equations after 10 seconds are:\n\n");
printf("%f, %f, %f, %f", x1, v1, x2, v2);