Hello,

I want to simplify a fraction by cancelind the powers of ten. When the factors in front of the powers are short, it works fine with 'simplify', but with such long decimal numbers, Maple does not want to cancel down. How can I force it?

Example:

Regards, Yader

I do not understand why I am getting permission denied trying to create a file.

restart;
fd := FileTools[Text][Open]("C:\\foo.txt", create=true,overwrite=true);

     Error, (in FileTools:-Text:-Open) permission denied

I am running the above on my PC, from Maple, using the same user which is me. I can ofcourse create a file on my C drive outside of Maple. So I do not understand why it says permission denied. Also tried "C:/foo.txt". Also tried

         fd := FileTools[Text][Open]("C:\\foo.txt", create,overwrite);

Any idea what can be the problem? and how to determine why it is giving this error? I am on windows 7, using Maple 2015.

http://www.maplesoft.com/support/help/maple/view.aspx?path=FileTools%2fText%2fOpen

http://www.maplesoft.com/support/help/maple/view.aspx?path=IO_errors

"permission denied" - An attempt was made to read from a file on which the user does not have
read permissions, or write to a file on which the user does not have write permissions."

Well, I do have write permission to make a file on my C drive !

I tried using fopen, and got same error

restart;
  try     
     fd := fopen("C:\\foo.txt",APPEND,TEXT);
  catch:
     print(`Unable to open file, error is`);
     print(StringTools:-FormatMessage(lastexception[2]));
  end try:

                 Unable to open file, error is
                      "permission denied"

So I have no idea why I keep getting permission denied.

plz help me, I m trying to solve homotopy analysis method for lane emden, what additional steps I have to taken in above programming?

Download LE_EQ.mw

Hi,

  Excuse me, I have the following code

###########################

with(combinat);
with(LinearAlgebra);
f:=([1,2,3]);
permute(f);

n:=3;
g:=array(1..n);

for i from 1 to n do
  g[i]:=i;
end do;

permute(g);

###########################

 The output ends up with

Error, (in combinat:-permute) 1st argument must be a list, set or a non-negative integer

 It seems that maple will not generate permutation of g... How to generate the permutation of array g? i.e. by assigning all variables automatically.

Thank you very much!

Hi, 

     I was just curious about the difference between defining tensors as arrays/matricies/TensorArrays vs defining them as algebraic symbols. I found that defining them as an expression lead to the wrong answer, and I was forced to define a tensor (LKh) as an array. I've attached a worksheet demonstrating my problem.

I apologize for the amount of tensors needed to find this problem, but it is the only one I have reproduced the issue. I basically define the metric
Metric = g_
auxillary tensor = h
Killing vector = K
LieDerivative of h, wrt K = LKh (not a tensor array)
LieDerivative of h, wrt K = LKh2 (tensor array)
Then I compare two expressions, rho1 and rho2 computed from LKh and LKh2 and they disagree. 

MapleQuestion.mw

Thanks!

I ran into a problem mentioned here before: With lots of code in one execution group, the standard GUI editor becomes very slow at least in older versions of Maple (I have reason to believe that this is mostly fixed in Maple 17 and, presumably, newer). With one of my packages this has become enough of a problem to force me to deal with it (and Modules and Packages have to be all in one execution group). I have some machines that cannot go beyond Maple 15 so I keep my stuff compatible with it (so far not a big sacrifice).

The fix is rather easy: In the GUI, export the code as "Maple Input" which creates a file with extension .mpl. Then, in the GUI have a small file that uses "read" to read-in and execute the code (of the module). In my case, I have the necessary LibTools statement to save as a .mla right after the package code ends so, bingo, I have just updated the package. I edit the .mpl file in Emacs.

So I installed Joe Riel's Emacs mode maplev-mode.el. This works rather well. However, it took not long for it to show its age: it does not recognize the multi-line comments (*...*). This is kind-of bad as this forces me to not use quotes in the comments since they mess-up the keyword highlighting esp. when using single quotes as in "I'll" and so on. Likewise, maplev-indent-buffer messes up the whole buffer as it interprets the multiline comments as code.

So my question is: is someone still working on this and using it?? I guess I can work around some of these things, but it would be nice if this worked.

Incidentally, my version of maplev-mode.el seems to be dated June 2005.

I know Joe is here on occasion, maybe he can chime in. It is not absolutely fatal; I could forgo using maplev, but if I am forced to go this route of readin in the package; using maplev-mode at least gives me something in return. Unfortunately my elisp skills are not quite up to fixing this myself; at least not in a reasonable time.

