phil2

30 Reputation

4 Badges

7 years, 13 days
Maplesoft
Technical Support Analyst
Waterloo, Ontario, Canada

MaplePrimes Activity


These are answers submitted by phil2

pdsolve([PDE1, PDE2], [w1(x, t), w2(x, t)]) generates some output, but I'm not sure if it's what you're looking for; you may need to give Maple more specific input (maybe some ICs).

The message says the independent variable in the equations does not match the 'continuoustimevar' (t) specified by the DynamicSystems[SystemOptions], so double-check the terms in your equation and make sure that you used the same independant variable as what your command is expecting you to use.  If you need further assistance, please post your worksheet.

I assume that method = bvp[middefer] is used to address the singularity that Maple finds when no method is specified, so the error message Error, (in dsolve/numeric/bvp) initial Newton iteration is not converging suggests that the problem was either set up incorrectly to begin with, or it cannot be solved (at least not with this method/approach).

A couple of suggestions:

  1. In the first DE, there is a fourth derivative for f(eta), but the ICs only go up to the second derivative. Try adding an IC for the third derivative; the first derivative has two (which may explain the singularity), so perhaps one of them was supposed to be for the third.
  2. Since the first two DEs do not use phi(eta), you may want to concentrate on those first, then use the solution for f(eta) to find phi(eta) in a separate call to dsolve().

In general, this error message occurs when either the worksheet or the Maple installation itself has been corrupted, or when Maple senses that some loop will not terminate, some computation would be too demanding (for example 2^(2^32-1)), or Maple was consuming too much of the RAM or CPU during a calculation.

Considering that the code works elsewhere, the worksheet itself is not likely the issue (unless you copied the code rather than transfering the file).  So, if other code works fine on that installation of Maple, the 100x100 Matrix is likely too demanding for the machine you experienced the issue on. This could be a hardware limitation, but you might also have been using 32-bit Maple (which has limited memory).

Hi Lilian,The source of the issue for list_1 appears to be rooted in the feasibilitytolerance option.  Without any constraints, there is no need for the option (only optimalitytolerance would be relevant), so when Maple executes the call to LSSolve, Maple's algorithms seem to find the minimum without any issue.  However, by adding constraints at zero, Maple's algorithms seem to change, now looking at values even closer to 0 and (for some reason) ignoring the ones further away.

 

 

As you may have guessed by the mention of the feasibilitytoleranceoption, Maple's algorithms are designed to violate constraints in order to achieve a better solution.  However, they will only do so if the violations are within the set feasibility tolerance. The LSSolve help page does not seem to specify what the default feasibility tolerance actually is, but you will notice that the negative value that you saw was a very small number (i.e. reasonably close zero).  If you incorporate a suitably small tolerance, the negatives go away (see attached).

Admittedly, this still doesn't quite explain why the new minimum is actually larger.  However, considering that, as you pointed out, the constraints were redundant to begin with, the notion of feasibility tolerance makes those constraints inadvisable anyway.

For list_2, the initial call to LSSolve did yeild a negative value, so incorporating a constraint to avoid it is perfectly reasonable.  However, based on what we now know about Maple's feasibility tolerance in problems like these, only one constraint should be used.  Once it is added, another negative does come up, but after a constraint is added for that one, no others appear (see attached).

On a side note, you may want to consider upgrading to 2016.2.

