以c++为基础编写一个弹簧质点模型

好的，以下是一个简单的弹簧质点模型的 C 代码：

#include <stdio.h>
#include <math.h>

#define GRAVITY 9.81
#define TIME_STEP 0.01
#define MASS 1.0
#define SPRING_CONSTANT 10.0
#define REST_LENGTH 1.0

int main() {
double position = 0.0;
double velocity = 0.0;
double acceleration = 0.0;
double force = 0.0;
double time = 0.0;

while (1) {
// Calculate force
force = -SPRING_CONSTANT * (position - REST_LENGTH);

// Calculate acceleration
acceleration = force / MASS;

// Update velocity and position
velocity += acceleration * TIME_STEP;
position += velocity * TIME_STEP;

// Print results
printf("Time: %f, Position: %f, Velocity: %f\n", time, position, velocity);

// Check if the spring has reached equilibrium
if (fabs(position - REST_LENGTH) < 0.001 && fabs(velocity) < 0.001) {
printf("Equilibrium reached!\n");
break;
}

// Update time
time += TIME_STEP;
}

return 0;
}