Here is my code:
#include int main() { float dt = 0.0002, t0; int it0; for ( t0 = 0.0512; t0 <= 0.0544; t0 += 0.0004