list_1 := [-0.777253031437943e-1+6.57999999999998*D2+(276.680000000000-Phi12_17)*D17, .978053142243381-633.614000000000*D2+(-6.58400000000000-Phi23_17)*D17, -0.178064883733416e-1+253.738100000000*D2+(-5.34190000000000-Phi34_17)*D17, -0.353869747788979e-2+39.4444000000000*D2+(-1.06160000000000-Phi45_17)*D17, -0.884707703096015e-3+9.86159000000000*D2+(-.265410000000000-Phi56_17)*D17, -0.226985309680582e-1+190.480500000000*D2+(-6.80950000000000-Phi67_17)*D17, -0.862540843423853e-3+2.15544000000000*D2+(-.258760000000000-Phi78_17)*D17, -0.141061228211511e-2+2.10972000000000*D2+(-.423180000000000-Phi89_17)*D17, -0.163814759658455e-2+2.45006000000000*D2+(-.491440000000000-Phi910_17)*D17, -0.525504575536290e-2+6.99450000000000*D2+(-1.57650000000000-Phi1011_17)*D17, -0.305389325672901e-1+28.4964000000000*D2+(-9.16160000000000-Phi1112_17)*D17, -0.987175261941468e-1+25.7100000000000*D2+(-29.6150000000000-Phi1213_17)*D17, -0.493604297782517e-2+1.28550000000000*D2+(-1.48080000000000-Phi1314_17)*D17, -0.522437882168112e-1+6.41700000000000*D2+(-15.6730000000000-Phi1415_17)*D17, -.141614566362726+13.9970000000000*D2+(-42.4840000000000-Phi1516_17)*D17, 0.777253031360147e-1-42.4900000000000*D15+(155.450000000000-Phi12_17)*D28, 0.219470855957043e-1-13.9960000000000*D15+(43.8940000000000-Phi23_17)*D28, 0.178065694473007e-1-12.8320000000000*D15+(35.6130000000000-Phi34_17)*D28, 0.353886380190268e-2-2.92230000000000*D15+(7.07770000000000-Phi45_17)*D28, 0.884703450426918e-3-.730700000000000*D15+(1.76940000000000-Phi56_17)*D28, 0.226985885266366e-1-20.8670000000000*D15+(45.3970000000000-Phi67_17)*D28, 0.862553364039491e-3-.965300000000000*D15+(1.72510000000000-Phi78_17)*D28, 0.141060550149453e-2-1.62670000000000*D15+(2.82120000000000-Phi89_17)*D28, 0.163815638896446e-2-1.88900000000000*D15+(3.27630000000000-Phi910_17)*D28, 0.525502049507568e-2-6.09000000000000*D15+(10.5100000000000-Phi1011_17)*D28, 0.305386191034954e-1-35.8000000000000*D15+(61.0770000000000-Phi1112_17)*D28, 0.987153849993141e-1-117.980000000000*D15+(197.430000000000-Phi1213_17)*D28, 0.493586925035572e-2-5.89830000000000*D15+(9.87170000000000-Phi1314_17)*D28, 0.522452037612234e-1-62.6800000000000*D15+(104.490000000000-Phi1415_17)*D28, -.858384447685986+609.990000000000*D15+(283.230000000000-Phi1516_17)*D28, 0.133722502645145e-2-5.34000000000003*D3+(264.760000000000-Phi12_19)*D19, 0.492898645263975e-1-253.740000000000*D3+(373.290000000000-Phi23_19)*D19, -.792605146598487+1386.78000000000*D3+(1127.70000000000-Phi34_19)*D19, .132715839620078-165.954000000000*D3+(-206.460000000000-Phi45_19)*D19, 0.331294375945810e-1-41.4870000000000*D3+(-51.6140000000000-Phi56_19)*D19, .512857041219517-749.640000000000*D3+(-946.930000000000-Phi67_19)*D19, 0.345803892354315e-2-6.58880000000000*D3+(-9.00300000000000-Phi78_19)*D19, 0.253797811562394e-2-5.18510000000000*D3+(-7.71800000000000-Phi89_19)*D19, 0.294701900455624e-2-6.02140000000000*D3+(-8.96290000000000-Phi910_19)*D19, 0.754969755588153e-2-15.8980000000000*D3+(-24.4690000000000-Phi1011_19)*D19, 0.181312716670633e-1-44.7810000000000*D3+(-82.4390000000000-Phi1112_19)*D19, 0.133647518064223e-1-35.1810000000000*D3+(-90.5060000000000-Phi1213_19)*D19, 0.667987641718040e-3-1.75900000000000*D3+(-4.52530000000000-Phi1314_19)*D19, 0.217080361770671e-2-6.75400000000000*D3+(-28.8440000000000-Phi1415_19)*D19, 0.365299883394169e-2-12.8320000000000*D3+(-69.3130000000000-Phi1516_19)*D19, -0.130020251659236e-2+6.81000000000000*D6+(263.430000000000-Phi12_19)*D20, -0.493576878414119e-1+190.470000000000*D6+(323.980000000000-Phi23_19)*D20, -.207482316974689+749.650000000000*D6+(920.260000000000-Phi34_19)*D20, -.132570648903326+410.334000000000*D6+(461.-Phi45_19)*D20, -0.329843798422150e-1+102.583000000000*D6+(115.250000000000-Phi56_19)*D20, .487420161727992-1880.19000000000*D6+(-1459.80000000000-Phi67_19)*D20, -0.345553822678983e-2+26.0430000000000*D6+(-12.4610000000000-Phi78_19)*D20, -0.254139584204717e-2+18.7200000000000*D6+(-10.2570000000000-Phi89_19)*D20, -0.295045955688278e-2+21.7400000000000*D6+(-11.9110000000000-Phi910_19)*D20, -0.755117615405593e-2+55.1470000000000*D6+(-32.0200000000000-Phi1011_19)*D20, -0.181178219908242e-1+124.980000000000*D6+(-100.570000000000-Phi1112_19)*D20, -0.133520796896221e-1+89.6200000000000*D6+(-103.870000000000-Phi1213_19)*D20, -0.668104062372097e-3+4.48110000000000*D6+(-5.19340000000000-Phi1314_19)*D20, -0.217033804692734e-2+13.3910000000000*D6+(-31.0150000000000-Phi1415_19)*D20, -0.365556938318863e-2+20.8670000000000*D6+(-72.9670000000000-Phi1516_19)*D20, 0.135489886901528e-2-.259999999999991*D7+(256.360000000000-Phi12_21)*D22, 0.680675384195741e-2-2.16000000000000*D7+(131.350000000000-Phi23_21)*D22, 0.167104193845212e-1-6.59000000000000*D7+(164.020000000000-Phi34_21)*D22, 0.655190381659506e-2-2.76600000000000*D7+(47.9000000000000-Phi45_21)*D22, 0.163878244157080e-2-.692000000000000*D7+(11.9750000000000-Phi56_21)*D22, 0.603897781618217e-1-26.0400000000000*D7+(394.350000000000-Phi67_21)*D22, -.981708865268294+89.7700000000000*D7+(51.2660000000000-Phi78_21)*D22, .286945758784479-5.67000000000000*D7+(-34.6470000000000-Phi89_21)*D22, .333288992029305-6.58400000000000*D7+(-40.2350000000000-Phi910_21)*D22, .115166403866294-14.3630000000000*D7+(-101.530000000000-Phi1011_21)*D22, 0.796486977999666e-1-12.5100000000000*D7+(-238.060000000000-Phi1112_21)*D22, 0.513893785319347e-1-8.19999999999999*D7+(-201.690000000000-Phi1213_21)*D22, 0.257108190144080e-2-.409500000000000*D7+(-10.0840000000000-Phi1314_21)*D22, 0.486795807939043e-2-.850999999999999*D7+(-45.2570000000000-Phi1415_21)*D22, 0.524862014258993e-2-.965000000000003*D7+(-94.7990000000000-Phi1516_21)*D22, -0.136101769174909e-2+1.57999999999998*D10+(255.450000000000-Phi12_21)*D23, -0.680508845874531e-2+6.99000000000001*D10+(126.790000000000-Phi23_21)*D23, -0.167210744986313e-1+15.9000000000000*D10+(152.820000000000-Phi34_21)*D23, -0.654955044315159e-2+6.04100000000000*D10+(43.5110000000000-Phi45_21)*D23, -0.163599881722489e-2+1.51050000000000*D10+(10.8780000000000-Phi56_21)*D23, -0.603847441196420e-1+55.1500000000000*D10+(353.890000000000-Phi67_21)*D23, -0.182876336377057e-1+14.3600000000000*D10+(39.0120000000000-Phi78_21)*D23, -.286946970814718+35.7040000000000*D10+(83.0820000000000-Phi89_21)*D23, -.333290206967530+41.4630000000000*D10+(96.4820000000000-Phi910_21)*D23, .884841237755273-359.510000000000*D10+(-178.690000000000-Phi1011_21)*D23, -0.796889746454705e-1+94.6100000000000*D10+(-291.440000000000-Phi1112_21)*D23, -0.513575859600819e-1+60.8700000000000*D10+(-236.110000000000-Phi1213_21)*D23, -0.256926809156710e-2+3.04300000000000*D10+(-11.8060000000000-Phi1314_21)*D23, -0.486633264478440e-2+5.69300000000000*D10+(-48.5180000000000-Phi1415_21)*D23, -0.524686208104891e-2+6.08500000000001*D10+(-98.3150000000000-Phi1516_21)*D23, 0.740538677273994e-2-9.16000000000000*D11+(244.710000000000-Phi12_23)*D25, 0.642783571664902e-2-28.4970000000000*D11+(91.3030000000000-Phi23_23)*D25, 0.879545937424008e-2-44.7820000000000*D11+(92.1380000000000-Phi34_23)*D25, 0.267263958813606e-2-14.5080000000000*D11+(22.9620000000000-Phi45_23)*D25, 0.668159897034014e-3-3.62700000000000*D11+(5.74050000000000-Phi56_23)*D25, 0.224071702943972e-1-124.980000000000*D11+(173.760000000000-Phi67_23)*D25, 0.204975705587393e-2-12.5020000000000*D11+(12.1500000000000-Phi78_23)*D25, 0.398190796934755e-2-24.6980000000000*D11+(22.6800000000000-Phi89_23)*D25, 0.462414151228524e-2-28.6810000000000*D11+(26.3380000000000-Phi910_23)*D25, 0.152190448689061e-1-94.6100000000000*D11+(86.2100000000000-Phi1011_23)*D25, -.905930086863230+912.110000000000*D11+(526.060000000000-Phi1112_23)*D25, .687841196392981-376.270000000000*D11+(-673.250000000000-Phi1213_23)*D25, 0.343937963386998e-1-18.8130000000000*D11+(-33.6620000000000-Phi1314_23)*D25, 0.305873475333330e-1-34.1060000000000*D11+(-88.3170000000000-Phi1415_23)*D25, 0.294965405620479e-1-35.8000000000000*D11+(-140.200000000000-Phi1516_23)*D25, -0.740003311264899e-2+15.6700000000000*D14+(213.610000000000-Phi12_23)*D26, -0.643002877220721e-2+6.41700000000000*D14+(64.3070000000000-Phi23_23)*D26, -0.879503935483085e-2+6.75400000000000*D14+(55.1990000000000-Phi34_23)*D26, -0.267001194740175e-2+1.73800000000000*D14+(11.7380000000000-Phi45_23)*D26, -0.668002989087778e-3+.434300000000000*D14+(2.93440000000000-Phi56_23)*D26, -0.224051002552570e-1+13.3910000000000*D14+(79.6550000000000-Phi67_23)*D26, -0.205000917309872e-2+.851000000000000*D14+(3.54140000000000-Phi78_23)*D26, -0.398201781818492e-2+1.50890000000000*D14+(5.95680000000000-Phi89_23)*D26, -0.462402069093096e-2+1.75230000000000*D14+(6.91760000000000-Phi910_23)*D26, -0.152200681046646e-1+5.69300000000000*D14+(22.2930000000000-Phi1011_23)*D26, -0.941004210676048e-1+34.1030000000000*D14+(130.980000000000-Phi1112_23)*D26, -.687853077910222+122.340000000000*D14+(437.750000000000-Phi1213_23)*D26, -0.344056792915157e-1+6.11800000000000*D14+(21.8880000000000-Phi1314_23)*D26, .969399863074721-383.950000000000*D14+(-216.780000000000-Phi1415_23)*D26, -0.295001320031278e-1+62.6800000000000*D14+(-264.080000000000-Phi1516_23)*D26]:

