Changeset 198 for branches/newlanguage/sample
- 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)) = (l0-ll)*xl/ll53 eps^2*M/2 * diff(q(1))*'s' = (l0-ll)*xl/ll 54 54 +lambda(1)*xb +2*lambda(1)*(xl-xr); 55 55 56 eps^2*M/2 * diff(q(2)) = (l0-ll)*yl/ll56 eps^2*M/2 * diff(q(2))*'s' = (l0-ll)*yl/ll 57 57 +lambda(1)*yb +2*lambda(2)*(yr-yl) - eps^2*M/2; 58 58 59 59 60 eps^2*M/2 * diff(q(3)) = (l0-lr)*(xr-xb)/lr60 eps^2*M/2 * diff(q(3))*'s' = (l0-lr)*(xr-xb)/lr 61 61 -2*lambda(2)*(xl-xr); 62 62 63 eps^2*M/2 * diff(q(4)) = (l0-lr)*(yr-yb)/lr63 eps^2*M/2 * diff(q(4))*'s' = (l0-lr)*(yr-yb)/lr 64 64 -2*lambda(2)*(yl-yr) - eps^2*M/2; 65 65 … … 70 70 xb = sqrt(l^2-yb^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 = 1e-5; 92 92 AbsoluteAccuracy = 1e-5; 93 #indVarAccuracy = 1e-3;94 93 end
Note: See TracChangeset
for help on using the changeset viewer.