Changeset 198
 Timestamp:
 Mar 10, 2007, 2:08:45 PM (16 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

branches/newlanguage/sample/miscellaneous/sample_car.mso
r176 r198 49 49 50 50 EQUATIONS 51 diff([xl, yl, xr, yr]) = q;51 diff([xl, yl, xr, yr])*'s' = q; 52 52 53 eps^2*M/2 * diff(q(1)) = (l0ll)*xl/ll53 eps^2*M/2 * diff(q(1))*'s' = (l0ll)*xl/ll 54 54 +lambda(1)*xb +2*lambda(1)*(xlxr); 55 55 56 eps^2*M/2 * diff(q(2)) = (l0ll)*yl/ll56 eps^2*M/2 * diff(q(2))*'s' = (l0ll)*yl/ll 57 57 +lambda(1)*yb +2*lambda(2)*(yryl)  eps^2*M/2; 58 58 59 59 60 eps^2*M/2 * diff(q(3)) = (l0lr)*(xrxb)/lr60 eps^2*M/2 * diff(q(3))*'s' = (l0lr)*(xrxb)/lr 61 61 2*lambda(2)*(xlxr); 62 62 63 eps^2*M/2 * diff(q(4)) = (l0lr)*(yryb)/lr63 eps^2*M/2 * diff(q(4))*'s' = (l0lr)*(yryb)/lr 64 64 2*lambda(2)*(ylyr)  eps^2*M/2; 65 65 … … 70 70 xb = sqrt(l^2yb^2); 71 71 72 yb = r*sin(w*time );72 yb = r*sin(w*time*'rad/s'); 73 73 74 74 ll = sqrt(xl^2 + yl^2); … … 86 86 TimeStep = 0.01; 87 87 TimeEnd = 3; 88 #integration ="index0"; # "original";89 #DAESolver = "mebdf";90 88 Integration = "index0";# "index0"; # "original"; 89 DAESolver = "mebdf"; 90 NLASolver = "sundials"; # "nlasolver"; 91 91 RelativeAccuracy = 1e5; 92 92 AbsoluteAccuracy = 1e5; 93 #indVarAccuracy = 1e3;94 93 end
Note: See TracChangeset
for help on using the changeset viewer.