with(Optimization):

[0.182130325886275e-7, [D10 = HFloat(0.002000093347416389), D11 = HFloat(6.666205093053039e-4), D14 = HFloat(0.0022221520824599844), D15 = HFloat(0.0012820279138359672), D17 = HFloat(0.004998861403449358), D19 = HFloat(3.029255266744349e-4), D2 = HFloat(0.0010000234934204896), D20 = HFloat(1.4284944659735357e-4), D22 = HFloat(0.01111211271354064), D23 = HFloat(0.0022222805411430922), D25 = HFloat(7.142936214972088e-4), D26 = HFloat(8.333263499191028e-4), D28 = HFloat(2.1739653171991575e-4), D3 = HFloat(4.002175679009348e-4), D6 = HFloat(1.6687886220206148e-4), D7 = HFloat(0.009999698285428371), Phi1011_17 = HFloat(-1.228502431180965), Phi1011_19 = HFloat(-20.533519373654585), Phi1011_21 = HFloat(-104.09096431277415), Phi1011_23 = HFloat(19.214449939212358), Phi1112_17 = HFloat(-9.57006175328916), Phi1112_19 = HFloat(-81.68486302365807), Phi1112_21 = HFloat(-242.1498491748034), Phi1112_23 = HFloat(109.00135135008516), Phi1213_17 = HFloat(-44.219748932875795), Phi1213_19 = HFloat(-92.82671959308529), Phi1213_21 = HFloat(-204.44416589041828), Phi1213_23 = HFloat(-61.44476127504735), Phi12_17 = HFloat(262.4476511854361), Phi12_19 = HFloat(262.14919240665415), Phi12_21 = HFloat(256.2484052767522), Phi12_23 = HFloat(246.52117286304335), Phi1314_17 = HFloat(-2.2110615918866032), Phi1314_19 = HFloat(-4.64255435896982), Phi1314_21 = HFloat(-10.221215875683743), Phi1314_23 = HFloat(-3.077984004772912), Phi1415_17 = HFloat(-24.84038315298048), Phi1415_19 = HFloat(-30.594450771890468), Phi1415_21 = HFloat(-45.58470252594737), Phi1415_23 = HFloat(-77.32976800445691), Phi1516_17 = HFloat(-68.01327120138), Phi1516_19 = HFloat(-74.2023324472602), Phi1516_21 = HFloat(-95.1952296374022), Phi1516_23 = HFloat(-132.32846708126223), Phi23_17 = HFloat(62.31653108816923), Phi23_19 = HFloat(200.80422545113467), Phi23_21 = HFloat(130.01879159887292), Phi23_23 = HFloat(73.70432624302144), Phi34_17 = HFloat(41.85634777007731), Phi34_19 = HFloat(343.4099879320013), Phi34_21 = HFloat(159.59399606138695), Phi34_23 = HFloat(62.656475770029026), Phi45_17 = HFloat(6.121364130365998), Phi45_19 = HFloat(12.383917194570708), Phi45_21 = HFloat(46.000528179934825), Phi45_23 = HFloat(13.166579651635805), Phi56_17 = HFloat(1.530420688109261), Phi56_19 = HFloat(3.166146874145713), Phi56_21 = HFloat(11.499811489254649), Phi56_23 = HFloat(3.2909339469103527), Phi67_17 = HFloat(26.75542459204723), Phi67_19 = HFloat(-244.28897794294045), Phi67_21 = HFloat(376.3514935402879), Phi67_23 = HFloat(88.48304651883907), Phi78_17 = HFloat(-1.1056208585550253e-4), Phi78_19 = HFloat(-6.280613803912376), Phi78_21 = HFloat(43.70358459671294), Phi78_23 = HFloat(3.3512347369068123), Phi89_17 = HFloat(-0.28331633046782223), Phi89_19 = HFloat(-6.188115079131574), Phi89_21 = HFloat(-13.925822440131276), Phi89_23 = HFloat(5.203255725370639), Phi910_17 = HFloat(-0.32900739184223554), Phi910_19 = HFloat(-7.185809707858909), Phi910_21 = HFloat(-16.1669897157201), Phi910_23 = HFloat(6.042912241714565)]]

 