If it is relevant: I use Emacs 22.1 on Mac OS X. Yes, I could upgrade to at least Emace 23.4, if that were the issue.

Thanks in advance,

Mac Dude

I have the following paper:

http://projecteuclid.org/download/pdf_1/euclid.dmj/1077490637

Now I wanted to check Fritz John's claim in the proof of Theorem 1.1, he says that equation (7) can be easily verified for case i=1,k=2.

Now at first I tried to calculate by hand, but it's just a lengthy calculation, so now I turned to maple to check its validity, I get that this claim is false, am I wrong in my code? if yes, then how to change it?

P.S

I changed between xi and x and eta and y.

In the following is the code:

Download ultrhyperbolic_pde.mw

Download ultrhyperbolic_pde.mw

As part of a project, I am numerically estimating the roots of many large polynomials. Occasionally, "fsolve" fails with strange errors related to "fsolve/refine2". Searches for these error messages have turned up nothing.

I've inluded the code below that causes the error on Maple 18.02. I apologize for the polynomial in question being so long, it's the shortest example I have. The error it generates is:

    Error, (in fsolve/refine2) invalid input: evalf expects its 2nd argument, n, to be of type posint, but received undefined

Is this a bug? Or am I missing some fsolve option to prevent this? Note that it only happens when the "complex" flag is used.

=========

