==SASDATA B34SEXEC DATA NOHEAD NOOB = 200 HEADING=('DATA FROM SAS ')$ INPUT X1 X2 X3 X4 Y ERROR PROBIT PROBITY TOBITY YGLS $ COMMENT=('FILE CREATED 9/ 9/91 AT TIME 10:34:16 BY B34S ')$ DATACARDS CARD72 $ 0.8400258296E+00 0.3141183718E+00 0.3874750235E+00 0.2927201628E+00 0.2065688681E+03 0.5054276064E+00 0.3397669941E+00 0.1000000000E+01 0.3397669941E+00 0.2051219440E+03 0.7477757656E+00 0.8672924652E+00 0.5844622203E+00 0.5653683236E-01 0.2244966336E+03 -0.3159042581E+00 -0.5983697468E-01 0.0000000000E+00 0.0000000000E+00 0.2239912060E+03 0.2145414525E+00 0.7981916829E+00 0.2076145210E+00 0.3772549569E+00 0.1952612828E+03 -0.5948131871E+00 -0.9972105739E+00 0.0000000000E+00 0.0000000000E+00 0.1955771870E+03 0.5240600728E+00 0.8776442721E+00 0.5672804364E+00 0.2822945802E+00 0.2116087488E+03 -0.4617912047E+00 -0.2105118431E+00 0.0000000000E+00 0.0000000000E+00 0.2122035620E+03 0.5250102156E+00 0.8466933235E+00 0.3746873834E+00 0.3708530902E+00 0.3248915721E+03 0.6456009732E+00 0.7628449858E+00 0.1000000000E+01 0.7628449858E+00 0.3253533633E+03 0.9278868758E+00 0.9947217717E+00 0.2888166853E+00 0.1420302182E+00 0.5260827413E+03 0.2295684846E+01 0.2649140397E+01 0.1000000000E+01 0.2649140397E+01 0.5254371403E+03 0.1018764549E+00 0.2375773950E+00 0.9632785716E+00 0.8229530393E+00 0.2166159852E+03 0.1053879550E+01 0.1179565011E+01 0.1000000000E+01 0.1179565011E+01 0.2143203003E+03 0.3717317604E+00 0.6956978156E+00 0.5931871780E+00 0.6969009726E+00 0.1562910111E+03 -0.6454332588E+00 -0.2879155320E+00 0.0000000000E+00 0.0000000000E+00 0.1552371315E+03 0.8146470063E+00 0.7722354992E+00 0.9620355465E+00 0.9314307808E+00 0.1494074056E+03 -0.7647578976E+00 0.7155909354E+00 0.1000000000E+01 0.7155909354E+00 0.1500528388E+03 0.5571328127E+00 0.7311827554E+00 0.9885706254E+00 0.9065015703E+00 0.3359920771E+03 0.1221602269E+01 0.2404990033E+01 0.1000000000E+01 0.2404990033E+01 0.3367568350E+03 0.5718912059E+00 0.7754967376E+00 0.7736682692E+00 0.4259979120E-01 0.6173194779E+02 -0.1660315752E+01 -0.1496659748E+01 0.0000000000E+00 0.0000000000E+00 0.6051034552E+02 0.9746906185E+00 0.6252244427E+00 0.1472082385E+00 0.1288642376E+00 0.6464767304E+02 -0.1644584886E+01 -0.1768597349E+01 0.0000000000E+00 0.0000000000E+00 0.6630798879E+02 0.8212408894E+00 0.5956277116E+00 0.7149492004E+00 0.1512110471E+00 0.1133179969E+03 -0.8403319393E+00 -0.5573030908E+00 0.0000000000E+00 0.0000000000E+00 0.1149625818E+03 0.4040682637E+00 0.1753082289E+00 0.4054037281E+00 0.6204587012E+00 0.3271214537E+02 -0.9206011625E+00 -0.1315362241E+01 0.0000000000E+00 0.0000000000E+00 0.3355247731E+02 0.4939039284E-01 0.1043325230E+00 0.5167133012E+00 0.4004530250E+00 0.3006079973E+02 -0.7108320196E+00 -0.1639942778E+01 0.0000000000E+00 0.0000000000E+00 0.3098140089E+02 0.4139913344E+00 0.9523569657E+00 0.2635233053E+00 0.3619218757E-01 0.3608267508E+03 0.7736598551E+00 0.4397236481E+00 0.1000000000E+01 0.4397236481E+00 0.3615375828E+03 0.2820965522E+00 0.1967527731E+00 0.8238580325E+00 0.5819517015E+00 0.1730189080E+03 0.6336352473E+00 0.5182943067E+00 0.1000000000E+01 0.5182943067E+00 0.1722452481E+03 0.8622479164E+00 0.8007310047E+00 0.8859956809E+00 0.9294094234E+00 0.4050661091E+03 0.1697481961E+01 0.3175865987E+01 0.1000000000E+01 0.3175865987E+01 0.4044324738E+03 0.5841787034E+00 0.2914675637E+00 0.6953433552E+00 0.6357705533E+00 0.1727444132E+03 0.3533002497E+00 0.5600604252E+00 0.1000000000E+01 0.5600604252E+00 0.1710469313E+03 0.3956885957E+00 0.3382281973E+00 0.6013120660E+00 0.2518926351E+00 0.1803809492E+03 0.3337537944E+00 -0.7912471153E-01 0.0000000000E+00 0.0000000000E+00 0.1800276489E+03 0.5595186823E+00 0.8304934963E+00 0.1041917303E+00 0.1504105567E+00 0.2396431988E+03 -0.2832249565E+00 -0.6386104910E+00 0.0000000000E+00 0.0000000000E+00 0.2393094450E+03 0.9502264583E+00 0.4560838952E+00 0.4020261713E+00 0.8538615494E+00 0.3093332854E+03 0.1216320594E+01 0.1878518668E+01 0.1000000000E+01 0.1878518668E+01 0.3096165103E+03 0.8510608249E+00 0.7792846723E+00 0.4374881187E+00 0.8628101884E+00 0.3439921507E+03 0.9412254613E+00 0.1871869266E+01 0.1000000000E+01 0.1871869266E+01 0.3427758301E+03 0.2508370900E+00 0.8189709499E+00 0.4447542080E+00 0.9839733928E+00 0.2565282479E+03 0.4430251046E-01 0.5428381510E+00 0.1000000000E+01 0.5428381510E+00 0.2555870225E+03 0.6408119672E+00 0.1267326935E+00 0.9963799813E+00 0.1583454250E+00 0.1577939282E+03 0.6809188192E+00 0.6031888862E+00 0.1000000000E+01 0.6031888862E+00 0.1577496257E+03 0.3115585768E+00 0.3649997129E+00 0.5501748868E+00 0.7893228917E+00 0.2216114573E+03 0.6520683467E+00 0.6681244150E+00 0.1000000000E+01 0.6681244150E+00 0.2209305385E+03 0.1498414450E+00 0.3851666613E+00 0.4960772791E+00 0.5708297431E+00 0.1496429891E+03 -0.1041158333E+00 -0.5022007048E+00 0.0000000000E+00 0.0000000000E+00 0.1489909208E+03 0.9354915046E+00 0.8057180359E+00 0.7030294188E+00 0.8154422729E+00 0.3755191036E+03 0.1311736350E+01 0.2571417582E+01 0.1000000000E+01 0.2571417582E+01 0.3756232194E+03 0.1382810944E+00 0.9035317739E-01 0.5658524730E+00 0.2825141191E+00 0.1319160713E+03 0.3538171058E+00 -0.5691820303E+00 0.0000000000E+00 0.0000000000E+00 0.1306043349E+03 0.2147989594E+00 0.1261112704E+00 0.5521215149E+00 0.5063012629E+00 0.3913121942E+01 -0.1022474029E+01 -0.1623141021E+01 0.0000000000E+00 0.0000000000E+00 0.3559304836E+01 0.4053252512E+00 0.3014963432E+00 0.2490405240E+00 0.6240870296E+00 0.1135089734E+03 -0.4325524037E+00 -0.8526032557E+00 0.0000000000E+00 0.0000000000E+00 0.1145314474E+03 0.3070652718E-01 0.8460226795E-01 0.9103173860E+00 0.7043070882E+00 0.9734267927E+01 -0.7187113878E+00 -0.9887781184E+00 0.0000000000E+00 0.0000000000E+00 0.1016682033E+02 0.2892307513E+00 0.1012377604E+00 0.5030388930E+00 0.5746744513E+00 0.3837923420E+02 -0.6600332525E+00 -0.1191851397E+01 0.0000000000E+00 0.0000000000E+00 0.3909794559E+02 0.5535022938E+00 0.7130513801E+00 0.2545457431E+00 0.1503044753E+00 0.8603107969E+02 -0.1519202747E+01 -0.1847798855E+01 0.0000000000E+00 0.0000000000E+00 0.8669111295E+02 0.1673169308E+00 0.9565653047E-01 0.6993075409E+00 0.2618398606E+00 0.1156598661E+03 0.2361641569E+00 -0.5397149803E+00 0.0000000000E+00 0.0000000000E+00 0.1171790689E+03 0.7425368452E+00 0.8167564585E+00 0.2257977971E+00 0.9835756510E+00 0.2589216033E+03 -0.7063629786E-01 0.6980304538E+00 0.1000000000E+01 0.6980304538E+00 0.2586854392E+03 0.9559656256E+00 0.9142692550E+00 0.1233689287E+00 0.4615842341E+00 0.3016916394E+03 0.1227507487E+00 0.5779387922E+00 0.1000000000E+01 0.5779387922E+00 0.3017622757E+03 0.8462232673E+00 0.4744528546E+00 0.1291272589E+00 0.2418403416E+00 0.2362400563E+03 0.3723052315E+00 0.6394895388E-01 0.1000000000E+01 0.6394895388E-01 0.2361173056E+03 0.6106211715E+00 0.7100300356E+00 0.4748080173E+00 0.9834604850E-01 0.3297873388E+03 0.1016001345E+01 0.9098066177E+00 0.1000000000E+01 0.9098066177E+00 0.3294150336E+03 0.9020371460E+00 0.5383134200E+00 0.4336503395E+00 0.3612557738E+00 0.3230712755E+03 0.1232289058E+01 0.1467545737E+01 0.1000000000E+01 0.1467545737E+01 0.3220552742E+03 0.6257904046E+00 0.6593304737E+00 0.3672722676E+00 0.7450012802E+00 0.3105309737E+03 0.8447467601E+00 0.1242141186E+01 0.1000000000E+01 0.1242141186E+01 0.3092986846E+03 0.2365165674E+00 0.1339485063E+00 0.2725457498E+00 0.6764166742E+00 0.1162537829E+03 -0.4563700178E-01 -0.7262095041E+00 0.0000000000E+00 0.0000000000E+00 0.1154090361E+03 0.5350426608E+00 0.4619996711E+00 0.8284716406E+00 0.1228634058E+00 0.1182174194E+03 -0.4452297791E+00 -0.4968524008E+00 0.0000000000E+00 0.0000000000E+00 0.1182630564E+03 0.9652621173E+00 0.1604059651E+00 0.9430547007E+00 0.9203552692E+00 0.2059893601E+03 0.1002051217E+01 0.1991129270E+01 0.1000000000E+01 0.1991129270E+01 0.2064345899E+03 0.4110097759E+00 0.8413042291E+00 0.8001785972E+00 0.6016834020E+00 0.2302840463E+03 -0.1068763458E+00 0.5472996584E+00 0.1000000000E+01 0.5472996584E+00 0.2292819951E+03 0.4929368817E+00 0.7901711478E+00 0.4064809033E+00 0.7245412855E+00 0.2398101340E+03 -0.9297491943E-01 0.3211552989E+00 0.1000000000E+01 0.3211552989E+00 0.2399170104E+03 0.3653860861E+00 0.4394824943E-01 0.6382280984E+00 0.6996497451E+00 0.1229237289E+03 0.3442577770E+00 0.9146995601E-01 0.1000000000E+01 0.9146995601E-01 0.1230167039E+03 0.1326644653E-01 0.9691669061E+00 0.7881905683E+00 0.1188819497E+00 0.1921975221E+03 -0.6847073893E+00 -0.7952015186E+00 0.0000000000E+00 0.0000000000E+00 0.1918532643E+03 0.4892842704E-01 0.3400733336E+00 0.6125176051E+00 0.5833884119E+00 0.2512843472E+02 -0.1199542013E+01 -0.1614634235E+01 0.0000000000E+00 0.0000000000E+00 0.2581314211E+02 0.9038897697E-02 0.9167535915E+00 0.8776124069E+00 0.3172198823E-01 0.1209499738E+03 -0.1249124100E+01 -0.1413997215E+01 0.0000000000E+00 0.0000000000E+00 0.1221495158E+03 0.1514561349E+00 0.5232585832E+00 0.4070070782E+00 0.5679628788E+00 0.7724308759E+01 -0.1837804793E+01 -0.2188120117E+01 0.0000000000E+00 0.0000000000E+00 0.8973432858E+01 0.7521048211E+00 0.6257290556E+00 0.6282376298E+00 0.7898436277E+00 -0.8324092235E+02 -0.2928427817E+01 -0.2132512683E+01 0.0000000000E+00 0.0000000000E+00 -0.8140311756E+02 0.9018510454E+00 0.4105205997E+00 0.6197189786E+00 0.6168741225E+00 -0.2750313487E+02 -0.1950622198E+01 -0.1401657452E+01 0.0000000000E+00 0.0000000000E+00 -0.2457470705E+02 0.8033773028E+00 0.3623286500E+00 0.6576199480E+00 0.6184665947E+00 0.2265635357E+03 0.7124935752E+00 0.1154286071E+01 0.1000000000E+01 0.1154286071E+01 0.2285141579E+03 0.5680574947E+00 0.3423134848E+00 0.2627382674E+00 0.8420606195E+00 0.9595516336E+02 -0.7110066613E+00 -0.6958367950E+00 0.0000000000E+00 0.0000000000E+00 0.9524266979E+02 0.5128319287E+00 0.1662252261E+00 0.7473755726E+00 0.1412481158E+00 0.1592136148E+03 0.5227115934E+00 0.9039243662E-01 0.1000000000E+01 0.9039243662E-01 0.1599246215E+03 0.9570823186E+00 0.6825288399E+00 0.2622124054E+00 0.3898346331E-02 0.2662218681E+03 0.3140091868E+00 0.2197310971E+00 0.1000000000E+01 0.2197310971E+00 0.2656991566E+03 0.5195067846E+00 0.3505291357E+00 0.3431832143E+00 0.8802829175E+00 0.3038625640E+03 0.1389171327E+01 0.1482673379E+01 0.1000000000E+01 0.1482673379E+01 0.3035485549E+03 0.9149948531E+00 0.3184957667E+00 0.9583513536E+00 0.1120077958E-01 -0.7038705771E+02 -0.2020830553E+01 -0.1817787800E+01 0.0000000000E+00 0.0000000000E+00 -0.7177622903E+02 0.2515024566E+00 0.1788230147E-02 0.5478408420E-01 0.7561032264E+00 -0.5155058758E+01 -0.1094525299E+01 -0.2030347301E+01 0.0000000000E+00 0.0000000000E+00 -0.3134228205E+01 0.8269269341E+00 0.1609808668E+00 0.6054275730E+00 0.4212188958E+00 -0.6070336080E+01 -0.1226084876E+01 -0.1211530607E+01 0.0000000000E+00 0.0000000000E+00 -0.4975810782E+01 0.4259820306E+00 0.4799878611E+00 0.1559818932E+00 0.5876787691E+00 0.6093951188E+02 -0.1355490531E+01 -0.1705859977E+01 0.0000000000E+00 0.0000000000E+00 0.6216559675E+02 0.1170728975E+00 0.6441882416E+00 0.8717760243E+00 0.9396411716E+00 -0.9352888760E+02 -0.2907491017E+01 -0.2334812682E+01 0.0000000000E+00 0.0000000000E+00 -0.9217339706E+02 0.5491710925E+00 0.9185508802E+00 0.8464285219E-01 0.5924167722E+00 0.2402778887E+03 -0.4824643947E+00 -0.3376827978E+00 0.0000000000E+00 0.0000000000E+00 0.2431853797E+03 0.7486895793E+00 0.2257594927E+00 0.3397941687E+00 0.9205935835E+00 0.3009136657E+03 0.1582830527E+01 0.1817667351E+01 0.1000000000E+01 0.1817667351E+01 0.3013961301E+03 0.4163577282E+00 0.7243376388E+00 0.9426956722E+00 0.8861619621E+00 0.1579695656E+03 -0.5495643535E+00 0.4199886477E+00 0.1000000000E+01 0.4199886477E+00 0.1563867351E+03 0.7240963154E+00 0.8867737218E+00 0.5941967014E-02 0.8666396081E+00 0.3435896599E+03 0.5491625893E+00 0.1032614202E+01 0.1000000000E+01 0.1032614202E+01 0.3441392242E+03 0.6118925237E+00 0.7764598867E-01 0.9961315640E+00 0.9831965654E+00 0.4222638988E+02 -0.4150405867E+00 0.2538260552E+00 0.1000000000E+01 0.2538260552E+00 0.4167722729E+02 0.5846749418E+00 0.6317475303E+00 0.7807418802E+00 0.9287802218E+00 0.2101075504E+03 0.6839294668E-01 0.9943375209E+00 0.1000000000E+01 0.9943375209E+00 0.2105225910E+03 0.9188564960E-02 0.4322112813E+00 0.1750042387E+00 0.2962403088E+00 -0.3753325627E+02 -0.2180234174E+01 -0.3267589781E+01 0.0000000000E+00 0.0000000000E+00 -0.3760164921E+02 0.9108697944E+00 0.9886350199E+00 0.9887788692E+00 0.4064539014E+00 0.1545495910E+03 -0.1118008890E+01 0.1767286945E+00 0.1000000000E+01 0.1767286945E+00 0.1567298252E+03 0.2707214379E+00 0.1520636817E-01 0.5734299072E+00 0.6364504187E+00 0.6655761640E+02 -0.1771563769E+00 -0.6813482450E+00 0.0000000000E+00 0.0000000000E+00 0.6767562529E+02 0.8221877092E+00 0.5088277920E+00 0.8687003949E+00 0.2475366234E+00 0.1387430706E+03 -0.3512793104E+00 0.9597320901E-01 0.1000000000E+01 0.9597320901E-01 0.1389202270E+03 0.3480292284E+00 0.3272416010E+00 0.9495872804E+00 0.7134216218E+00 0.2248391481E+03 0.9317568055E+00 0.1270036537E+01 0.1000000000E+01 0.1270036537E+01 0.2251904274E+03 0.4771980220E+00 0.2671549483E+00 0.7321641827E-01 0.5483418643E+00 0.7677942142E+02 -0.8101695179E+00 -0.1444258265E+01 0.0000000000E+00 0.0000000000E+00 0.7584766462E+02 0.9817125225E+00 0.6423648312E+00 0.2257175692E+00 0.6351861446E+00 0.2147348222E+03 -0.1702534447E+00 0.3147276228E+00 0.1000000000E+01 0.3147276228E+00 0.2155449918E+03 0.5735325688E+00 0.3618830903E+00 0.1690994003E+00 0.5362060061E-01 0.7280212631E+02 -0.9830664624E+00 -0.1824930802E+01 0.0000000000E+00 0.0000000000E+00 0.7297237976E+02 0.2014344060E+00 0.5080614104E+00 0.9881244847E+00 0.4082137137E+00 0.2269252531E+03 0.6374691123E+00 0.7433031270E+00 0.1000000000E+01 0.7433031270E+00 0.2279083196E+03 0.8478856855E+00 0.4147163196E+00 0.1371843019E+00 0.6565613090E+00 0.1407960336E+03 -0.4800996870E+00 -0.4237520710E+00 0.0000000000E+00 0.0000000000E+00 0.1401585645E+03 0.8259204649E+00 0.2452541902E+00 0.9871742376E+00 0.4374115556E+00 0.9075108678E+01 -0.1079734996E+01 -0.5839745474E+00 0.0000000000E+00 0.0000000000E+00 0.9555208365E+01 0.5760142081E+00 0.7079623779E-01 0.8723685839E+00 0.8987899660E+00 0.1904080423E+03 0.1035065519E+01 0.1453034515E+01 0.1000000000E+01 0.1453034515E+01 0.1914877773E+03 0.9629587852E+00 0.4483035221E+00 0.6372955390E+00 0.2612439683E-01 0.2828004473E+03 0.1107832412E+01 0.1182514655E+01 0.1000000000E+01 0.1182514655E+01 0.2817653818E+03 0.7273747260E-01 0.4987019463E+00 0.6836122074E+00 0.4703700482E+00 0.2003600717E+03 0.2693578243E+00 -0.5220501154E-02 0.0000000000E+00 0.0000000000E+00 0.1992522393E+03 0.5093999528E+00 0.4850072577E+00 0.5169795998E+00 0.8761344641E+00 0.2607096931E+03 0.7646369251E+00 0.1152158200E+01 0.1000000000E+01 0.1152158200E+01 0.2604403353E+03 0.1919373148E+00 0.8904500469E+00 0.7939376333E+00 0.7098025321E+00 0.2511554376E+03 0.1736360613E-01 0.6034911332E+00 0.1000000000E+01 0.6034911332E+00 0.2503908006E+03 0.6511570537E+00 0.9966007485E+00 0.8687801310E+00 0.5876614838E+00 0.3080271781E+03 0.3661469611E+00 0.1470346378E+01 0.1000000000E+01 0.1470346378E+01 0.3080098145E+03 0.8265580707E+00 0.9614946409E+00 0.8404298224E+00 0.1040258943E+00 0.2273933886E+03 -0.3755276741E+00 0.3569807543E+00 0.1000000000E+01 0.3569807543E+00 0.2270272416E+03 0.3632056736E+00 0.3977568462E+00 0.9931334341E-01 0.1593626826E+00 0.2713333056E+03 0.9162353998E+00 -0.6412605442E-01 0.0000000000E+00 0.0000000000E+00 0.2717088333E+03 0.4086065763E+00 0.4507273301E+00 0.3742362359E+00 0.7884161114E+00 0.2750348931E+03 0.9295343893E+00 0.9515206429E+00 0.1000000000E+01 0.9515206429E+00 0.2741186577E+03 0.9095838987E+00 0.3765859022E+00 0.2792579480E+00 0.4883327871E+00 0.1271213277E+03 -0.4772526392E+00 -0.4234921032E+00 0.0000000000E+00 0.0000000000E+00 0.1261917933E+03 0.4091521674E+00 0.6204772159E+00 0.3605668345E+00 0.4678665986E-01 0.1670905895E+03 -0.4582593475E+00 -0.1021276470E+01 0.0000000000E+00 0.0000000000E+00 0.1675678421E+03 0.3433923149E+00 0.3946357681E+00 0.6433551558E+00 0.8701040237E+00 0.3194817670E+03 0.1604344418E+01 0.1855831680E+01 0.1000000000E+01 0.1855831680E+01 0.3199400263E+03 0.8383266800E+00 0.7565108999E+00 0.6786940792E+00 0.8113885163E+00 0.5896781347E+02 -0.1755907305E+01 -0.6709871295E+00 0.0000000000E+00 0.0000000000E+00 0.5736346905E+02 0.6794007961E-02 0.1868917956E+00 0.9040867770E-01 0.4986461191E+00 0.2404106427E+03 0.1043586860E+01 -0.1736725393E+00 0.0000000000E+00 0.0000000000E+00 0.2421665500E+03 0.7453244501E+00 0.6680324837E+00 0.6219529309E+00 0.1629101573E+00 0.1406811023E+03 -0.7444921371E+00 -0.5462721150E+00 0.0000000000E+00 0.0000000000E+00 0.1396375154E+03 0.3101456398E-01 0.2617768361E+00 0.6832847246E+00 0.9663662999E+00 0.1821826902E+03 0.5406658897E+00 0.4831083143E+00 0.1000000000E+01 0.4831083143E+00 0.1829271823E+03 0.7184023135E+00 0.1876835688E+00 0.3977408507E+00 0.8304775595E+00 0.6894833217E+02 -0.6282193595E+00 -0.4939150670E+00 0.0000000000E+00 0.0000000000E+00 0.6840766628E+02 0.8363420689E+00 0.4011512413E+00 0.1489120075E+00 0.7641092925E+00 0.3222385864E+03 0.1362275873E+01 0.1512790483E+01 0.1000000000E+01 0.1512790483E+01 0.3228668058E+03 0.3848796063E+00 0.6715432568E+00 0.6275177994E+00 0.6916541993E+00 0.2112125933E+03 -0.3419859789E-01 0.3413962639E+00 0.1000000000E+01 0.3413962639E+00 0.2098503174E+03 0.6321271085E+00 0.1603128626E+00 0.3782817372E+00 0.7811573236E+00 0.1395922621E+03 0.1356874658E+00 0.8756649778E-01 0.1000000000E+01 0.8756649778E-01 0.1396264607E+03 0.9111381643E+00 0.4991281494E+00 0.8468062379E+00 0.2724410533E+00 0.1480941552E+03 -0.2579239337E+00 0.2715896712E+00 0.1000000000E+01 0.2715896712E+00 0.1479584677E+03 0.9167827689E+00 0.3679970517E+00 0.9264476248E+00 0.8052301197E+00 0.1577929997E+03 0.1083685895E+00 0.1124826155E+01 0.1000000000E+01 0.1124826155E+01 0.1580509236E+03 0.5026223236E+00 0.5733923479E+00 0.5191127304E-02 0.2472765981E+00 0.1262781005E+03 -0.9443975682E+00 -0.1615915171E+01 0.0000000000E+00 0.0000000000E+00 0.1261697319E+03 0.9777846360E+00 0.6263766601E+00 0.5125258812E+00 0.2248490741E-01 0.1607906268E+03 -0.5233636322E+00 -0.3841915476E+00 0.0000000000E+00 0.0000000000E+00 0.1617350244E+03 0.9038389199E+00 0.8207274540E+00 0.9663199284E+00 0.9390368103E+00 0.1568459612E+03 -0.7948134593E+00 0.8351096534E+00 0.1000000000E+01 0.8351096534E+00 0.1573693248E+03 0.3916707956E+00 0.8110620393E+00 0.5196952173E+00 0.5175164132E+00 0.1668790082E+03 -0.7949079538E+00 -0.5549634884E+00 0.0000000000E+00 0.0000000000E+00 0.1676738216E+03 0.8983572894E+00 0.6909625693E+00 0.7902602199E-02 0.8190351664E+00 0.2710614362E+03 0.2012998551E+00 0.6175574824E+00 0.1000000000E+01 0.6175574824E+00 0.2718563442E+03 0.5240418338E+00 0.5711010410E+00 0.4951968440E+00 0.7733568638E+00 0.1590491627E+03 -0.4298478387E+00 -0.6615125606E-01 0.0000000000E+00 0.0000000000E+00 0.1588478629E+03 0.8088092678E+00 0.6573642953E+00 0.3217116096E+00 0.7022642996E-02 0.2571926721E+03 0.3143020635E+00 0.1092098792E+00 0.1000000000E+01 0.1092098792E+00 0.2576225199E+03 0.2956083650E-01 0.8289790167E+00 0.6503344723E+00 0.1714759856E+00 0.1148382183E+03 -0.1241461910E+01 -0.1561111599E+01 0.0000000000E+00 0.0000000000E+00 0.1145239162E+03 0.9968891912E+00 0.7166367336E+00 0.5135819116E+00 0.7711877166E+00 0.1452176305E+03 -0.8985052456E+00 0.9979030742E-01 0.1000000000E+01 0.9979030742E-01 0.1464590924E+03 0.3519533896E+00 0.2806188180E+00 0.3604743131E+00 0.4917810052E+00 0.2624098829E+03 0.1158080758E+01 0.6429082842E+00 0.1000000000E+01 0.6429082842E+00 0.2633083881E+03 0.3633538955E+00 0.8889214857E+00 0.1034109690E+00 0.2815636387E-01 0.3152892943E+03 0.3817734804E+00 -0.2343838055E+00 0.0000000000E+00 0.0000000000E+00 0.3141312135E+03 0.2240075790E+00 0.8953807964E+00 0.6650450774E+00 0.4126166492E+00 0.3019226313E+03 0.4714025129E+00 0.6684526150E+00 0.1000000000E+01 0.6684526150E+00 0.3015408578E+03 0.8480227166E+00 0.7177976533E+00 0.2515874291E-01 0.8429920613E+00 0.4332886620E+02 -0.2118440260E+01 -0.1684469086E+01 0.0000000000E+00 0.0000000000E+00 0.4285746369E+02 0.1675739908E+00 0.4160628035E+00 0.7675379761E+00 0.1076514414E-01 0.2478801186E+03 0.9594212524E+00 0.3213611669E+00 0.1000000000E+01 0.3213611669E+00 0.2499985588E+03 0.9297774904E+00 0.7702818378E+00 0.1268476942E+00 0.9291957854E+00 0.2551542863E+03 -0.7391384196E-01 0.6821889659E+00 0.1000000000E+01 0.6821889659E+00 0.2541948651E+03 0.9935657107E+00 0.8589001908E+00 0.5355067202E+00 0.2614462166E+00 0.2650292719E+03 0.5033134499E-01 0.6997501833E+00 0.1000000000E+01 0.6997501833E+00 0.2651031857E+03 0.1265625530E+00 0.1368281408E+00 0.6705627837E+00 0.1487060129E+00 0.5171265896E+02 -0.4882792510E+00 -0.1405619761E+01 0.0000000000E+00 0.0000000000E+00 0.5166232762E+02 0.3019583249E+00 0.1356580993E-01 0.5675205030E-03 0.5383170932E+00 0.6403663779E+02 -0.4436328953E+00 -0.1589224147E+01 0.0000000000E+00 0.0000000000E+00 0.6452491704E+02 0.4953847402E+00 0.9313281639E+00 0.8324514725E+00 0.1189872017E-01 0.2780744287E+03 0.2259086818E+00 0.4969717786E+00 0.1000000000E+01 0.4969717786E+00 0.2785180616E+03 0.9817899191E+00 0.9431700124E+00 0.8583990880E+00 0.1134725004E+00 0.8963280928E+02 -0.1724752941E+01 -0.8279214212E+00 0.0000000000E+00 0.0000000000E+00 0.8940690060E+02 0.1323133879E+00 0.7911110087E+00 0.2027237215E+00 0.1775873839E+00 0.3101208873E+03 0.5840473478E+00 -0.1122171502E+00 0.0000000000E+00 0.0000000000E+00 0.3118456402E+03 0.7111610794E+00 0.4842606499E+00 0.9687426337E+00 0.6574446647E+00 0.1811006381E+03 0.1550560721E+00 0.9766650998E+00 0.1000000000E+01 0.9766650998E+00 0.1805165907E+03 0.6724803322E+00 0.3769432767E+00 0.2856519916E+00 0.9530229806E+00 0.1522074515E+03 -0.2238808640E+00 0.6421771710E-01 0.1000000000E+01 0.6421771710E-01 0.1520523955E+03 0.4572348916E+00 0.7468230738E+00 0.8554018685E+00 0.7392033129E+00 0.1841843051E+03 -0.3666639478E+00 0.4319991990E+00 0.1000000000E+01 0.4319991990E+00 0.1844081860E+03 0.7900799089E+00 0.8730290695E+00 0.9995712438E+00 0.7938939868E+00 0.4070887216E+03 0.1635942022E+01 0.3092516231E+01 0.1000000000E+01 0.3092516231E+01 0.4074553855E+03 0.9762369418E+00 0.6142801358E+00 0.2062430033E+00 0.3261556714E+00 0.1125771787E+03 -0.1128035471E+01 -0.1005119719E+01 0.0000000000E+00 0.0000000000E+00 0.1109412366E+03 0.6983689540E+00 0.4870099886E+00 0.1768789389E+00 0.8043263792E+00 0.2161132035E+03 0.1531167874E+00 0.3197010482E+00 0.1000000000E+01 0.3197010482E+00 0.2172412390E+03 0.3134553158E+00 0.2434921875E+00 0.3731947490E+00 0.2841464668E+00 0.1565361308E+03 0.1932978204E+00 -0.5924134606E+00 0.0000000000E+00 0.0000000000E+00 0.1563830140E+03 0.6496668034E+00 0.9499651357E+00 0.6403503800E-01 0.2368836342E+00 0.2945970540E+03 -0.3235526608E-02 -0.1026849153E+00 0.0000000000E+00 0.0000000000E+00 0.2944037562E+03 0.3032405996E+00 0.5647570803E+00 0.8722483166E+00 0.8774579256E+00 0.9127649321E+02 -0.9158794086E+00 -0.2981754865E+00 0.0000000000E+00 0.0000000000E+00 0.9127972873E+02 0.4353552931E+00 0.1641111216E-01 0.8215620019E+00 0.9925667285E+00 0.1321287146E+03 0.5485727170E+00 0.8144678527E+00 0.1000000000E+01 0.8144678527E+00 0.1330445940E+03 0.6900587911E-01 0.7818101658E+00 0.8834574171E+00 0.2688097744E+00 0.3568115421E+03 0.1364040702E+01 0.1367123939E+01 0.1000000000E+01 0.1367123939E+01 0.3562629694E+03 0.8858787590E+00 0.9643020197E+00 0.2404479170E-01 0.1208141386E+00 0.3562703622E+03 0.5498102599E+00 0.5448499689E+00 0.1000000000E+01 0.5448499689E+00 0.3549063215E+03 0.5232267403E+00 0.8718242957E+00 0.7509382614E+00 0.1935923984E-01 0.1101229458E+03 -0.1372806317E+01 -0.1207457780E+01 0.0000000000E+00 0.0000000000E+00 0.1095731355E+03 0.3707440370E+00 0.9502945980E-01 0.1601308682E+00 0.3195018410E+00 0.4285899906E+02 -0.7456621515E+00 -0.1800255946E+01 0.0000000000E+00 0.0000000000E+00 0.4423180537E+02 0.8674412691E+00 0.8540982148E-01 0.4828695604E+00 0.5887010484E+00 0.1007582833E+03 -0.7178207882E-01 -0.4736037948E-01 0.0000000000E+00 0.0000000000E+00 0.1015039454E+03 0.2985204631E+00 0.2334225341E+00 0.1325301924E+00 0.4349433123E+00 0.3374490174E+03 0.1913033876E+01 0.1012450378E+01 0.1000000000E+01 0.1012450378E+01 0.3375207994E+03 0.9225005144E-01 0.4466145947E+00 0.2514925250E+00 0.8348674267E+00 0.2530151881E+03 0.6940961007E+00 0.3193206985E+00 0.1000000000E+01 0.3193206985E+00 0.2511021542E+03 0.6168398731E+00 0.2277478441E+00 0.7580152252E+00 0.9618895291E+00 0.3002630501E+03 0.1763302896E+01 0.2327795368E+01 0.1000000000E+01 0.2327795368E+01 0.2995689540E+03 0.4773161733E+00 0.2529239669E+00 0.8931108662E+00 0.5143281012E+00 0.9711720739E+02 -0.2240862097E+00 -0.8640710223E-01 0.0000000000E+00 0.0000000000E+00 0.9535390449E+02 0.3123966032E+00 0.4497094913E+00 0.2674209275E+00 0.5435289680E+00 0.2850854455E+03 0.1008010362E+01 0.5810663524E+00 0.1000000000E+01 0.5810663524E+00 0.2853095317E+03 0.9136522049E-01 0.5752607896E+00 0.4080907949E+00 0.7819902985E+00 0.2148207171E+03 0.1249285962E+00 -0.1836430020E-01 0.0000000000E+00 0.0000000000E+00 0.2138127067E+03 0.9109474891E+00 0.2944490775E+00 0.8056447604E+00 0.4714876055E+00 0.7399285986E+02 -0.6172114385E+00 -0.1346825061E+00 0.0000000000E+00 0.0000000000E+00 0.7386793127E+02 0.2921852791E+00 0.7579854145E+00 0.4608610221E+00 0.6911987453E+00 0.4181987393E+03 0.1800408339E+01 0.2002638800E+01 0.1000000000E+01 0.2002638800E+01 0.4188159508E+03 0.9773118668E+00 0.6805457040E+00 0.9316466847E+00 0.1858289181E+00 0.2416699573E+03 0.3491936070E+00 0.1124526781E+01 0.1000000000E+01 0.1124526781E+01 0.2398695490E+03 0.2266269349E+00 0.9188956930E+00 0.8799127587E+00 0.6937356680E+00 -0.3372943286E+02 -0.2854072705E+01 -0.2134901651E+01 0.0000000000E+00 0.0000000000E+00 -0.3407862647E+02 0.6153715723E+00 0.5500159725E+00 0.1184491613E+00 0.7750538825E+00 -0.4271642609E+02 -0.2576552918E+01 -0.2517662329E+01 0.0000000000E+00 0.0000000000E+00 -0.3986235339E+02 0.3306032741E+00 0.4492269768E+00 0.1577997156E+00 0.1398207932E+00 0.2684462899E+03 0.8138114558E+00 -0.1087377845E+00 0.0000000000E+00 0.0000000000E+00 0.2710228428E+03 0.9680709415E+00 0.3683139845E+00 0.2531371512E+00 0.4761009752E+00 0.2358634474E+03 0.6102433376E+00 0.6758663900E+00 0.1000000000E+01 0.6758663900E+00 0.2350496360E+03 0.8290898823E+00 0.5136521175E+00 0.9511381038E+00 0.7781113475E+00 0.7981165888E+02 -0.9420128170E+00 0.1299786340E+00 0.1000000000E+01 0.1299786340E+00 0.7920141554E+02 0.7174168103E+00 0.6243313416E+00 0.1368581154E+00 0.1743461727E+00 0.2362186000E+03 0.9191331665E-01 -0.2551342433E+00 0.0000000000E+00 0.0000000000E+00 0.2371606128E+03 0.2361248457E+00 0.5502814665E+00 0.5806066343E+00 0.2557020574E+00 0.2702677184E+03 0.8153775160E+00 0.4380925198E+00 0.1000000000E+01 0.4380925198E+00 0.2701758050E+03 0.5844785147E+00 0.3303971454E+00 0.9848220148E+00 0.9036029517E+00 0.9022302359E+02 -0.4387185875E+00 0.3645820390E+00 0.1000000000E+01 0.3645820390E+00 0.8940764607E+02 0.8548087501E+00 0.7706628754E+00 0.5309465178E+00 0.6181251495E+00 0.1957342120E+03 -0.4720637607E+00 0.3024795321E+00 0.1000000000E+01 0.3024795321E+00 0.1961729306E+03 0.8293878906E+00 0.5222779973E+00 0.9263013103E+00 0.3461222886E+00 0.2461000551E+03 0.7145092162E+00 0.1338598703E+01 0.1000000000E+01 0.1338598703E+01 0.2465721189E+03 0.2773046815E+00 0.6597816958E+00 0.9509614655E+00 0.8093503504E+00 0.2900807280E+03 0.9219593329E+00 0.1619357526E+01 0.1000000000E+01 0.1619357526E+01 0.2893662188E+03 0.7513392740E+00 0.7591779991E+00 0.5046309629E+00 0.3325926601E+00 0.2755726464E+03 0.3625982196E+00 0.7103391156E+00 0.1000000000E+01 0.7103391156E+00 0.2746506871E+03 0.8848382258E+00 0.4760607793E+00 0.1535176882E+00 0.1717848866E+00 0.1219216184E+03 -0.7639658359E+00 -0.1077764256E+01 0.0000000000E+00 0.0000000000E+00 0.1215590202E+03 0.1885895007E+00 0.6237387522E+00 0.1772084334E+00 0.3421400427E+00 0.2753373581E+03 0.5461297508E+00 -0.1221935202E+00 0.0000000000E+00 0.0000000000E+00 0.2761013239E+03 0.3476974689E+00 0.7513604774E+00 0.1155438968E+00 0.9462734847E+00 0.2145574994E+03 -0.3895455059E+00 -0.2286701780E+00 0.0000000000E+00 0.0000000000E+00 0.2140113697E+03 0.1845790679E-01 0.2220393582E+00 0.8154936829E+00 0.2327912953E-02 0.7631202214E+02 -0.3322583977E+00 -0.1273939537E+01 0.0000000000E+00 0.0000000000E+00 0.7670156765E+02 0.1252329969E+00 0.7909786417E+00 0.9780313023E+00 0.7720983713E+00 0.1052971930E+03 -0.1159560112E+01 -0.4932187993E+00 0.0000000000E+00 0.0000000000E+00 0.1056294514E+03 0.6573263281E+00 0.6835964187E+00 0.2050091722E+00 0.5891572475E+00 0.3405727821E+03 0.1031498433E+01 0.1166587599E+01 0.1000000000E+01 0.1166587599E+01 0.3417323423E+03 0.9658578956E+00 0.1736506238E+00 0.5460336157E+00 0.1869792729E+00 0.1783239410E+03 0.5647978644E+00 0.4373192723E+00 0.1000000000E+01 0.4373192723E+00 0.1772924426E+03 0.5606391006E+00 0.6613640352E+00 0.5453398624E+00 0.5270665227E+00 0.1708016217E+03 -0.4626329491E+00 -0.1682234281E+00 0.0000000000E+00 0.0000000000E+00 0.1702368238E+03 0.4070477509E+00 0.2515492751E+00 0.7886670479E+00 0.1270734160E+00 0.3275080842E+03 0.2064050676E+01 0.1638388166E+01 0.1000000000E+01 0.1638388166E+01 0.3279707171E+03 0.7229030603E+00 0.8317338093E+00 0.9501332142E+00 0.8889317154E+00 0.3717591996E+03 0.1345944768E+01 0.2739646567E+01 0.1000000000E+01 0.2739646567E+01 0.3696951489E+03 0.2753406488E+00 0.6502849474E+00 0.3391103097E+00 0.4269755273E+00 0.2073722216E+03 -0.1299130868E+00 -0.4382016536E+00 0.0000000000E+00 0.0000000000E+00 0.2060262768E+03 0.1776875430E+00 0.3945348426E+00 0.9470991888E+00 0.8960663020E+00 0.2347070258E+03 0.9026811541E+00 0.1318069030E+01 0.1000000000E+01 0.1318069030E+01 0.2348369388E+03 0.1863371753E+00 0.7689054920E+00 0.9946046122E+00 0.3197178333E+00 0.1567851529E+03 -0.5768990810E+00 -0.3073339681E+00 0.0000000000E+00 0.0000000000E+00 0.1558824717E+03 0.4976243847E+00 0.5730337366E+00 0.9780102470E+00 0.4182212590E+00 0.2001927313E+03 0.2057307443E+00 0.6726203716E+00 0.1000000000E+01 0.6726203716E+00 0.2007696303E+03 0.4470041627E-01 0.2798962607E+00 0.2164536663E+00 0.9367700606E+00 0.1022483386E+03 -0.4955426035E+00 -0.1017722200E+01 0.0000000000E+00 0.0000000000E+00 0.1020426079E+03 0.2944079214E+00 0.1139345928E+00 0.8987013497E+00 0.4735850631E+00 0.5669915832E+02 -0.3275560673E+00 -0.5469271403E+00 0.0000000000E+00 0.0000000000E+00 0.5719470092E+02 0.5441559576E+00 0.6291801049E+00 0.6300236830E+00 0.8080394357E+00 0.2474345569E+03 0.3920779756E+00 0.1003477157E+01 0.1000000000E+01 0.1003477157E+01 0.2477621130E+03 0.7187949889E+00 0.7873790129E+00 0.4790699580E+00 0.7287836908E+00 0.2801946185E+03 0.3248695579E+00 0.1038897209E+01 0.1000000000E+01 0.1038897209E+01 0.2798025405E+03 0.6674918854E+00 0.5361173763E+00 0.5247433730E+00 0.3618706851E+00 0.2448899689E+03 0.5174618637E+00 0.6076851835E+00 0.1000000000E+01 0.6076851835E+00 0.2445650993E+03 0.9606049233E+00 0.8869460984E+00 0.9030760326E+00 0.9988806187E+00 0.1312287838E+03 -0.1219286188E+01 0.5302214850E+00 0.1000000000E+01 0.5302214850E+00 0.1307113219E+03 0.1865587854E+00 0.4935059263E+00 0.3541039475E+00 0.4250460356E+00 0.1080396373E+03 -0.7942589622E+00 -0.1335044267E+01 0.0000000000E+00 0.0000000000E+00 0.1092589235E+03 0.7487195604E+00 0.7296524509E+00 0.2687420347E+00 0.7473768233E+00 0.2110183677E+03 -0.3458029471E+00 0.1486879223E+00 0.1000000000E+01 0.1486879223E+00 0.2118126267E+03 0.1622697436E+00 0.2675811066E+00 0.2356591173E+00 0.7227852135E+00 -0.1091417754E+03 -0.2577612782E+01 -0.3189317600E+01 0.0000000000E+00 0.0000000000E+00 -0.1087959724E+03 0.8510832363E+00 0.1559521719E+00 0.8815292502E-01 0.5862107812E+00 0.1157701877E+03 -0.2307155721E+00 -0.5493164577E+00 0.0000000000E+00 0.0000000000E+00 0.1183478004E+03 0.4446003460E+00 0.3980157687E+00 0.4510250934E+00 0.3787454455E+00 0.2150579747E+03 0.4850916932E+00 0.1574783469E+00 0.1000000000E+01 0.1574783469E+00 0.2152886903E+03 0.5747027867E+00 0.2973686298E-01 0.7874561646E+00 0.7757584690E+00 0.6121387492E+02 -0.2049870281E+00 -0.3733274475E-01 0.0000000000E+00 0.0000000000E+00 0.6072878323E+02 0.1725886785E+00 0.6979202642E+00 0.9458805867E+00 0.4150214582E+00 0.4899658433E+02 -0.1537155974E+01 -0.1305744986E+01 0.0000000000E+00 0.0000000000E+00 0.4920157135E+02 0.2656475478E+00 0.7383356257E+00 0.2068617144E+00 0.7248333632E+00 0.2193152610E+03 -0.2573745277E+00 -0.3216962766E+00 0.0000000000E+00 0.0000000000E+00 0.2208524169E+03 0.2743353271E+00 0.7538424446E+00 0.8299668082E+00 0.2521454780E+00 0.2096880147E+03 -0.9395982108E-01 0.1633023688E-01 0.1000000000E+01 0.1633023688E-01 0.2099453893E+03 0.8090493553E+00 0.6925150617E+00 0.1006420352E+00 0.4906864404E+00 0.2739416442E+03 0.2922231357E+00 0.3851160284E+00 0.1000000000E+01 0.3851160284E+00 0.2740356040E+03 0.9670033776E+00 0.4257676615E+00 0.8770874687E+00 0.2090872317E+00 0.2682563153E+03 0.1101020742E+01 0.1579966482E+01 0.1000000000E+01 0.1579966482E+01 0.2679640922E+03 0.1291025453E+00 0.8264792725E+00 0.6371320882E+00 0.2790061771E+00 0.3460182438E+03 0.1054330127E+01 0.9260502105E+00 0.1000000000E+01 0.9260502105E+00 0.3449172230E+03 0.2568184161E+00 0.3471187665E+00 0.2510877931E-01 0.3253945615E-02 0.9833447185E+02 -0.7259405783E+00 -0.2093640671E+01 0.0000000000E+00 0.0000000000E+00 0.9728014172E+02 0.6890639537E+00 0.9786985954E-01 0.8987293341E+00 0.9439185499E+00 0.1201137602E+03 0.2757491737E+00 0.9053308710E+00 0.1000000000E+01 0.9053308710E+00 0.1208397008E+03 0.4390685284E+00 0.4247561111E+00 0.8759594508E+00 0.2504903465E+00 0.1882527650E+03 0.3532466211E+00 0.3435210579E+00 0.1000000000E+01 0.3435210579E+00 0.1879770158E+03 0.9912529164E+00 0.9877659245E+00 0.3818937340E+00 0.4879870354E+00 0.4177036950E+03 0.1242194763E+01 0.2091094374E+01 0.1000000000E+01 0.2091094374E+01 0.4173504484E+03 0.5981043004E+00 0.3389775469E+00 0.1956306943E+00 0.9650786608E+00 0.1423770294E+03 -0.2781279645E+00 -0.1803367621E+00 0.0000000000E+00 0.0000000000E+00 0.1411348346E+03 0.7705126936E-01 0.6840690039E-03 0.4971477487E+00 0.5622126102E+00 0.1351933858E+03 0.5285234946E+00 -0.3343808082E+00 0.0000000000E+00 0.0000000000E+00 0.1354715138E+03 0.1073388183E+00 0.4351842543E-01 0.4141762205E+00 0.5973861506E-01 -0.6240018059E+02 -0.1546663694E+01 -0.2921891615E+01 0.0000000000E+00 0.0000000000E+00 -0.6292870408E+02 0.2690334851E-01 0.1645783945E+00 0.6907633925E-01 0.9660337507E+00 0.1514528246E+03 0.1640822606E+00 -0.6093259064E+00 0.0000000000E+00 0.0000000000E+00 0.1529994883E+03 0.1292476677E+00 0.2655515983E+00 0.1257122406E+00 0.8456274359E+00 -0.3032374257E+02 -0.1835490497E+01 -0.2469351555E+01 0.0000000000E+00 0.0000000000E+00 -0.3048782483E+02 B34SRETURN$ B34SEEND$ /$ B34SEXEC OPTIONS OUTPUT=NORMAL$ B34SEEND$ B34SEXEC OPTIONS DUMPSET$ B34SEEND $ /$ New Command b34sexec reg$ model y= x1{0 to 4} x2{0 to 4} x3{0 to 3} x4 $ test x1{0 to 4}$ b34seend$ b34sexec reg$ model y= x1 x2 x3 x4 $ b34seend$ b34sexec robust$ model y= x1{0 to 4} x2{0 to 4} x3{0 to 3} x4 $ b34seend$ b34sexec robust$ model y= x1 x2 x3 x4 $ b34seend$ B34SEXEC REGRESSION RESIDUALA BLUS=BOTH $ MODEL Y = X1 X2 X3 X4 $ COMMENT=('TEST OF COMMENT IN NEW SYNTEXT') $ ORDER CONSTANT X4 $ RA VARS=(X1,X2) RESID=ALL PERIOD=36 GRAPH $ B34SEEND$ /$ .............. /$ IN CORE BLUS B34SEXEC REGRESSION RESIDUALA$ MODEL Y = X1 X2 X3 X4$ RA VARS=(X1,X2) RESID=ALLBLUS$ B34SEEND$ /$ B34SEXEC REGRESSION MAXGLS=2$ MODEL YGLS = X1 X2 X3 X4 $ B34SEEND$ /$ B34SEXEC LIST IEND=50 $ VAR X1 X2 X3 X4 Y PROBITY TOBITY YGLS$ B34SEEND$ /$ B34SEXEC PLOT GRAPH $ TITLE=('TEST OF THE NEW GRAPH OPTION')$ VAR X1 Y $ B34SEEND$ /$ B34SEXEC PLOT LPPLOT IBEGIN=10 IEND=30 CHARACTER=(Y X Z) $ VAR Y X1 X2 $ B34SEEND $ /$ B34SEXEC PROBIT $ MODEL PROBITY = X1 X2 X3 X4$ /$ B34SEXEC PROBIT NADJ=2$ MODEL PROBITY = X1 X2 X3 X4$ MEANS X1(.301,.444), X2(.8118,.77712) $ * THIS WILL TEST THE MEANS OPTION IN PROBIT $ B34SEEND $ /$ B34SEXEC TOBIT LOWER LIMIT=0.0 $ MODEL TOBITY = X1 X2 X3 X4$ B34SEEND$ /$ B34SEXEC LOGLIN KTEST KCHCK $ MODEL PROBITY = X1 X2 X3 X4 $ B34SEEND$ /$ B34SEXEC ECOMP REGFIRST BOTHP NREG=25 NPER=8$ MODEL Y = X1 X2 X3 X4$ B34SEEND $ /$ /$ DEBUB SEE PROBLEM .......... B34SEXEC AUTOC BOTH DIF=1 $ VAR YGLS X1$ B34SEEND$ /$ B34SEXEC RR IRB=1 IRRPT=PLOTO IBCPT=PLOTO IBCLS=LIST IOBSCH=50 ICUM=PLOTO ICUMSQ=PLOTO ICHOW=CALPRINT IQUANT=PLOT NTEST=2 IMORR=MTEST NBEGIN=20 MMRNOB=24 INCR=2$ MODEL Y = X1 X2 X3 X4 $ * TESTS RR OPTIONS $ B34SEEND $ B34SEXEC RR WHITE$ MODEL Y = X1 X2 X3 X4 $ B34SEEND $ B34SEXEC RR $ MODEL Y = X1 X2 X3 X4 $ BISPEC ITURNO IAUTO lm(10) df adf(10) adft(10) pp app(10) appt(10)$ B34SEEND$ /$ B34SEXEC QR PLOT IPCC=PCREG$ MODEL Y = X1 X2 X3 X4 $ TITLE=('THIS IS A TEST OF THE COMMENT FOR QR') $ B34SEEND$ /$ /$ Test data for B34SII commands B34SEXEC MPROBIT $ TITLE1=('SAMPLE DATA FROM SAS')$ TITLE2=(' TESTING TO SEE IF LIKE PROBIT COMMAND')$ VVALUES=(0.0 1.0)$ MODEL PROBITY = X1 X2 X3 X4$ /$ /$ B34SEXEC MLOGLIN IP=1 ICOV=PRINT BASIS=DCBASIS COV=ANAL$ LEVEL PROBITY(TEST) $ MODEL PROBITY =X1 X2 X3 X4$ B34SEEND $ /$B34SEXEC OPTIONS DEBUGSUBS=(SP04SA,SP04AA) $ B34SEEND $ B34SEXEC MLOGLIN IP=1 ICOV=PRINT BASIS=DCBASIS COV=ANAL VER=OLD $ * USING MORE OPTIONS $ VALUES X1(.301) X2(.8118)$ ELAST VALUES(1) EVARS(X3,X4)$ LEVEL PROBITY(TEST) $ MODEL PROBITY =X1 X2 X3 X4$ B34SEEND $ B34SEXEC MLOGLIN IP=1 ICOV=PRINT BASIS=DCBASIS COV=ANAL VER=NEW $ * USING MORE OPTIONS $ VALUES X1(.301) X2(.8118)$ ELAST VALUES(1) EVARS(X3,X4)$ LEVEL PROBITY(TEST) $ MODEL PROBITY =X1 X2 X3 X4$ SCALE X1(1.0) X2(1.0) $ B34SEEND $ /$ PROBLEM FROM KMENTA (1971) PAGE 565 - 582 B34SEXEC DATA NOOB=20 HEADING=('KMENTA TEST DATA')$ COMMENT=(' TEST DATA FROM PAGE 565 ') $ INPUT Q P D F A $ DATACARDS $ 98.485 100.323 87.4 98.0 1 99.187 104.264 97.6 99.1 2 102.163 103.435 96.7 99.1 3 101.504 104.506 98.2 98.1 4 104.240 98.001 99.8 110.8 5 103.243 99.456 100.5 108.2 6 103.993 101.066 103.2 105.6 7 99.900 104.763 107.8 109.8 8 100.350 96.446 96.6 108.7 9 102.820 91.228 88.9 100.6 10 95.435 93.085 75.1 81.0 11 92.424 98.801 76.9 68.6 12 94.535 102.908 84.6 70.9 13 98.757 98.756 90.6 81.4 14 105.797 95.119 103.1 102.3 15 100.225 98.451 105.1 105.0 16 103.522 86.498 96.4 110.5 17 99.929 104.016 104.4 92.5 18 105.223 105.769 110.7 89.3 19 106.232 113.490 127.1 93.0 20 B34SRETURN $ B34SEEND$ B34SEXEC LIST $ VAR Q P D F A $ B34SEEND $ /$ B34SEXEC SIMEQ PRINTSYS REDUCED OLS LIML LS2 LS3 KCOV=DIAG $ HEADING=('TEST CASE FROM KMENTA 1971 PAGE 565 - 582 ' ) $ EXOGENOUS CONSTANT D F A $ ENDOGENOUS P Q $ MODEL LVAR=Q RVAR=(CONSTANT P D) NAME=('DEMAND EQ.') $ MODEL LVAR=Q RVAR=(CONSTANT P F A) NAME=('SUPPLY EQ.') $ B34SEEND $ /$ B34SEXEC DATA NOOB =5 HEADING=('TEST DATA FROM LEE-JUDGE-ZELLNER') FILEF=FIXED $ FORMAT=('(F3.0,F5.0,F6.0)') $ INPUT VAR1 VAR2 SAMPSIZE $ PGMCARDS$ .50 .50 10.0 .75 .25 10.0 .88 .12 10.0 .94 .06 10.0 .97 .03 10.0 B34SRETURN$ B34SEEND$ B34SEXEC LIST $ VAR VAR1 VAR2 SAMPSIZE $ /$ /$ B34SEXEC TRANSPROB KDROP=VAR2 KW=ML KV=PRINT1 KP=PREDICT KEY1=TRANS SS=SAMPSIZE KEY3=PRINT KEY4=PRINTP KEY5=2 KEY7=RESA KEY8=QP $ VAR VAR1 VAR2 $ HEADING=('TEST DATA FROM TEXT ' ) $ PRIORI KALT=MBETA PARM=(1,1,99.0, 1.0 ) (1,2, 1.0,99.0 ) (2,1,50.0 50.0 ) (2,2,50.0 50.0 ) $ RECURSIVE NREC=6 NST=3 KP=PREDICTA KEY1=TRANS KEY4=PRINTP KEY5=2 KEY7=RESA KEY8=QP $ FINISH$ B34SEEND$ B34SEXEC SORT$ BY VAR2$ B34SEEND$ B34SEXEC LIST$ B34SEEND$ /$ DATA FROM BJ BOOK ************* B34SEXEC DATA NOOB=296$ INPUT TIME GASIN GASOUT$ * DATA COMES FROM B-J TIME SERIES BOOK $ DATACARDS$ 1 -0.109 53.8 2 0.000 53.6 3 0.178 53.5 4 0.339 53.5 5 0.373 53.4 6 0.441 53.1 7 0.461 52.7 8 0.348 52.4 9 0.127 52.2 10 -0.180 52.0 11 -0.588 52.0 12 -1.055 52.4 13 -1.421 53.0 14 -1.520 54.0 15 -1.302 54.9 16 -0.814 56.0 17 -0.475 56.8 18 -0.193 56.8 19 0.088 56.4 20 0.435 55.7 21 0.771 55.0 22 0.866 54.3 23 0.875 53.2 24 0.891 52.3 25 0.987 51.6 26 1.263 51.2 27 1.775 50.8 28 1.976 50.5 29 1.934 50.0 30 1.866 49.2 31 1.832 48.4 32 1.767 47.9 33 1.608 47.6 34 1.265 47.5 35 0.790 47.5 36 0.360 47.6 37 0.115 48.1 38 0.088 49.0 39 0.331 50.0 40 0.645 51.1 41 0.960 51.8 42 1.409 51.9 43 2.670 51.7 44 2.834 51.2 45 2.812 50.0 46 2.483 48.3 47 1.929 47.0 48 1.485 45.8 49 1.214 45.6 50 1.239 46.0 51 1.608 46.9 52 1.905 47.8 53 2.023 48.2 54 1.815 48.3 55 0.535 47.9 56 0.122 47.2 57 0.009 47.2 58 0.164 48.1 59 0.671 49.4 60 1.019 50.6 61 1.146 51.5 62 1.155 51.6 63 1.112 51.2 64 1.121 50.5 65 1.223 50.1 66 1.257 49.8 67 1.157 49.6 68 0.913 49.4 69 0.620 49.3 70 0.255 49.2 71 -0.280 49.3 72 -1.080 49.7 73 -1.551 50.3 74 -1.799 51.3 75 -1.825 52.8 76 -1.456 54.4 77 -0.944 56.0 78 -0.570 56.9 79 -0.431 57.5 80 -0.577 57.3 81 -0.960 56.6 82 -1.616 56.0 83 -1.875 55.4 84 -1.891 55.4 85 -1.746 56.4 86 -1.474 57.2 87 -1.201 58.0 88 -0.927 58.4 89 -0.524 58.4 90 0.040 58.1 91 0.788 57.7 92 0.943 57.0 93 0.930 56.0 94 1.006 54.7 95 1.137 53.2 96 1.198 52.1 97 1.054 51.6 98 0.595 51.0 99 -0.080 50.5 100 -0.314 50.4 101 -0.288 51.0 102 -0.153 51.8 103 -0.109 52.4 104 -0.187 53.0 105 -0.255 53.4 106 -0.229 53.6 107 -0.007 53.7 108 0.254 53.8 109 0.330 53.8 110 0.102 53.8 111 -0.423 53.3 112 -1.139 53.0 113 -2.275 52.9 114 -2.594 53.4 115 -2.716 54.6 116 -2.510 56.4 117 -1.790 58.0 118 -1.346 59.4 119 -1.081 60.2 120 -0.910 60.0 121 -0.876 59.4 122 -0.885 58.4 123 -0.800 57.6 124 -0.544 56.9 125 -0.416 56.4 126 -0.271 56.0 127 0.000 55.7 128 0.403 55.3 129 0.841 55.0 130 1.285 54.4 131 1.607 53.7 132 1.746 52.8 133 1.683 51.6 134 1.485 50.6 135 0.993 49.4 136 0.648 48.8 137 0.577 48.5 138 0.577 48.7 139 0.632 49.2 140 0.747 49.8 141 0.900 50.4 142 0.993 50.7 143 0.968 50.9 144 0.790 50.7 145 0.399 50.5 146 -0.161 50.4 147 -0.553 50.2 148 -0.603 50.4 149 -0.424 51.2 150 -0.194 52.3 151 -0.049 53.2 152 0.060 53.9 153 0.161 54.1 154 0.301 54.0 155 0.517 53.6 156 0.566 53.2 157 0.560 53.0 158 0.573 52.8 159 0.592 52.3 160 0.671 51.9 161 0.933 51.6 162 1.337 51.6 163 1.460 51.4 164 1.353 51.2 165 0.772 50.7 166 0.218 50.0 167 -0.237 49.4 168 -0.714 49.3 169 -1.099 49.7 170 -1.269 50.6 171 -1.175 51.8 172 -0.676 53.0 173 0.033 54.0 174 0.556 55.3 175 0.643 55.9 176 0.484 55.9 177 0.109 54.6 178 -0.310 53.5 179 -0.697 52.4 180 -1.047 52.1 181 -1.218 52.3 182 -1.183 53.0 183 -0.873 53.8 184 -0.336 54.6 185 0.063 55.4 186 0.084 55.9 187 0.000 55.9 188 0.001 55.2 189 0.209 54.4 190 0.556 53.7 191 0.782 53.6 192 0.858 53.6 193 0.918 53.2 194 0.862 52.5 195 0.416 52.0 196 -0.336 51.4 197 -0.959 51.0 198 -1.813 50.9 199 -2.378 52.4 200 -2.499 53.5 201 -2.473 55.6 202 -2.330 58.0 203 -2.053 59.5 204 -1.739 60.0 205 -1.261 60.4 206 -0.569 60.5 207 -0.137 60.2 208 -0.024 59.7 209 -0.050 59.0 210 -0.135 57.6 211 -0.276 56.4 212 -0.534 55.2 213 -0.871 54.5 214 -1.243 54.1 215 -1.439 54.1 216 -1.422 54.4 217 -1.175 55.5 218 -0.813 56.2 219 -0.634 57.0 220 -0.582 57.3 221 -0.625 57.4 222 -0.713 57.0 223 -0.848 56.4 224 -1.039 55.9 225 -1.346 55.5 226 -1.628 55.3 227 -1.619 55.2 228 -1.149 55.4 229 -0.488 56.0 230 -0.160 56.5 231 -0.007 57.1 232 -0.092 57.3 233 -0.620 56.8 234 -1.086 55.6 235 -1.525 55.0 236 -1.858 54.1 237 -2.029 54.3 238 -2.024 55.3 239 -1.961 56.4 240 -1.952 57.2 241 -1.794 57.8 242 -1.302 58.3 243 -1.030 58.6 244 -0.918 58.8 245 -0.798 58.8 246 -0.867 58.6 247 -1.047 58.0 248 -1.123 57.4 249 -0.876 57.0 250 -0.395 56.4 251 0.185 56.3 252 0.662 56.4 253 0.709 56.4 254 0.605 56.0 255 0.501 55.2 256 0.603 54.0 257 0.943 53.0 258 1.223 52.0 259 1.249 51.6 260 0.824 51.6 261 0.102 51.1 262 0.025 50.4 263 0.382 50.0 264 0.922 50.0 265 1.032 52.0 266 0.866 54.0 267 0.527 55.1 268 0.093 54.5 269 -0.458 52.8 270 -0.748 51.4 271 -0.947 50.8 272 -1.029 51.2 273 -0.928 52.0 274 -0.645 52.8 275 -0.424 53.8 276 -0.276 54.5 277 -0.158 54.9 278 -0.033 54.9 279 0.102 54.8 280 0.251 54.4 281 0.280 53.7 282 0.000 53.3 283 -0.493 52.8 284 -0.759 52.6 285 -0.824 52.6 286 -0.740 53.0 287 -0.528 54.3 288 -0.204 56.0 289 0.034 57.0 290 0.204 58.0 291 0.253 58.6 292 0.195 58.5 293 0.131 58.3 294 0.017 57.8 295 -0.182 57.3 296 -0.262 57.0 B34SRETURN$ B34SEEND $ b34sexec describe$ b34srun$ B34SEXEC BJIDEN TYPE=ARIMAI LIST=(GASIN,GASOUT) PLOT=(GASIN,GASOUT) NAC=24$ VAR GASIN GASOUT$ SERIESN VAR=GASIN NAME=('SERIES # 1 FROM BJ -- GAS INPUT')$ SERIESN VAR=GASOUT NAME=('SERIES # 2 FROM BJ -- GAS OUTPUT')$ RTRANS DIF=(2 1)(1 12)$ RAUTO GASIN GASOUT $ B34SEEND$ /$ B34SEXEC BJIDEN TYPE=TRANSI NCC=50 IRWCNT=5$ VAR GASIN GASOUT$ SERIESN VAR=GASIN NAME=('SERIES # 1 FROM BJ -- GAS INPUT')$ SERIESN VAR=GASOUT NAME=('SERIES # 2 FROM BJ -- GAS OUTPUT')$ MODELPRE VAR=GASIN P=(1,2,3) AR=(1.97494,-1.3732,.3424) AVEPA=-.061$ OUTPUT AVEPA=53.50799$ RCROSS=(GASIN,GASOUT) SPECTRAL$ PCROSS=(GASIN,GASOUT) SPECTRAL$ SPECTRAL NFREQ=50 PLOT=(ATRANSXY ATRANSYX) $ B34SEEND$ B34SEXEC BJIDEN TYPE=STARTVALUE$ STARTVALUE NTERMS=(1,2,3) DTERMS=(1,2,3) RUFILE=SDATA5 PRINT=(HA,RW)$ B34SEEND$ /$ B34SEXEC BJEST$ MODEL GASIN$ MODELN P=(1,2,3) AR=(.97,-.94,.5) AVEPA=.5$ FORECAST NF=24 NT=(296) $ /$ B34SEXEC BJEST $ MODEL GASOUT=GASIN$ MODELN P=(1,2) AR=(.8,.76) AVEPA=53.509 $ MODELI VAR=GASIN LAG=3 NUM=(0,1,2) DENOM=(1,2) NUMC=(-.53, .4 .6) DENOMC=(.6,.2)$ MODELPRE VAR=GASIN P=(1,2,3) AR=(1.97494,-1.3732,.3424) AVEPA=-.061 $ FORECAST NF=24 NT=(296) $ B34SEEND$ /$ B34SEXEC BTIDEN$ TITLE=('IDENTIFICATION RUN WITH GAS DATA') $ SERIESN VAR=GASIN NAME=('B-J GAS INPUT DATA') DIF=(1 1)(1 12) $ * DIF USED FOR TESTING ONLY ----------------------------$ SERIESN VAR=GASOUT NAME=('B-J GAS OUTPUT DATA') $ IDEN LAGRHO=36 ISACF ISCCF $ BISPEC IAUTO IMART NBD=1 ITURNO $ B34SEEND $ /$ B34SEXEC BTIDEN$ TITLE=('IDENTIFICATION RUN WITH GAS DATA') $ SERIESN VAR=GASIN NAME=('B-J GAS INPUT DATA') $ SERIESN VAR=GASOUT NAME=('B-J GAS OUTPUT DATA') $ ESTVAR P=6 OUTPUT=NORMAL NUMIRF=12 GRANGER$ BISPEC IAUTO IMART ITURNO NBD=1$ B34SEEND$ /$ B34SEXEC BTEST$ TITLE=('ESTIMATION RUN WITH GAS DATA') $ SERIESN VAR=GASIN NAME=('B-J GAS INPUT DATA') $ SERIESN VAR=GASOUT NAME=('B-J GAS OUTPUT DATA') $ AR(1,1,1)=1.9 $ AR(1,2,2)=1.5 $ AR(2,1,1)=-1.2 $ AR(2,2,2)=-.59 $ AR(3,1,1)=.17 $ AR(3,2,1)=-.44 $ AR(3,2,2)=-.17 $ AR(4,2,1)=.15 $ AR(4,2,2)=.13 $ AR(5,2,1)=-.12 $ AR(6,2,1)=.11 $ OUTPUT IPRINT LAGRHO=24 $ CONSTANT=(YES,NO) $ COVARIANCE(1,1)=.3$ COVARIANCE(2,1)= -.1 $ COVARIANCE(2,2)= .06 $ FORECAST NT=(296,250) NF=(24,20) SE ACTUAL OUTPUT=BOTH $ BISPEC INEWSP INAUTO IAUTO IMART ITURNO $ B34SEEND$ /$ B34SEXEC FORECAST K=2 MINVAR=3 MAXVAR=6 DMOD(4,5) ISEAR=2$ PGMCARDS$ TITLE('TEST RUN OF VAR SEARCH')$ SERIESN VAR=GASIN NAME('B-J GASIN SERIES')$ SERIESN VAR=GASOUT NAME('B-J GASOUT SERIES')$ OUTPUT LAGRHO=24$ FORECAST NT=250 NF=20 OUTPUT=NEW$ B34SRETURN$ B34SEEND$ /$ B34SEXEC VARFREQ DATAP=4 YEARS=(1900,1) YEARE=(1973,4) RUN=BATCH$ VAR TIME GASIN GASOUT CONSTANT$ VARF NLAGS=6 VAR=(GASIN,GASOUT) FEEDX=(GASIN) FEEDY=(GASOUT) DUMMY=CONS NREP=4 TABLE=('TEST OF PROGRAM') FREQ=(1.0 .9 .85 .8 .7 .6 .5 .4 .3 .2 .1 0.0)$ B34SEEND$ /$ B34SEXEC SPECTRAL LIST(ALL) PLOT(P S A K PH)$ WEIGHTS(1 2 3 4 3 2 1)$ VAR GASIN GASOUT$ B34SEEND$ B34SEXEC POLYSOLV$ POLYDIVIDE TOP(1.0) BOT(1.0 -.5) $ B34SEEND$ B34SEXEC POLYSOLV$ POLYDIVIDE TOP(1.0 -.4 0.0, -.3) BOT(1.0 -.5) NTERMS=40$ B34SEXEC POLYSOLV REAL IORDER=3$ PGMCARDS$ 0 1.0 1 -3. 2 4. 3 -2. B34SRETURN$ B34SEEND$ /$ B34SEXEC KFILTER $ SERIESN VAR=GASIN $ SERIESN VAR=GASOUT $ ESTIMATE K=8 NSTAR=5 FORECAST$ B34SEEND$ b34sexec data heading('Weighted Data fron Glass(1954)')$ input father son count income$ label father ='Father'$ label son ='Son '$ label count ='Count or weight'$ label income ='Mean income '$ datacards$ 1 1 50 361 1 2 45 512 1 3 8 634 1 4 18 684 1 5 8 718 2 1 28 431 2 2 174 631 2 3 84 652 2 4 154 793 2 5 55 815 3 1 11 492 3 2 78 670 3 3 110 732 3 4 223 831 3 5 96 842 4 1 14 496 4 2 150 683 4 3 185 763 4 4 714 852 4 5 447 852 5 1 3 582 5 2 42 696 5 3 72 762 5 4 320 864 5 5 411 875 b34sreturn$ b34seend$ b34sexec expand weight=count$ b34srun$ b34sexec describe$ b34srun$ b34sexec freq verthist horizhist$ table(father son)$ b34srun$ b34sexec data set$ gen if(father.eq.1)father=2$ gen if(son.eq.1)son=2$ gen if(father.eq.5)father=4$ gen if(son.eq.5)son=4$ b34srun$ b34sexec freq $ table(father son) $ b34srun$ b34sexec options open('LP.dmf') disp=unknown unit=60 form=formatted$ b34srun$ /$ Formatted DMF with files ASIM, BSIM and CSIM b34sexec options clean(60)$ b34srun$ b34sexec lpmax n=2 m1=4 m2=0$ * Sample problem from IMSL $ amatrix(1.0 0. 0. 1. 1. 1. -1. -1.)$ bvector(1. 1. 1.5 -.5)$ cvector(1 3.)$ b34srun$ b34sexec lpmax n=4 m1=3 m2=0$ * problem from Hadley(1962) page 3$ * max should be 127370588235294 $ * primal = 294.118 1500 0.0 58.8235 $ * Dual = 1.95353 .242353 1.37824 $ amatrix(1.5 1.0 2.4 1. 1. 5. 1. 3.5 1.5 3. 3.5 1.)$ bvector(2000 8000 5000)$ cvector(5.24 7.30 8.34 4.18)$ b34srun$ b34sexec lpmax n=2 m1=2 m2=0$ * problem from Hadley(1962) page 135$ * max should be 12.37. x1=1.053 x2=2.368$ amatrix(3. 5. 5. 2.)$ bvector(15. 10.)$ cvector(5. 3.)$ b34srun$ b34sexec lpmax n=4 m1=3 m2=0$ * problem from Hadley (1962) page 138$ * answer should be 40 x1=7.273 x3 = 6.36364$ amatrix(1 3 2 5 -2 -16 -1 -1 3 -1 -5 10)$ bvector(20 -4 -10)$ cvector(2 1 4 5)$ b34srun$ b34sexec lpmax n=2 m1=4 m2=0$ amatrix(1.0 0. 0. 1. 1. 1. -1. -1.)$ bvector(1. 1. 1.5 -.5)$ cvector(1 3.)$ simulate=a i=2 j=1 begin=0. end=3.5 nsteps=10 dmfunit=60 fdmf dmfname=asim rewinddmf$ b34srun$ b34sexec lpmax n=2 m1=4 m2=0$ amatrix(1.0 0. 0. 1. 1. 1. -1. -1.)$ bvector(1. 1. 1.5 -.5)$ cvector(1 3.)$ simulate=b i=2 begin=.5 end=7.5 nsteps=10 dmfunit=60 fdmf dmfname=bsim $ b34srun$ b34sexec lpmax n=2 m1=4 m2=0$ amatrix(1.0 0. 0. 1. 1. 1. -1. -1.)$ bvector(1. 1. 1.5 -.5)$ cvector(1 3.)$ simulate=c i=2 begin=.5 end=1.5 nsteps=10 dmfunit=60 fdmf dmfname=csim $ b34srun$ B34SEXEC DATA NOOB=50 HEADING('MARS TEST DATA') NOHEAD$ LABEL X1 = 'INPUT SERIES # 1 FROM FRIEDMAN'$ LABEL X2 = 'INPUT SERIES # 2 FROM FRIEDMAN'$ LABEL X3 = 'INPUT SERIES # 3 FROM FRIEDMAN'$ LABEL X4 = 'INPUT SERIES # 4 FROM FRIEDMAN'$ LABEL X5 = 'INPUT SERIES # 5 FROM FRIEDMAN'$ LABEL Y = 'OUTPUT SERIES FROM FRIEDMAN '$ INPUT X1 X2 X3 X4 X5 Y$ DATACARDS$ 0.500 0.333 0.200 0.143 0.091 9.653 0.250 0.667 0.400 0.286 0.182 8.765 0.750 0.111 0.600 0.429 0.273 8.292 0.125 0.444 0.800 0.571 0.364 13.796 0.625 0.778 0.040 0.714 0.455 24.134 0.375 0.222 0.240 0.857 0.545 17.111 0.875 0.556 0.440 0.020 0.636 13.420 0.062 0.889 0.640 0.163 0.727 7.951 0.562 0.037 0.840 0.306 0.818 9.939 0.312 0.370 0.080 0.449 0.909 15.615 0.812 0.704 0.280 0.592 0.008 16.956 0.188 0.148 0.480 0.735 0.099 8.924 0.688 0.481 0.680 0.878 0.190 19.789 0.438 0.815 0.880 0.041 0.281 14.596 0.938 0.259 0.120 0.184 0.372 13.713 0.031 0.593 0.320 0.327 0.463 6.862 0.531 0.926 0.520 0.469 0.554 17.263 0.281 0.074 0.720 0.612 0.645 11.168 0.781 0.407 0.920 0.755 0.736 21.977 0.156 0.741 0.160 0.898 0.826 18.631 0.656 0.185 0.360 0.061 0.917 7.909 0.406 0.519 0.560 0.204 0.017 9.866 0.906 0.852 0.760 0.347 0.107 10.142 0.094 0.296 0.960 0.490 0.198 12.161 0.594 0.630 0.008 0.633 0.289 22.449 0.344 0.963 0.208 0.776 0.380 19.795 0.844 0.012 0.408 0.918 0.471 11.399 0.219 0.346 0.608 0.082 0.562 5.924 0.719 0.679 0.808 0.224 0.653 14.938 0.469 0.123 0.048 0.367 0.744 13.153 0.969 0.457 0.248 0.510 0.835 21.515 0.016 0.790 0.448 0.653 0.926 11.851 0.516 0.235 0.648 0.796 0.025 12.167 0.266 0.568 0.848 0.939 0.116 17.204 0.766 0.901 0.088 0.102 0.207 14.095 0.141 0.049 0.288 0.245 0.298 3.925 0.641 0.383 0.488 0.388 0.388 14.185 0.391 0.716 0.688 0.531 0.479 14.717 0.891 0.160 0.888 0.673 0.570 18.291 0.078 0.494 0.128 0.816 0.661 15.741 0.578 0.827 0.328 0.959 0.752 24.748 0.328 0.272 0.528 0.122 0.843 11.031 0.828 0.605 0.728 0.265 0.934 20.199 0.203 0.938 0.928 0.408 0.033 13.252 0.703 0.086 0.168 0.551 0.124 10.041 0.453 0.420 0.368 0.694 0.215 15.079 0.953 0.753 0.568 0.837 0.306 18.463 0.047 0.198 0.768 0.980 0.397 15.470 0.547 0.531 0.968 0.003 0.488 16.365 0.297 0.864 0.016 0.146 0.579 15.694 B34SRETURN$ B34SEEND$ B34SEXEC REGRESSION RESIDUALA$ MODEL Y = X1 X2 X3 X4 X5$ B34SEEND$ B34SEXEC MARS NK=15 MI=2 M=LINEAR GRAPHRES IRES RABASIS FABASIS RABASIST$ MODEL Y = X1 X2 X3 X4 X5$ * BISPEC IAUTO ITURNO IMART$ * USES THE FIRST THREE OBS $ FORECAST X1(.5 .250 .750) X2(.333 .667 .111) X3(.20 .4 .6) X4(.143 .286 .429) X5(.091 .182 .273)$ B34SEEND$ /$ B34SEXEC MARS NK=15 MI=2 M=CUBIC GRAPHRES IRES RABASIS FABASIS$ MODEL Y = X1 X2 X3 X4 X5$ * BISPEC IAUTO ITURNO IMART$ * USES THE FIRST THREE OBS $ FORECAST X1(.5 .250 .750) X2(.333 .667 .111) X3(.20 .4 .6) X4(.143 .286 .429) X5(.091 .182 .273)$ B34SEEND$ /$ B34SEXEC MARS NK=15 MI=2 M=SEARCH PLOTRES IRES $ MODEL Y = X1 X2 X3 X4 X5$ * USES THE FIRST THREE OBS $ FORECAST X1(.5 .250 .750) X2(.333 .667 .111) X3(.20 .4 .6) X4(.143 .286 .429) X5(.091 .182 .273)$ B34SEEND$ /$ B34SEXEC MARS NK=15 MI=2 M=SEARCH IX=2 IRES$ MODEL Y = X1 X2 X3 X4 X5$ * USES THE FIRST THREE OBS $ FORECAST X1(.5 .250 .750) X2(.333 .667 .111) X3(.20 .4 .6) X4(.143 .286 .429) X5(.091 .182 .273)$ B34SEEND$ /$ B34SEXEC MARS NK=15 MI=2 M=SEARCH IX=10 IRES$ MODEL Y = X1 X2 X3 X4 X5$ * USES THE FIRST THREE OBS $ FORECAST X1(.5 .250 .750) X2(.333 .667 .111) X3(.20 .4 .6) X4(.143 .286 .429) X5(.091 .182 .273)$ B34SEEND$ /$ B34SEXEC MARS NK=15 MI=2 M=SEARCH IX=-2 IRES$ MODEL Y = X1 X2 X3 X4 X5$ * USES THE FIRST THREE OBS $ FORECAST X1(.5 .250 .750) X2(.333 .667 .111) X3(.20 .4 .6) X4(.143 .286 .429) X5(.091 .182 .273)$ B34SEEND$ B34SEXEC PISPLINE IRES CENTER=0.0$ MODEL Y = X1 X2 X3 X4 X5$ * USES THE FIRST THREE OBS $ * BISPEC IAUTO ITURNO IMART$ FORECAST X1(.5 .250 .750) X2(.333 .667 .111) X3(.20 .4 .6) X4(.143 .286 .429) X5(.091 .182 .273)$ B34SEEND$ B34SEXEC DATA NOOB=88 HEADING('PIMPLE DATA') NOHEAD$ INPUT Y E_RATIO C_RATIO$ DATACARDS$ 3.741 .907 12. 2.295 .761 12. 1.498 1.108 12. 2.881 1.016 12. .76 1.189 12. 3.12 1.001 9. .638 1.231 9. 1.17 1.123 9. 2.358 1.042 12. .606 1.215 12. 3.669 .93 12. 1. 1.152 12. .981 1.138 15. 1.192 .601 18. .926 .696 7.5 1.59 .686 12. 1.806 1.072 12. 1.962 1.074 15. 4.028 .934 15. 3.148 .808 9. 1.836 1.071 9. 2.845 1.009 7.5 1.013 1.142 7.5 .414 1.229 18. .812 1.175 18. .374 .568 15. 3.623 .977 15. 1.869 .767 7.5 2.836 1.006 7.5 3.567 .893 9. .866 1.152 15. 1.369 .693 15. .542 1.232 15. 2.739 1.036 15. 1.2 1.125 15. 1.719 1.081 9. 3.423 .868 9. 1.634 .762 7.5 1.021 1.144 7.5 2.157 1.045 7.5 3.361 .797 18. 1.39 1.115 18. 1.947 1.07 18. .962 1.219 18. .571 .637 9. 2.219 .733 9. 1.419 .715 9. 3.519 .872 9. 1.732 .765 7.5 3.206 .878 7.5 2.471 .811 7.5 1.777 .676 15. 2.571 1.045 18. 3.952 .968 18. 3.931 .846 15. 1.587 .684 15. 1.397 .729 7.5 3.536 .911 7.5 2.202 .808 7.5 .756 1.168 7.5 1.62 .749 7.5 3.656 .892 7.5 2.964 1.002 7.5 3.76 .812 18. .672 1.23 18. 3.677 .804 18. 3.517 .813 12. 3.29 1.002 12. 1.139 .696 9. .727 1.199 9. 2.581 1.03 9. .923 .602 15. 1.527 .694 15. 3.388 .816 15. 2.085 1.037 15. .966 1.181 15. 3.488 .899 7.5 .754 1.227 7.5 .797 1.18 9. 2.064 .795 7.5 3.732 .99 18. .586 1.201 18. .561 .629 7.5 .563 .608 9. .678 .584 12. .37 .562 15. .53 .535 18. 1.9 .655 18. B34SRETURN$ B34SEEND$ B34SEXEC PISPLINE CENTER=2.526 IRES PLOTRES PMODEL$ BISPEC IAUTO ITURNO$ FORECAST C_RATIO(12. 12.) E_RATIO(.907 .761)$ MODEL Y = E_RATIO C_RATIO$ B34SEEND$ B34SEXEC PISPLINE CENTER=2.526 IRES NG=200$ BISPEC IAUTO ITURNO$ FORECAST E_RATIO(.907 .761 .5) C_RATIO(12. 12. 50.)$ MODEL Y = E_RATIO C_RATIO$ B34SEEND$ B34SEXEC PISPLINE CENTER=2.526 IRES GRAPHRES$ BISPEC IAUTO ITURNO$ FORECAST C_RATIO(12. 12.) E_RATIO(.907 .761) $ MODEL Y = E_RATIO C_RATIO$ B34SEEND$ B34SEXEC MARS NK=20 MI=3 IRES$ FORECAST C_RATIO(12. 12. .5) E_RATIO(.907 .761 50.)$ MODEL Y = E_RATIO C_RATIO$ B34SEEND$ b34sexec matrix; * Math with matrix and vectors ; * For bigger problems, change n; n=3; right=integers(1,((n*n)-1))+10; call print('Right ',right); x=matrix(n,n:right,-7); x2=x*2.; v=vector(n:integers(1,n)); call print('X, 2*x v' ,x,x2,v) ; call print('Inverse of x (INV)' ,(1./x)) ; call print('X*inv' ,x*(1./x)); call print('Vector times matrix (v*x)' ,v*x) ; call print('Matrix times vector (x*v)' ,x*v) ; call print('Matrix times matrix (x*x2)' ,x*x2) ; call print('Matrix times scaler (x*2.)' ,x*2.) ; call print('Scaler times Matrix (3.*x)' ,3.*x) ; call print('Vector plus matrix (v+x)' ,v+x) ; call print('Matrix plus vector (x+v)' ,x+v) ; call print('Matrix plus matrix (x+x2)' ,x+x2) ; call print('Matrix plus scaler (x+2.)' ,x+2.) ; call print('Scaler plus matrix (3.+x)' ,3.+x) ; call print('Vector minus matrix (v-x)' ,v-x) ; call print('Matrix minus vector (x-v)' ,x-v) ; call print('Matrix minus matrix (x-x2)' ,x-x2) ; call print('Matrix minus scaler (x-2.)' ,x-2.) ; call print('Scaler minus matrix (3.-x)' ,3.-x) ; call print('Array Math Here '); x=afam(x); v=afam(v); x2=x*2.; call print('X, 2*x v' ,x,x2,v) ; call print('Inverse of x (INV)' ,(1./x)) ; call print('X*inv' ,x*(1./x)); call print('Array(1) times Array(2) (v*x)' ,v*x) ; call print('Array(2) times Array(1) (x*v)' ,x*v) ; call print('Array(2) times Array(2) (x*x2)' ,x*x2) ; call print('Array(2) times Scaler (x*2.)' ,x*2.) ; call print('Scaler times Array(2) (3.*x)' ,3.*x) ; call print('Array(1) plus Array(2) (v+x)' ,v+x) ; call print('Array(2) plus Array(1) (x+v)' ,x+v) ; call print('Array(2) plus Array(2) (x+x2)' ,x+x2) ; call print('Array(2) plus Scaler (x+2.)' ,x+2.) ; call print('Scaler plus Array(2) (3.+x)' ,3.+x) ; call print('Array(1) minus Array(2) (v-x)' ,v-x) ; call print('Array(2) minus Array(1) (x-v)' ,x-v) ; call print('Array(2) minus Array(2) (x-x2)' ,x-x2) ; call print('Array(2) minus scaler (x-2.)' ,x-2.) ; call print('Scaler minus Array(2) (3.-x)' ,3.-x) ; call print(' Complex Results ' '++++++++++++++++++++++++++++++++++++++++'); x=mfam(complex(x,x2)); v=vfam(complex(v,v+8.0)); x2=mfam(complex(x2)); call print('X, x2 v' ,x,x2,v) ; call print('Inverse of x (INV)' , (complex(1.)/x)) ; call print('X*inv' , x*(complex(1.)/x)); call print('Vector times matrix (v*x)' ,v*x) ; call print('Matrix times vector (x*v)' ,x*v) ; call print('Matrix times matrix (x*x2)' ,x*x2) ; call print('Matrix times scaler (x*2.)',x*complex(2.)) ; call print('Scaler times Matrix (3.*x)',complex(3.)*x) ; call print('Vector plus matrix (v+x)',v+x) ; call print('Matrix plus vector (x+v)',x+v) ; call print('Matrix plus matrix (x+x2)',x+x2) ; call print('Matrix plus scaler (x+2.)',x+complex(2.)) ; call print('Scaler plus matrix (3.+x)',complex(3.)+x) ; call print('Vector minus matrix (v-x)',,v-x) ; call print('Matrix minus vector (x-v)' ,x-v) ; call print('Matrix minus matrix (x-x2)',x-x2) ; call print('Matrix minus scaler (x-2.)',x-complex(2.)) ; call print('Scaler minus matrix (3.-x)',complex(3.)-x) ; call print('Array Math Here '); x=afam(x); v=afam(v); x2=afam(x2); call print('X, 2*x v' ,x,x2,v) ; call print('Inverse of x (INV)', (complex(1.)/x)) ; call print('X*inv' , x*(complex(1.)/x)); call print('Array(1) times Array(2) (v*x)' ,v*x) ; call print('Array(2) times Array(1) (x*v)' ,x*v) ; call print('Array(2) times Array(2) (x*x2)' ,x*x2) ; call print('Array(2) times Scaler (x*complex(2.))',x*complex(2.)); call print('Scaler times Array(2) (complex(3.)*x)',complex(3.)*x); call print('Array(1) plus Array(2) (v+x)' ,v+x) ; call print('Array(2) plus Array(1) (x+v)' ,x+v) ; call print('Array(2) plus Array(2) (x+x2)' ,x+x2) ; call print('Array(2) plus Scaler (x+complex(2.))',x+complex(2.)); call print('Scaler plus Array(2) (complex(3.)+x)',complex(3.)+x); call print('Array(1) minus Array(2) (v-x)' ,v-x) ; call print('Array(2) minus Array(1) (x-v)' ,x-v) ; call print('Array(2) minus Array(2) (x-x2)' ,x-x2) ; call print('Array(2) minus scaler (x-complex(2.))',x-complex(2.)); call print('Scaler minus Array(2) (complex(3.)-x)',complex(3.)-x); b34srun; b34sexec matrix; * Test case for Real Matrix from IMSL Math (10) pp 295-297; * Eigenvectors have NOT been normalized but are tested below; a=matrix(3,3:8.,-1.,-5.,-4., 4.,-2.,18.,-5.,-7.); call print('A Matrix',a); e=eigenval(a); call print('Eigenvalues of a', e, 'Sum of the eigenvalues of General Martix A',sum(e), 'Trace of General Matrix A',trace(a), 'Product of the eigenvalues of Martix A',prod(e), 'Determinant of Matrix A',det(a)); ee=eigenval(a,evec:); call print('Non scaled Eigenvectors',evec); ee=eigenval(a,evec); call print('Scaled Eigenvectors',evec); call print('Test transpose(evec)*evec ', transpose(evec)*evec , ' ' 'Test evec*transpose(evec) ', evec*transpose(evec)) ; * Complex Case See IMSL Math (10) pp 302-304 ; r=matrix(4,4:5., 5.,-6.,-7., 3., 6.,-5.,-6., 2., 3.,-1.,-5., 1., 2.,-3.,0.0); i=matrix(4,4:9., 5.,-6.,-7., 3.,10.,-5.,-6., 2., 3., 3.,-5., 1., 2.,-3., 4.); ca=complex(r,i); call print('CA Complex Matrix',ca); ce=eigenval(ca,cevec:); call print('Non scaled Eigenvectors of CA',cevec); ee=eigenval(ca,cevec); call print('Scaled Eigenvectors of CA',cevec); call print('Eigenvalues of ca', ce, 'Sum of the eigenvalues of General Martix CA',sum(ce), 'Trace of General Matrix CA',trace(ca), 'Product of the eigenvalues of Martix CA',prod(ce), 'Determinant of Matrix CA',det(ca)); b34srun; b34sexec matrix; * Example from Limdep 7.0 Manual page 376 ; * Eigenalysis of Klein Model 1 ; r=matrix(3,3:.172,-.051,-.008,1.511,.848,.743,-.287,-.161,.818); call print(r,eigenval(r)); b34srun; b34sexec matrix; * Test Problem of FFT from MATLAB page 6-32 ; x=array(8:4., 3., 7., -9., 1., 0., 0., 0.); call print(x,fft(x)); b34srun;