Warning, limiting number of major iterations has been reached

 

[0.241435876715391214e-2, [D10 = HFloat(0.001988900799741625), D11 = HFloat(6.700826330455615e-4), D14 = HFloat(0.0022108174875126335), D15 = HFloat(0.001298523707970689), D17 = HFloat(0.0013459282038351902), D19 = HFloat(1.7662427832980496e-4), D2 = HFloat(0.0012416955747892842), D20 = HFloat(5.0643622452156546e-5), D22 = HFloat(0.009580185320006616), D23 = HFloat(0.0022983853492699385), D25 = HFloat(7.03200964138363e-4), D26 = HFloat(8.495269430887406e-4), D28 = HFloat(1.481515200815612e-4), D3 = HFloat(4.443903053331545e-4), D6 = HFloat(2.2200591124570294e-4), D7 = HFloat(0.01005915629472498), Phi1011_17 = HFloat(0.886883773931914), Phi1011_19 = HFloat(-4.609980737722373), Phi1011_21 = HFloat(-104.60164775805397), Phi1011_23 = HFloat(18.5690397972046), Phi1112_17 = HFloat(-6.025408928913167), Phi1112_19 = HFloat(-20.634426393263894), Phi1112_21 = HFloat(-242.95585164839656), Phi1112_23 = HFloat(107.58403867034505), Phi1213_17 = HFloat(-80.57275733334271), Phi1213_19 = HFloat(-22.882545413535468), Phi1213_21 = HFloat(-204.98204126827963), Phi1213_23 = HFloat(-53.249651948908294), Phi12_17 = HFloat(225.74078276804795), Phi12_19 = HFloat(60.952879518061174), Phi12_21 = HFloat(256.2282457438418), Phi12_23 = HFloat(245.86494164793925), Phi1314_17 = HFloat(-4.028668938494458), Phi1314_19 = HFloat(-0.4143116598740555), Phi1314_21 = HFloat(-10.248057903128945), Phi1314_23 = HFloat(-2.6679364992389254), Phi1415_17 = HFloat(-49.255106866735815), Phi1415_19 = HFloat(-6.820689129082128), Phi1415_21 = HFloat(-45.64604945651936), Phi1415_23 = HFloat(-75.03143763880095), Phi1516_17 = HFloat(-135.5797611549774), Phi1516_19 = HFloat(-17.480501715400774), Phi1516_21 = HFloat(-95.26809008906258), Phi1516_23 = HFloat(-134.46896755575636), Phi23_17 = HFloat(135.06987821521912), Phi23_19 = HFloat(20.64805582778886), Phi23_21 = HFloat(129.79718467316934), Phi23_23 = HFloat(73.32903011011902), Phi34_17 = HFloat(214.64764738975026), Phi34_19 = HFloat(65.92781053509862), Phi34_21 = HFloat(158.86989360837936), Phi34_23 = HFloat(62.196984715694455), Phi45_17 = HFloat(32.55413149833871), Phi45_19 = HFloat(28.26026735166782), Phi45_21 = HFloat(45.691055952991924), Phi45_23 = HFloat(13.035682295568991), Phi56_17 = HFloat(8.138977809310253), Phi56_19 = HFloat(7.623185198712813), Phi56_21 = HFloat(11.422405767530362), Phi56_23 = HFloat(3.2585259211177986), Phi67_17 = HFloat(151.28346131545362), Phi67_19 = HFloat(-34.288465707205596), Phi67_21 = HFloat(373.42266757833), Phi67_23 = HFloat(87.41340860226272), Phi78_17 = HFloat(1.0724693768209936), Phi78_19 = HFloat(-0.8681124789323473), Phi78_21 = HFloat(43.074821505878106), Phi78_23 = HFloat(3.262784200274733), Phi89_17 = HFloat(0.4517463493486114), Phi89_19 = HFloat(-0.8499083612028935), Phi89_21 = HFloat(-10.660554777884979), Phi89_23 = HFloat(5.034270143793161), Phi910_17 = HFloat(0.5246496411398992), Phi910_19 = HFloat(-1.1110703087821467), Phi910_21 = HFloat(-12.476983437436498), Phi910_23 = HFloat(5.846600456480841)]]