P := x^14-22702264347017701018473605850972699930097274504938699916055555261201515180511538865331807292689345943133521696082918467714371257277276696385067641909170155322906230250853577229812913946663078548646992393337618113886746876557117483839533553328895358682670189394678910311793504505447628428181885141769168591937690303328913335175451328463754619536253583902806843310134957600949886784187209785783810122275010505534415815566439121541947044486358488039865870455952098827525405324562601732796858645293515431747164008309785658410612354201118685855495413079021176507985235094746401708925593687656572387531020719291601076812080687859808747213536777976702071405128537760507468013438105233313663196919816564525291458692028177366393652501832447863872200682143768513389322886600569382594287138458765510827267842205096062437750804878586024353928794905249283675708441066101095406513448522689302522442783437142289641259057413952301148939774149714785/3195755849586795631956816504521213454239300164039404772924331154185577854140658969534719471406093912112781063157828311505891258148739680804289213024862131311540960306206602785748866445362483281617891374949555209869677857419473553982132073059025609434698683760348542259396937054082293168625919023158753310878489047944378154369352523436731294817697449949932655665007647918855300664365159027040571937825740235967492228453331261542499260943085539271304638576578246276634403350307801994081681247214869084246168101721298760198550961832560608341435638093413744839736250679074198753022491225840288065341597851066663786665723409977381822591654466626645542917017628998630902708076612502066607817250779545511895971357711983287763127653752300554550391349040027472903180009282594974618980021621163037989247901106508257414514187962209356325857887950302223210328647697948055097831009797738621154319922212951316644741457327450027692469090867369598*x^13+19088859498864751331345860430721481446264521641744903691362655800372349704990331481604867685549645823662978708926030541236042951546550966879612333115396628902751820387999904934599090760358886795430484312266737008386396041896213971568537362589851823779617430207078749426022658232278071527361264481524611089324107754031784837527081637219350016169914382322455035364613935875393571579561406195287363628553419822536428710010055920488818415526206620047517917895155637033562338042275152771173240104076821411360366799172066699958543868065037999702280159896040588223787643434915579465270491451613199185385049196526456210057933748521047167538262357063585093474544299142560492581751607753970282443057122762426600024763892341448332834018680513343674283251162037067303651651086278409136799357849452879897251530675098741236156640469815784447341282424004221641529187217962536022784563163918511210513153785881467158114512281634789894107114727680109/1065251949862265210652272168173737818079766721346468257641443718061859284713552989844906490468697970704260354385942770501963752716246560268096404341620710437180320102068867595249622148454161093872630458316518403289892619139824517994044024353008536478232894586782847419798979018027431056208639674386251103626163015981459384789784174478910431605899149983310885221669215972951766888121719675680190645941913411989164076151110420514166420314361846423768212858859415425544801116769267331360560415738289694748722700573766253399516987277520202780478546031137914946578750226358066251007497075280096021780532617022221262221907803325793940863884822208881847639005876332876967569358870834022202605750259848503965323785903994429254375884584100184850130449680009157634393336427531658206326673873721012663082633702169419138171395987403118775285962650100741070109549232649351699277003265912873718106640737650438881580485775816675897489696955789866*x^12-127688837609696458957114129756229560761957972259253280819996356067917173759565012901801561924391178368568146719627801670086606489531437386224078360185442651606983719684283163392876990522586784115059551865746707609765679864632874671595399416688286257053075135779925094175440416074968471245768830366824397599424731191899057489251725430472639828977416853808059394673266682604308077331301860791811476274942568803494246399367164616630866928631772760003749091917886558963952047434319195736393271420111064778587861639539510320744497931007588784407172972776901653630399291814617861650330433072614870207218474263898528043868017109168847074788133295715653324601280999334137328493510780499508083274179117783232296907665583279993325725716354393277745170409349317876378784871325009748734263290375761397883657890413900529632709410443413043575189427898559331856967020187201932742096158736566419271039506140015010172468151681141071869870925420155369/6391511699173591263913633009042426908478600328078809545848662308371155708281317939069438942812187824225562126315656623011782516297479361608578426049724262623081920612413205571497732890724966563235782749899110419739355714838947107964264146118051218869397367520697084518793874108164586337251838046317506621756978095888756308738705046873462589635394899899865311330015295837710601328730318054081143875651480471934984456906662523084998521886171078542609277153156492553268806700615603988163362494429738168492336203442597520397101923665121216682871276186827489679472501358148397506044982451680576130683195702133327573331446819954763645183308933253291085834035257997261805416153225004133215634501559091023791942715423966575526255307504601109100782698080054945806360018565189949237960043242326075978495802213016514829028375924418712651715775900604446420657295395896110195662019595477242308639844425902633289482914654900055384938181734739196*x^11+1941963889284143967630503461641384772639246155223080045213834947146248802376220874585881737054605033501866020180580996666713652795880193379326169002867732637532200447194534846339338413543240801939477478241099683412186213038204202290858666453453417846899586996164402928265510496311234255565224399736489137714957014062613618467711330149177700905620590617841256796029108309659216987574764436490494934048714919670190684046029243176314833939438755957046700890497104507387346236050782952954331963365980088907386124713398155694605596359102946980291494749083679360161061519037964676115079433007443504537411864172301459893087256329861985994656612965817883540871319790509064913361633903111901088284675188114992527367188875256164648035419067179258498467050438971237757123876227876902374176109894916835963212462977715488403210262610643862278435267351697867431486692646214503482828458653994117921039913207766285237066798400775441411774079610837/177541991643710868442045361362289636346627786891078042940240619676976547452258831640817748411449661784043392397657128416993958786041093378016067390270118406196720017011477932541603691409026848978771743052753067214982103189970752999007337392168089413038815764463807903299829836337905176034773279064375183937693835996909897464964029079818405267649858330551814203611535995491961148020286612613365107656985568664860679358518403419027736719060307737294702143143235904257466852794877888560093402623048282458120450095627708899919497879586700463413091005189652491096458371059677708501249512546682670296755436170370210370317967220965656810647470368146974606500979388812827928226478472337033767625043308083994220630983999071542395980764016697475021741613334859605732222737921943034387778978953502110513772283694903189695232664567186462547660441683456845018258205441558616546167210985478953017773456275073146930080962636112649581616159298311*x^10-8863297827898165839415750496524113595646716762121322844115735229732707054220168863570410233048901583522983420448394277638284549018035513758104914498710641607628947697242244219841860993879788358381735454419316105861594264938271360192839234405577377200072500528683390799739067094807744149646139901716484318789624752809347475833611242793680607929117314808486761679788190886930919706651326101755691947512243418216460237241769969836429891299366314558409971924025105037908119564861578530568458046073973392717591352587749276788404403969612641258746844907150027624513801294297807549646148499291417732004938891897241268904655286649579390754923163271693343760744341965406109121853934340673760793157027856466007930099451628400810185712964495614693463085689905135333950374874743558356849483930901197948196695521171654728867080261994068527012182974920030328883744764966738189985551194012023576203112120540148697524837372858043563451876509542849/3195755849586795631956816504521213454239300164039404772924331154185577854140658969534719471406093912112781063157828311505891258148739680804289213024862131311540960306206602785748866445362483281617891374949555209869677857419473553982132073059025609434698683760348542259396937054082293168625919023158753310878489047944378154369352523436731294817697449949932655665007647918855300664365159027040571937825740235967492228453331261542499260943085539271304638576578246276634403350307801994081681247214869084246168101721298760198550961832560608341435638093413744839736250679074198753022491225840288065341597851066663786665723409977381822591654466626645542917017628998630902708076612502066607817250779545511895971357711983287763127653752300554550391349040027472903180009282594974618980021621163037989247901106508257414514187962209356325857887950302223210328647697948055097831009797738621154319922212951316644741457327450027692469090867369598*x^9+2017798632508126214178032687554207106376957218619840832956662086046164665731840663616010195125899727795939076982732552766119787016500832211715921717230052496681979577432852310401087868502196643627732719648105718026988556567683519630485815400374479891004302010335606831102703848328584097945978508243431022893640963113656761657600197417501485533198377194139684052176343994472985299657953300469363286556184526802586584214755307697158247107887873091236237632599123600602905959806309348139033950376959229535155791803352079139955123963320105932227644780215619966766055108655288741020886209290486474711298006297607280189057416424593467478571219372613281308689321078436254551057742499063473864232694061267319983756694814398288262834174390168011042951364931686132193005284966871663010432330119034962795839890899090501481880553750048136836870884603174737258837536327019963547476613035420193242235958864072191310211266838717111661997427211841/12783023398347182527827266018084853816957200656157619091697324616742311416562635878138877885624375648451124252631313246023565032594958723217156852099448525246163841224826411142995465781449933126471565499798220839478711429677894215928528292236102437738794735041394169037587748216329172674503676092635013243513956191777512617477410093746925179270789799799730622660030591675421202657460636108162287751302960943869968913813325046169997043772342157085218554306312985106537613401231207976326724988859476336984672406885195040794203847330242433365742552373654979358945002716296795012089964903361152261366391404266655146662893639909527290366617866506582171668070515994523610832306450008266431269003118182047583885430847933151052510615009202218201565396160109891612720037130379898475920086484652151956991604426033029658056751848837425303431551801208892841314590791792220391324039190954484617279688851805266578965829309800110769876363469478392*x^8+595653463870329564049793978610653077784924147310081049725886229514496074424431896371139535290512513543370421669991943879022295484566837982193741714293063648577558130603631634766521328149129569969823197271545558633437311107596551672056665506079114364599202039865730405948002218769139583348447586421596752286092282439204649127363793862275710130284296695116082940485912870421340555598970364233037216088883415900477944075155274644931796851988347845481689494432787436457416213309975377477034435315885301330073984386960225656990818706334595431458471582550550165189482320607945492496296417100187628303765229872504055798768101834332864953805589825320224343693907336414200688255911370715097612333973794836212429385145521765825313810064172829130430532871221590834385381629229656091998869097712301634818490602700217107948854487937370243140870821241350717775059492175130737668904369695266161667460527553083982270988490551570900042454610585889/12783023398347182527827266018084853816957200656157619091697324616742311416562635878138877885624375648451124252631313246023565032594958723217156852099448525246163841224826411142995465781449933126471565499798220839478711429677894215928528292236102437738794735041394169037587748216329172674503676092635013243513956191777512617477410093746925179270789799799730622660030591675421202657460636108162287751302960943869968913813325046169997043772342157085218554306312985106537613401231207976326724988859476336984672406885195040794203847330242433365742552373654979358945002716296795012089964903361152261366391404266655146662893639909527290366617866506582171668070515994523610832306450008266431269003118182047583885430847933151052510615009202218201565396160109891612720037130379898475920086484652151956991604426033029658056751848837425303431551801208892841314590791792220391324039190954484617279688851805266578965829309800110769876363469478392*x^7-15113570599033390421079877152704282575760402634458053122854896745489277827307511104729555317925906225568214914439058228407909604233709303373800635755523163622556864164238666772567187619354483299953911709839860063588020162857990784248028866643232820765874877631054265552390768506568101602491000432124326227024341745780450739624272135251182271133507810838226637332261885858935760923125210938994349370924860536080811884391229918127119531006177386365479201266087848843332916814370814507988299037716701840225002363896766748196913932600780960329893902502167652386550281812383391521738502785942278676540705757260537410521405647255640494560579216704794939121721257870245415144651904709474211768371408204914862883986252875760943819790296200235921138992713190368600609580449924683146430563973244161071950255162230256773426228399168919888429747422610165410976294531879556935995391906949049683251359274465900646212494068580346654890616742327/25566046796694365055654532036169707633914401312315238183394649233484622833125271756277755771248751296902248505262626492047130065189917446434313704198897050492327682449652822285990931562899866252943130999596441678957422859355788431857056584472204875477589470082788338075175496432658345349007352185270026487027912383555025234954820187493850358541579599599461245320061183350842405314921272216324575502605921887739937827626650092339994087544684314170437108612625970213075226802462415952653449977718952673969344813770390081588407694660484866731485104747309958717890005432593590024179929806722304522732782808533310293325787279819054580733235733013164343336141031989047221664612900016532862538006236364095167770861695866302105021230018404436403130792320219783225440074260759796951840172969304303913983208852066059316113503697674850606863103602417785682629181583584440782648078381908969234559377703610533157931658619600221539752726938956784*x^6-56606911212746288531797425095670805062399424421101677717862600827392887966086570498590673119613746469165156337523381795636321677036991477302420753674709713686811882568078016975400881624969306258459062007649906773777809876859663876646222545569178638387083155806346202130570988684784830435390865363436508319229238598519437586054312503598853928438972636264615238650707643552964679259318243634317570183872115416681752259271751432389553718574295517207917784228757458885377537383478835886996256222286223854172678118278099547239528535908757605400918412542872093301940268914810937904361836991602742563582312675452630601354931494058266771375234065791092815667038632389781434772470040614489288147052306960907419599693648698532033457556243783702691074860973568369668485091260383950009098762596509973546964151952489551510210947116627078116642234793376887344063164400839769932275480710548974537201716359604788112690110055153211687060319915667/25566046796694365055654532036169707633914401312315238183394649233484622833125271756277755771248751296902248505262626492047130065189917446434313704198897050492327682449652822285990931562899866252943130999596441678957422859355788431857056584472204875477589470082788338075175496432658345349007352185270026487027912383555025234954820187493850358541579599599461245320061183350842405314921272216324575502605921887739937827626650092339994087544684314170437108612625970213075226802462415952653449977718952673969344813770390081588407694660484866731485104747309958717890005432593590024179929806722304522732782808533310293325787279819054580733235733013164343336141031989047221664612900016532862538006236364095167770861695866302105021230018404436403130792320219783225440074260759796951840172969304303913983208852066059316113503697674850606863103602417785682629181583584440782648078381908969234559377703610533157931658619600221539752726938956784*x^5+8803574970613871156806085396512138961742521948396413545070976399600915984183390805096731313217646981222548657773275592950616488668094770447693809833164459043946592009996136155728311094386907899766270966922695038949490311577603551489122618956774441280490375954836897874880872800831317477259323488993074897412948399841956182103363590355763668471507912443956969097316815558683861211173548132219642059999985419005273573000030449164050583320511714605312903276113074336198985650814714729654266427222607355408552351971127995723711943667875256201812876969068878200785868067258284670407922968682007482594553643052049155554750878922972932169341172726102387357444748432711229584365031020816215817387958789462262033872746673386146274342628377930132499351532652360977700457249013199924702013517630294342824985973844030539327838892887509452469632329879859377436506301311966132764081878323235415242489663285064434254301628330944024796246782645/25566046796694365055654532036169707633914401312315238183394649233484622833125271756277755771248751296902248505262626492047130065189917446434313704198897050492327682449652822285990931562899866252943130999596441678957422859355788431857056584472204875477589470082788338075175496432658345349007352185270026487027912383555025234954820187493850358541579599599461245320061183350842405314921272216324575502605921887739937827626650092339994087544684314170437108612625970213075226802462415952653449977718952673969344813770390081588407694660484866731485104747309958717890005432593590024179929806722304522732782808533310293325787279819054580733235733013164343336141031989047221664612900016532862538006236364095167770861695866302105021230018404436403130792320219783225440074260759796951840172969304303913983208852066059316113503697674850606863103602417785682629181583584440782648078381908969234559377703610533157931658619600221539752726938956784*x^4-523859084082833471446711755599121052886811016364975130659544060359420880174405546029952820794428450698789941159230333192978458494119123793478183362395347411263910209528686607499619732076210935923583757530528475415340759837176163831540648474842441088111430598419604346899579186230546208967364447721656454914767406049990213818512260948866882128853891871133892169918212832229057031510780243573302325692630372815883421975596513331009467248542883004381175675544983224650496679943101389848155364127861671314139401868335830005933182296640353260635533021407833519995789902759067978323853832610509380765976726389069758315717268633025536104645794437969501284514833324641003339055262908646594490675898481723072173169073269959496697492841280565500876977601579793290422081302241207973528760098915143648721380740535972297119460582597942660098036955825437773842971660604713909242648770501670101799996143752912217696962275386368069113806310933/25566046796694365055654532036169707633914401312315238183394649233484622833125271756277755771248751296902248505262626492047130065189917446434313704198897050492327682449652822285990931562899866252943130999596441678957422859355788431857056584472204875477589470082788338075175496432658345349007352185270026487027912383555025234954820187493850358541579599599461245320061183350842405314921272216324575502605921887739937827626650092339994087544684314170437108612625970213075226802462415952653449977718952673969344813770390081588407694660484866731485104747309958717890005432593590024179929806722304522732782808533310293325787279819054580733235733013164343336141031989047221664612900016532862538006236364095167770861695866302105021230018404436403130792320219783225440074260759796951840172969304303913983208852066059316113503697674850606863103602417785682629181583584440782648078381908969234559377703610533157931658619600221539752726938956784*x^3+2507277613118701025793992364166478343634113591260892474295326766388926595401380659553418488612372800473078429284383414164100295497618696922175584194648519640315288543669474665647481253042089841838880845698663503017597568567494036511285649062160850714939388491186949845790680376696816908809272981412907183537044694180946810281215322367828144494392069154586949731736984686662811664798321227156168742191864566896432033615316362529844704143667580907339601101101040983407725340793380960622046484034595588591168465069176540166697625520589666927661660285415510869031443183270437121086747446353056448389264226941586709914236711688346438937447603056458249768319677721039404318822571940087776493501534291441156924207226385372828419997653921508823109486892013281164698869726954506626852453903751862417478762888542286370388085783793727357089951716036021380037300853575906423073257325966464463919011630179868083502189793561658823882854823/5681343732598747790145451563593268363092089180514497374087699829663249518472282612506167949166389177089388556725028109343806681153314988096514156488643788998295040544367293841331318125088859167320695777688098150879427302079064095968234796549378861217242104462841852905594554762812965633112744930060005886006202751901116718878848930554188968564795466577658054515569151855742756736649171603627683445023538197275541739472588909408887575009929847593430468580583548936238939289436092433922988883937545038659854403060086684797423932146774414829218912166068879715086667873909686672039984401493845449496173957451846731850174951070901017940719051780703187408031340442010493703247311114785080564001385858687815060191487970289356671384448534319200695731626715507383431127613502177100408927326512067536440713078236902070247445266149966801525134133870619040584262574129875729477350751535326496568750600802340701762590804355604786611717097545952*x^2+409625955662068882951029126793688739518090138239218488984796795645801325262353406586076439055165056320023643100727785301142632216544301080743235707086824485890205950529700657617438269917459707919268974737411575302005751244657214481075183329087253965790410699412397427442839741206848483871370551438808586427029387351271518650463772184256590372936207478553171485915519400489218737299828282527927029698794411869349982447984333764576483854239683590209335608900236726346302994114690056427772581682505305971250477006390651418491364024408115184552516886617401858022820021648328348447629622760279509245058592997481694825738991388860444669425663371096552711454623191875274037228193231324209339530305072938805229134963657645497492583194974100090696050921013395443204011245764738501052755929663789410609492539391491688417079601608018345617616844409309669918907009896517735188492371103865855756475811783862809885417127411425999146175/25566046796694365055654532036169707633914401312315238183394649233484622833125271756277755771248751296902248505262626492047130065189917446434313704198897050492327682449652822285990931562899866252943130999596441678957422859355788431857056584472204875477589470082788338075175496432658345349007352185270026487027912383555025234954820187493850358541579599599461245320061183350842405314921272216324575502605921887739937827626650092339994087544684314170437108612625970213075226802462415952653449977718952673969344813770390081588407694660484866731485104747309958717890005432593590024179929806722304522732782808533310293325787279819054580733235733013164343336141031989047221664612900016532862538006236364095167770861695866302105021230018404436403130792320219783225440074260759796951840172969304303913983208852066059316113503697674850606863103602417785682629181583584440782648078381908969234559377703610533157931658619600221539752726938956784*x+19274602981612147290159524985559484093130831751480756525782715054667194973369495232408312165313759674110576759872392612089012272923859625516884268120491350321412487913144106351173483928355210063578761326082759699790669852993610922334255512904908344034497486576749204467479440080193445816812437124815705534871107288254095868559730608571611775333291508638854525168655723328657110979842065373668113720199606612577706540055748716804695554224139209569128003744375141059759616450648545383218233710479963101078863006113802889357950673977857921892481030024186104599766895886493039433820102404772715112824244120718747330555481069366856972350416806635004205386220265273799023113455715583543586661235789661857624566216488236144677299497354093058902257519792639035715396313400079080186425228931237501312897851917238249279942091650204339882709384865082513445071486731865811778705191906577495950612446972346017024882968602859856464/1597877924793397815978408252260606727119650082019702386462165577092788927070329484767359735703046956056390531578914155752945629074369840402144606512431065655770480153103301392874433222681241640808945687474777604934838928709736776991066036529512804717349341880174271129698468527041146584312959511579376655439244523972189077184676261718365647408848724974966327832503823959427650332182579513520285968912870117983746114226665630771249630471542769635652319288289123138317201675153900997040840623607434542123084050860649380099275480916280304170717819046706872419868125339537099376511245612920144032670798925533331893332861704988690911295827233313322771458508814499315451354038306251033303908625389772755947985678855991643881563826876150277275195674520013736451590004641297487309490010810581518994623950553254128707257093981104678162928943975151111605164323848974027548915504898869310577159961106475658322370728663725013846234545433684799:

