On Sat, Mar 12, 2016 at 5:27 AM, Doug Stewart <address@hidden> wrote:
yes with your numbers it is now at 5 min. and still going.try putting in a small number for your zeros
that did not change the time!!!!! wow this is the ode45 doing its work.
You could set the tolerances larger
--On Sat, Mar 12, 2016 at 4:18 AM, Bear <address@hidden> wrote:m_p = 0.003; %mass of pellet ( kg )
m_s = 0.050; %mass of spring piston + 1/3 spring mass ( kg )
r_piston = 0.023; %radius of spring piston in air gun ( m )
c = 1.7272; %caliber ( m )
R = 0; %twist rate of rifling ( 1/m )
f_e = 0; %experimentally quantified elastic compression friction force for
pellet-style ( N )
mu = 0.3; %kinetic coeficient of friction between pellet and barrel
materials
I_p = 0.00223741488; %mass moment of inertia of pellet along its
lognitudinal central axis ( kg*m^2 )
g = 9.801; %force of local gravitational feild (use standard earth value for
analysis) ( m/(kg*s^2) )
P_o = 1.013*10^5; %local atmospheric pressure (normally set at see level for
analysis) ( Pa )
V_o = 0.000006554826; %Volume of air in cylinder when the air gun is
"cocked" ( m^3 )
k = 3000; %spring constant ( N/m )
x_s_f = 0.3226; %compression of spring assembly after cocking ( m )
x_s_o = 1; %compression of spring assembly before cocking ( m )
barrel_length = 1; % ( m )
Windows 8.1 64bit, latest version (4.0.0)
--
View this message in context: http://octave.1599824.n4.nabble.com/How-Come-My-Program-Won-t-Work-tp4675378p4675431.html
Sent from the Octave - General mailing list archive at Nabble.com.
_______________________________________________
Help-octave mailing list
address@hidden
https://lists.gnu.org/mailman/listinfo/help-octave