(1)

rep3 := Optimization:-LSSolve(list_2);

Warning, limiting number of iterations reached

 

[0.173158517170151e-3, [D18 = HFloat(4.940399971912688e-4), D27 = HFloat(3.123022055695685e-4), D29 = HFloat(0.012053079546584859), D36 = HFloat(-0.02877258006663398), Phi1011_18 = HFloat(-10.047440876654184), Phi1011_24 = HFloat(-4.89440481171543), Phi1112_18 = HFloat(-58.20988683846057), Phi1112_24 = HFloat(-29.338754572925264), Phi1213_18 = HFloat(-83.08949138713488), Phi1213_24 = HFloat(-56.12110719411204), Phi12_18 = HFloat(260.26860758672166), Phi12_24 = HFloat(263.20617377185454), Phi1314_18 = HFloat(-3.4457440384526827), Phi1314_24 = HFloat(-2.5087129560312764), Phi1415_18 = HFloat(-32.549935720740386), Phi1415_24 = HFloat(-25.48888809044497), Phi1516_18 = HFloat(-72.10104065466892), Phi1516_24 = HFloat(-62.675427552811946), Phi23_18 = HFloat(89.12923249096826), Phi23_24 = HFloat(61.483196901345444), Phi34_18 = HFloat(125.54530604503496), Phi34_24 = HFloat(68.67880838477211), Phi45_18 = HFloat(13.68255794893495), Phi45_24 = HFloat(8.04368413990329), Phi56_18 = HFloat(2.919597731660753), Phi56_24 = HFloat(1.8129257306833229), Phi67_18 = HFloat(24.210838638941365), Phi67_24 = HFloat(20.01451125467167), Phi78_18 = HFloat(-1.4701283478035998), Phi78_24 = HFloat(-0.6483549804256254), Phi89_18 = HFloat(-2.559740721298145), Phi89_24 = HFloat(-1.2378238951622822), Phi910_18 = HFloat(-2.938954535723053), Phi910_24 = HFloat(-1.4230779501001942)]]

 