fsolve(P,complex);

I would like to get some sort of table (maybe also a plot) that shows me the effect on the function if I change a variable.

For example, I have f(x,y,z)=x+yz

Now I would like to get a list with the results for f if I run z from, say, -10 to +10.

Is it also possible to do this with all variables at the same time?

here a is constant.

Hi,

I am trying to solve a set of ode which depends on some parameters like A_0,0,  A_1,0, A_1,1, B_0,0,  B_1,0, B_1,1, C_0,0 and so on. Here is some part of my code:

restart;

sigma := 1; X := proc (m) if m <= 1 then 0 elif 2 <= m then 1 end if end proc;

lambda := proc (m, k) if k <= m and 0 <= k then 1 else 0 end if end proc;

Typesetting:-Settings(functionassign = false);

for m from 0  to 4 do    

for k  from 4 by -1 to 0 do  

A[m,k](r) :=diff(f[m,k](r),r$2)+((k+1)*(k+2))/(r^(2))*(lambda(m,k+2)*f[m,k+2](r)-f[m,k](r)):  

T[m,k]:=(k+1)*(lambda(m,k+1)*A[m,k+1](r) -lambda(m,k-1)*A[m,k-1](r)):      

S[m,k]:=(k+1)*(lambda(m,k+1)*f[m,k+1]( r)-lambda(m,k-1)*f[m,k-1](r)):               # There are also C[m,k], B[m,k] and E[m,k] definitions similar to A[m,k],T[m,k] and S[m,k]

