I am also trying to do somthing similar, I have learned how to solve differential equations in class, I am just trying to set up the differential equation, (I took calc but trig bassed physics) I have the air density formula, the force of drag formula, and basic contants (Drag coeficiant, frontal area, initial velocity, mass of projectile etc..) I need with the projectile I am using. So far I know that the velocity of time would be somthing like this V(t) = Vi - (Force of drag/mass of projectile)*t. Unless I have it wrong please correct me. No the problem is that Force of drag (Fd) requires the V(t) to calculate Fd. Ie V(t) = Vi - (905.8656 - ((.5 * p * (V(t)^2) * (4.56036731 × 10-5 ) *.507) / .0116638038)t). So I know I need calc here its just setting it up is getting me. Anyways I think this is my final formula sombody please check it I just tried solving using substitution but then again I think I nuked it, 912952.4078072/(p*time) = V(t)
Anyways one other side note, I used Vi as a constant but I will put as inputed variable into the program. I am debating on writing it as Ti-Basic (for 83/84 most likley) or as a .exe. In anycase I will post it when I am done onto the forums. If i do .exe I will have it input any mass, frontal area, and drag coeficient as a variable, also tied into the program will be manual entering of forcasted air pressure, humidity, and temperature, to further calculate the velocity of time t.
contants
p= density of fluid (air) for which I have a formula it is in kg/m^3 (I can take care of this for the most part, I will just leave it as a variable, I have a separate program all set for that bassed on temp and measured pressure and humidity. (takes altitude out of the equation)
V(t)= Velocity of time
Vi = Initial velocity (will be as variable in the program) = 905.8656 m/s
Mass of projectile if needed = 11.6638grams
Side note: I realized not everything has correct number of significant digits, I am going by conversions and ctrl + C...Ctrl + V for numbers, it should get me close enough to where the program will round to the nearest m/s anyhow.