Warning, limiting number of major iterations has been reached

 

[.131366159495864665, [D18 = HFloat(5.8286890720965e-4), D27 = HFloat(3.9234945275683984e-4), D29 = HFloat(-2.9726900078315097), D36 = HFloat(2.9754382616693835), Phi1011_18 = HFloat(1.0203797883189318), Phi1011_24 = HFloat(1.0079659723850967), Phi1112_18 = HFloat(1.1092824529457799), Phi1112_24 = HFloat(1.0312364864731682), Phi1213_18 = HFloat(1.1885498329489976), Phi1213_24 = HFloat(1.0066208797233582), Phi12_18 = HFloat(0.37042074900702604), Phi12_24 = HFloat(0.5964090885427208), Phi1314_18 = HFloat(1.0117497255740284), Phi1314_24 = HFloat(1.0030476511529876), Phi1415_18 = HFloat(1.0837133285911669), Phi1415_24 = HFloat(0.9954012128357053), Phi1516_18 = HFloat(1.2352965407494056), Phi1516_24 = HFloat(1.1798017873822157), Phi23_18 = HFloat(0.7097521011632247), Phi23_24 = HFloat(0.8837922498478785), Phi34_18 = HFloat(0.8337073271457043), Phi34_24 = HFloat(0.9624593994431938), Phi45_18 = HFloat(0.98503438366238), Phi45_24 = HFloat(0.999914375477929), Phi56_18 = HFloat(0.9980930649027827), Phi56_24 = HFloat(1.0021135395285548), Phi67_18 = HFloat(0.9845976365161451), Phi67_24 = HFloat(1.014443548053698), Phi78_18 = HFloat(1.0051330190684922), Phi78_24 = HFloat(1.0037813447548194), Phi89_18 = HFloat(1.00719778341775), Phi89_24 = HFloat(1.0042421054593569), Phi910_18 = HFloat(1.007965454569841), Phi910_24 = HFloat(1.00446826792656)]]