Q[m, k] := r^(4-sigma)*E[m, k]-2*lambda(m, k+2)*(k+1)*(k+2)*(r^2*(diff(f[m, k+2](r), `$`(r, 2)))-2*r*(diff(f[m, k+2](r), r)))-(k+1)*(k+2)*(k+4)*((k+3)*lambda(m, k+4)*f[m, k+4](r)-(2*(k+1))*lambda(m, k+2)*f[m, k+2](r));

ode[m, k] := r^4*(diff(f[m, k](r), `$`(r, 4)))-(k+1)*(k+2)*(2*r^2*(diff(f[m, k](r), `$`(r, 2)))-4*r*(diff(f[m, k](r), r))-(k-1)*(k+4)*f[m, k](r)) = Q[m, k];

soln[m, k] := rhs(dsolve(ode[m, k], f[m, k](r)))

After obtaining the solution, there are some additional parts in my code to find the coefficients of odes. The code is just working fine until m=3. When m=3, I am getting this error:

Error, (in solve) cannot solve expressions with Int(-(1/282240)*r*(3*(h . R)^2*(Int((6720*ln(r)*h*r^6-22176*ln(r)*h*r^5-50400*ln(r)*r^6-2814*h*r^6+33600*r^7+18144*ln(r)*h*r^4+90720*ln(r)*r^5+8238*h*r^5-95520*r^6+2352*h*ln(r)*r^3-7974*h*r^4+57780*r^5-2880*ln(r)*h*r^2-14400*ln(r)*r^3-2443*h*r^3+2100*r^4+6108*h*r^2+5340*r^3-1115*h-3300*r)/r^6, r))-4480*_C1), r) for _C1

I think the problem is due to declaration of f[m,k](r) values for the set of odes. For example, before solving ode[3,3], the code declares a value for f[3,3](r) which includes some integral definition in ode[3,3] although i want to find the solution for f[3,3](r). 

To illustrate,

ode[1,1]=r^4*(diff(f[1, 1](r), r, r, r, r))-12*r^2*(diff(f[1, 1](r), r, r))+24*r*(diff(f[1, 1](r), r)) = -Typesetting[delayDotProduct](r, h . R, true)*((2*(diff(f[0, 0](r), r))-4*f[0, 0](r)/r)*(diff(f[0, 0](r), r, r)-2*f[0, 0](r)/r^2)-(2*(-3/4-1/(4*r^2)+r))*A[0, 0]+(2*(diff(f[0, 0](r), r, r, r)+4*f[0, 0](r)/r^3-2*(diff(f[0, 0](r), r))/r^2))*f[0, 0](r))

soln[m, k] := rhs(dsolve(ode[m, k], f[m, k](r)))

where I already know the definition of f[0,0](r)=-(3/4)*r+1/(4*r)+(1/2)*r^2   and A[0,0]=1/(2*r^3)+1+(2*((3/4)*r-1/(4*r)-(1/2)*r^2))/r^2 

So, ode[1,1] can be solved with respect to f[1,1](r).