(2)

rep4 := Optimization:-LSSolve(list_2, {D29 >= 0, D36 >= 0});

Warning, limiting number of major iterations has been reached

 

[0.100051838928020805e-2, [D18 = HFloat(5.151431335854967e-4), D27 = HFloat(3.1711720791143086e-4), D29 = HFloat(0.005605722295400071), D36 = HFloat(0.0020169300763647198), Phi1011_18 = HFloat(-10.518110711035249), Phi1011_24 = HFloat(23.69087038554566), Phi1112_18 = HFloat(-27.770981377847882), Phi1112_24 = HFloat(31.78453691201447), Phi1213_18 = HFloat(-53.11317389788098), Phi1213_24 = HFloat(-81.22637858019229), Phi12_18 = HFloat(275.1761250113166), Phi12_24 = HFloat(204.95439428055727), Phi1314_18 = HFloat(-10.558267219470897), Phi1314_24 = HFloat(17.424914531663973), Phi1415_18 = HFloat(-33.08225878870024), Phi1415_24 = HFloat(-38.506622367871806), Phi1516_18 = HFloat(-88.70176908241132), Phi1516_24 = HFloat(-110.52221215837501), Phi23_18 = HFloat(116.90860096214605), Phi23_24 = HFloat(82.85584178444226), Phi34_18 = HFloat(87.86056251202156), Phi34_24 = HFloat(31.709695488160147), Phi45_18 = HFloat(23.344387968260033), Phi45_24 = HFloat(-23.939720201787527), Phi56_18 = HFloat(-0.4573382415630158), Phi56_24 = HFloat(11.85096445100273), Phi67_18 = HFloat(50.00293380451404), Phi67_24 = HFloat(43.15119566212074), Phi78_18 = HFloat(-1.746203841852997), Phi78_24 = HFloat(5.792224317605069), Phi89_18 = HFloat(-8.233370777799111), Phi89_24 = HFloat(21.321096735379633), Phi910_18 = HFloat(-8.345823587878009), Phi910_24 = HFloat(21.4610282086855)]]