I would be glad for any comments. 

Here is my code.

 test.mw

Int(exp(t), t= 0 ..x) +1 is just a way to write exp(x), use value(%).

ff:= x -> evalf(Int(exp(t), t= 0 .. x, method = _d01ajc) + 1) does
the same, but in a numerical way.

Now consider that as a differential equation, to be solved numerically:

de:=D(f)(x) = f(x);
dsolve({%, f(0)=1});
                         f(x) = exp(x)


ff:= x -> evalf(Int(exp(t), t= 0 .. x, method = _d01ajc) + 1);

D(f)(x) = ff(x);
sol:=dsolve({%, f(0)=1}, numeric);

  Warning, The use of global variables in numerical ODE problems is deprecated,
  and will be removed in a future release. Use the 'parameters' argument instead
  (see ?dsolve,numeric,parameters)

sol(1);

  Error, (in sol) parameter 't' must be assigned a numeric value before
  obtaining a solution

Question: how to write it down correctly?

I have been here before...  My head is full of cotton, as usual.

f := proc (t) 2*t^3+9*t^2-60*t+1 end proc;
deq:=diff(f(t),t);
df := t->deq;  ## this is most likely one of my problems.
plot(df(t),t=-10..10);
isneg := x -> if is(df(x) < 0) then df(x) else 0 end if;  ## and, another??
plot(isneg(t),t=-10..10);
Plt1:=plot(isneg(t),t=-10..10);
plots:-display(Plt1);
isneg(0.12017234);

The if statement is not fully evauated.

I am missing something.  What?

Tom Dean

Hi dear friends

How do I improve my ability to manipulate code?

Thanks

_________________________________________________________________________________________

restart:

Digits := 30: m := 20: p0 := 8.: g := -0.3: nu := 0.3: a := 1:

w := sum(b[n]*r^(1.02*n), n = 1 .. m):

ODE := r^3*(g*r+a)^3*(diff(w, r, r, r, r))+2*(diff(w, r, r, r))*r^2*(4*g*r+a)*(g*r+a)^2+(diff(w, r, r))*r*(3*g^3*nu*r^3+P*a^3*r^2+6*a*g^2*nu*r^2+11*g^3*r^3+3*a^2*g*nu*r+15*a*g^2*r^2+3*a^2*g*r-a^3)+(diff(w, r))*(6*g^3*nu*r^3+P*a^3*r^2+6*a*g^2*nu*r^2-2*g^3*r^3-3*a*g^2*r^2+a^3):

for P from p0 by 0.1e-1 to p0+1 do
W := simplify(subs(solve({seq(evalf(subs(r = j/(m+2), ODE)), j = 2 .. m-1), subs(r = 1, diff(w, r))}, {seq(b[j], j = 1 .. m)}), w/b[m])):
F := (int(((diff(W, `$`(r, 2))+(diff(W, r))/r)^2-(2*(1-nu))*(diff(W, `$`(r, 2)))*(diff(W, r))/r)*r*(1+g*r/a)^3, r = 0 .. a))/(int((diff(W, r))^2*r, r = 0 .. a)):
delta[trunc(100*(P-p0))] := abs(F-P) end do:

plot(W, r = 0 .. a);

Min := min(seq(delta[j], j = 0 .. 100)):

for j from 0 to 100 do if delta[j] = Min then print(p0+j/(100.)) else  end if end do;

Hi all,

I have this matrix for example:

A := <<1,2>|<3,4>>;

How to determine: 2^A

Thank you