(3)

``


 

Download worksheet_help_-_TS.mw

For most custom-made packages, the relevant files need to be saved in/extracted to

C:\Program Files\Maple 2016\lib

on Windows or

//Library/Frameworks/Maple.framework/Versions/2016/lib

on Macintosh.  If you are not familiar with using these directories (which can interfere with your installation of Maple if modified incorrectly), you can also install an add-on package to a specified user library folder:

If you have previously added a custom user library path to 'libname', please copy the library file “%PackageName%.mla” to that folder. If you have not, you may add that user library folder by doing the following:

  1. Copy “%PackageName%.mla”  to a user library folder. If you do not have a user library folder, simply create a folder for your Maple library files. This folder may be located anywhere on your computer, though a suggested path is to create a new folder called MapleLib in your "home directory" and place the files there. To find your home directory in Maple, use kernelopts( homedir ) ; For more information, see the kernopts help page.
     
  2. Open your Maple initialization file. See the initialization file help page for details on how to create and configure these initialization files.
     
  3. Add the line: libname := libname, "PathToUserLibraryFolder” : For example, to add the path to your "MapleLib" folder under your home directory use libname := libname, cat( kernelopts( homedir ), kernelopts( dirsep ), "MapleLib" ) :
     
  4. Save your Maple initialization file and restart Maple.

 

From there, you should be able to load the package using the with() command.  For more information, see the with command help page.

 

It looks like those particular values of m, N, and h are outside of the range that Maple can currently handle for this problem; the closest set of values I can get an answer for is m=2, N=17, h=0.41600. 

In most cases, even using N=16 results in an echoed call to fsolve(), so it might be difficult to find a way to get these values any closer to what you're hoping for, but hopefully this set gives you a good starting point. You may also want to consider trying to rewrite the equations in an different form.

You're right, this is an issue that our developers hope to have resolved for 2017.2.  In the meantime, please try launching Matlab via the 'matlab' script rather than the app bundle.

Thanks for the post.  Can you send us a copy of this worksheet?  Please also try some simple computations (such as 1+1), and test cmaple.  -support@maplesoft.com

1 2 Page 2 of 2