r - Character limit for ODE solver -


i'm trying solve odes in r.

currently i'm using package desolve, running problems there appears character limit on length of equations, error tells me there's unexpected symbol.

initially solved storing half of offending equation in vector combining... ideally not every equation.

does know way around character limit package, or of ode solver in r without character limit?

as requested, here equation causes error unexpected symbol:

dyx31 <- -c_2*yx14*(-6*c_6^2*y_0 + 4*c_6^2 - 12*c_6*y_0^2 + 4*c_6*y_0 + 4*c_6*yx16 + c_6 - 6*y_0^3)/(c_6^3 + 3*c_6^2*y_0 + 3*c_6*y_0^2 + y_0^3) + c_2*yx29*(-4*c_6^3*y_0 + 6*c_6^3 - 12*c_6^2*y_0^2 + 12*c_6^2*y_0 + 4*c_6^2 - 12*c_6*y_0^3 + 6*c_6*y_0^2 + 4*c_6*y_0 + 4*c_6*yx16 + c_6 - 4*y_0^4)/(c_6^4 + 4*c_6^3*y_0 + 6*c_6^2*y_0^2 + 4*c_6*y_0^3 + y_0^4) - c_2*yx4*(4*c_6*y_0 - 4*c_6*yx16 - c_6 + 4*y_0^2)/(c_6^2 + 2*c_6*y_0 + y_0^2) + yx16*(4*c_0*c_6^4 + 16*c_0*c_6^3*y_0 + 24*c_0*c_6^2*y_0^2 + 16*c_0*c_6*y_0^3 + 4*c_0*y_0^4 - 4*c_1*c_6^4*y_0 + 6*c_1*c_6^4 - 16*c_1*c_6^3*y_0^2 + 24*c_1*c_6^3*y_0 - 24*c_1*c_6^2*y_0^3 + 36*c_1*c_6^2*y_0^2 - 16*c_1*c_6*y_0^4 + 24*c_1*c_6*y_0^3 - 4*c_1*y_0^5 + 6*c_1*y_0^4 - 4*c_2*c_6^3*y_0*y_2 + 6*c_2*c_6^3*y_2 - 12*c_2*c_6^2*y_0^2*y_2 + 12*c_2*c_6^2*y_0*y_2 + 4*c_2*c_6^2*y_2 - 12*c_2*c_6*y_0^3*y_2 + 6*c_2*c_6*y_0^2*y_2 + 4*c_2*c_6*y_0*y_2 + 4*c_2*c_6*y_2*yx16 + c_2*c_6*y_2 - 4*c_2*y_0^4*y_2)/(c_6^4 + 4*c_6^3*y_0 + 6*c_6^2*y_0^2 + 4*c_6*y_0^3 + y_0^4) - yx31*(4*c_1*c_6^5 + 20*c_1*c_6^4*y_0 + 40*c_1*c_6^3*y_0^2 + 40*c_1*c_6^2*y_0^3 + 20*c_1*c_6*y_0^4 + 4*c_1*y_0^5 + 4*c_2*c_6^4*y_2 + 12*c_2*c_6^3*y_0*y_2 + 6*c_2*c_6^3*y_2 + 12*c_2*c_6^2*y_0^2*y_2 + 12*c_2*c_6^2*y_0*y_2 + 4*c_2*c_6^2*y_2 + 4*c_2*c_6*y_0^3*y_2 + 6*c_2*c_6*y_0^2*y_2 + 4*c_2*c_6*y_0*y_2 + 4*c_2*c_6*y_2*yx16 + c_2*c_6*y_2)/(c_6^5 + 5*c_6^4*y_0 + 10*c_6^3*y_0^2 + 10*c_6^2*y_0^3 + 5*c_6*y_0^4 + y_0^5) + yx6*(6*c_0*c_6^3 + 18*c_0*c_6^2*y_0 + 18*c_0*c_6*y_0^2 + 6*c_0*y_0^3 + 6*c_1*c_6^3*y_0 - 4*c_1*c_6^3 + 18*c_1*c_6^2*y_0^2 - 12*c_1*c_6^2*y_0 + 18*c_1*c_6*y_0^3 - 12*c_1*c_6*y_0^2 + 6*c_1*y_0^4 - 4*c_1*y_0^3 + 6*c_2*c_6^2*y_0*y_2 - 4*c_2*c_6^2*y_2 + 12*c_2*c_6*y_0^2*y_2 - 4*c_2*c_6*y_0*y_2 - 4*c_2*c_6*y_2*yx16 - c_2*c_6*y_2 + 6*c_2*y_0^3*y_2)/(c_6^3 + 3*c_6^2*y_0 + 3*c_6*y_0^2 + y_0^3) + (-4*c_0*c_6*yx16 + c_0*c_6 - 4*c_0*y_0*yx16 + c_0*y_0 + 4*c_1*c_6*y_0*yx16 + c_1*c_6*y_0 + 4*c_1*y_0^2*yx16 + c_1*y_0^2 + 4*c_2*y_0*y_2*yx16 + c_2*y_0*y_2)/(c_6 + y_0)

and here temporary solution problem 

temp <- 4*c_2*c_6*y_0^2*y_2*yx31 + 4*c_2*c_6*y_0^2*yx16*yx29 + c_2*c_6*y_0^2*yx29 - 4*c_2*c_6*y_0*y_2*yx16*yx31 - c_2*c_6*y_0*y_2*yx31 + 4*c_2*c_6*y_2*yx6^3*(20*y_0^2 + 24*yx6) + 6*c_2*c_6*y_2*yx6^3*(20*y_0^2 + 24*yx6)/y_0 + 4*c_2*c_6*y_2*yx6^3*(20*y_0^2 + 24*yx6)/y_0^2 + 4*c_2*c_6*y_2*yx16*yx6^3*(20*y_0^2 + 24*yx6)/y_0^3 + c_2*c_6*y_2*yx6^3*(20*y_0^2 + 24*yx6)/y_0^3 + 6*c_2*y_0^6*y_2*yx6 + c_2*y_0^6*y_2 + 6*c_2*y_0^6*yx14 - 4*c_2*y_0^6*yx29 - 4*c_2*y_0^6*yx4  dyx31 <- (6*c_0*c_6^6*yx6 + c_0*c_6^6 + 36*c_0*c_6^5*y_0*yx6 + 6*c_0*c_6^5*y_0 + 90*c_0*c_6^4*y_0^2*yx6 + 15*c_0*c_6^4*y_0^2 + 120*c_0*c_6^3*y_0^3*yx6 + 20*c_0*c_6^3*y_0^3 + 90*c_0*c_6^2*y_0^4*yx6 + 15*c_0*c_6^2*y_0^4 + 36*c_0*c_6*y_0^5*yx6 + 6*c_0*c_6*y_0^5 + 6*c_0*y_0^6*yx6 + c_0*y_0^6 + 6*c_1*c_6^6*y_0*yx6 + c_1*c_6^6*y_0 + 6*c_1*c_6^6*yx16 - 4*c_1*c_6^6*yx31 - 4*c_1*c_6^6*yx6 + 36*c_1*c_6^5*y_0^2*yx6 + 6*c_1*c_6^5*y_0^2 + 36*c_1*c_6^5*y_0*yx16 - 24*c_1*c_6^5*y_0*yx31 - 24*c_1*c_6^5*y_0*yx6 + 90*c_1*c_6^4*y_0^3*yx6 + 15*c_1*c_6^4*y_0^3 + 90*c_1*c_6^4*y_0^2*yx16 - 60*c_1*c_6^4*y_0^2*yx31 - 60*c_1*c_6^4*y_0^2*yx6 + 120*c_1*c_6^3*y_0^4*yx6 + 20*c_1*c_6^3*y_0^4 + 120*c_1*c_6^3*y_0^3*yx16 - 80*c_1*c_6^3*y_0^3*yx31 - 80*c_1*c_6^3*y_0^3*yx6 + 90*c_1*c_6^2*y_0^5*yx6 + 15*c_1*c_6^2*y_0^5 + 90*c_1*c_6^2*y_0^4*yx16 - 60*c_1*c_6^2*y_0^4*yx31 - 60*c_1*c_6^2*y_0^4*yx6 + 36*c_1*c_6*y_0^6*yx6 + 6*c_1*c_6*y_0^6 + 36*c_1*c_6*y_0^5*yx16 - 24*c_1*c_6*y_0^5*yx31 - 24*c_1*c_6*y_0^5*yx6 + 6*c_1*y_0^7*yx6 + c_1*y_0^7 + 6*c_1*y_0^6*yx16 - 4*c_1*y_0^6*yx31 - 4*c_1*y_0^6*yx6 + 6*c_2*c_6^5*y_0*y_2*yx6 + c_2*c_6^5*y_0*y_2 + 6*c_2*c_6^5*y_0*yx14 - 4*c_2*c_6^5*y_0*yx29 - 4*c_2*c_6^5*y_0*yx4 + 6*c_2*c_6^5*y_2*yx16 - 4*c_2*c_6^5*y_2*yx31 - 4*c_2*c_6^5*y_2*yx6 - 4*c_2*c_6^5*yx14 + 4*c_2*c_6^5*yx16*yx4 + 6*c_2*c_6^5*yx29 + c_2*c_6^5*yx4 + 30*c_2*c_6^4*y_0^2*y_2*yx6 + 5*c_2*c_6^4*y_0^2*y_2 + 30*c_2*c_6^4*y_0^2*yx14 - 20*c_2*c_6^4*y_0^2*yx29 - 20*c_2*c_6^4*y_0^2*yx4 + 24*c_2*c_6^4*y_0*y_2*yx16 - 16*c_2*c_6^4*y_0*y_2*yx31 - 16*c_2*c_6^4*y_0*y_2*yx6 - 16*c_2*c_6^4*y_0*yx14 + 16*c_2*c_6^4*y_0*yx16*yx4 + 24*c_2*c_6^4*y_0*yx29 + 4*c_2*c_6^4*y_0*yx4 - 4*c_2*c_6^4*y_2*yx16*yx6 + 4*c_2*c_6^4*y_2*yx16 - 6*c_2*c_6^4*y_2*yx31 - c_2*c_6^4*y_2*yx6 - 4*c_2*c_6^4*yx14*yx16 - c_2*c_6^4*yx14 + 4*c_2*c_6^4*yx29 + 4*c_2*c_6^4*y_2*yx6^3*(20*y_0^2 + 24*yx6)/y_0^3 + 60*c_2*c_6^3*y_0^3*y_2*yx6 + 10*c_2*c_6^3*y_0^3*y_2 + 60*c_2*c_6^3*y_0^3*yx14 - 40*c_2*c_6^3*y_0^3*yx29 - 40*c_2*c_6^3*y_0^3*yx4 + 36*c_2*c_6^3*y_0^2*y_2*yx16 - 24*c_2*c_6^3*y_0^2*y_2*yx31 - 24*c_2*c_6^3*y_0^2*y_2*yx6 - 24*c_2*c_6^3*y_0^2*yx14 + 24*c_2*c_6^3*y_0^2*yx16*yx4 + 36*c_2*c_6^3*y_0^2*yx29 + 6*c_2*c_6^3*y_0^2*yx4 - 12*c_2*c_6^3*y_0*y_2*yx16*yx6 + 12*c_2*c_6^3*y_0*y_2*yx16 - 18*c_2*c_6^3*y_0*y_2*yx31 - 3*c_2*c_6^3*y_0*y_2*yx6 - 12*c_2*c_6^3*y_0*yx14*yx16 - 3*c_2*c_6^3*y_0*yx14 + 12*c_2*c_6^3*y_0*yx29 + 4*c_2*c_6^3*y_2*yx16^2 + c_2*c_6^3*y_2*yx16 - 4*c_2*c_6^3*y_2*yx31 + 4*c_2*c_6^3*yx16*yx29 + c_2*c_6^3*yx29 + 12*c_2*c_6^3*y_2*yx6^3*(20*y_0^2 + 24*yx6)/y_0^2 + 6*c_2*c_6^3*y_2*yx6^3*(20*y_0^2 + 24*yx6)/y_0^3 + 60*c_2*c_6^2*y_0^4*y_2*yx6 + 10*c_2*c_6^2*y_0^4*y_2 + 60*c_2*c_6^2*y_0^4*yx14 - 40*c_2*c_6^2*y_0^4*yx29 - 40*c_2*c_6^2*y_0^4*yx4 + 24*c_2*c_6^2*y_0^3*y_2*yx16 - 16*c_2*c_6^2*y_0^3*y_2*yx31- 16*c_2*c_6^2*y_0^3*y_2*yx6 - 16*c_2*c_6^2*y_0^3*yx14 + 16*c_2*c_6^2*y_0^3*yx16*yx4 + 24*c_2*c_6^2*y_0^3*yx29 + 4*c_2*c_6^2*y_0^3*yx4 - 12*c_2*c_6^2*y_0^2*y_2*yx16*yx6 + 12*c_2*c_6^2*y_0^2*y_2*yx16 - 18*c_2*c_6^2*y_0^2*y_2*yx31 - 3*c_2*c_6^2*y_0^2*y_2*yx6 - 12*c_2*c_6^2*y_0^2*yx14*yx16 - 3*c_2*c_6^2*y_0^2*yx14 + 12*c_2*c_6^2*y_0^2*yx29 + 8*c_2*c_6^2*y_0*y_2*yx16^2 + 2*c_2*c_6^2*y_0*y_2*yx16 - 8*c_2*c_6^2*y_0*y_2*yx31 + 8*c_2*c_6^2*y_0*yx16*yx29 + 2*c_2*c_6^2*y_0*yx29 - 4*c_2*c_6^2*y_2*yx16*yx31 - c_2*c_6^2*y_2*yx31 + 12*c_2*c_6^2*y_2*yx6^3*(20*y_0^2 + 24*yx6)/y_0 + 12*c_2*c_6^2*y_2*yx6^3*(20*y_0^2 + 24*yx6)/y_0^2 + 4*c_2*c_6^2*y_2*yx6^3*(20*y_0^2 + 24*yx6)/y_0^3 + 30*c_2*c_6*y_0^5*y_2*yx6 + 5*c_2*c_6*y_0^5*y_2 + 30*c_2*c_6*y_0^5*yx14 - 20*c_2*c_6*y_0^5*yx29 - 20*c_2*c_6*y_0^5*yx4 + 6*c_2*c_6*y_0^4*y_2*yx16 - 4*c_2*c_6*y_0^4*y_2*yx31 - 4*c_2*c_6*y_0^4*y_2*yx6 - 4*c_2*c_6*y_0^4*yx14 + 4*c_2*c_6*y_0^4*yx16*yx4 + 6*c_2*c_6*y_0^4*yx29 + c_2*c_6*y_0^4*yx4 - 4*c_2*c_6*y_0^3*y_2*yx16*yx6 + 4*c_2*c_6*y_0^3*y_2*yx16 - 6*c_2*c_6*y_0^3*y_2*yx31 - c_2*c_6*y_0^3*y_2*yx6 - 4*c_2*c_6*y_0^3*yx14*yx16 - c_2*c_6*y_0^3*yx14 + 4*c_2*c_6*y_0^3*yx29 + 4*c_2*c_6*y_0^2*y_2*yx16^2 + c_2*c_6*y_0^2*y_2*yx16   - temp)/(c_6^6 + 6*c_6^5*y_0 + 15*c_6^4*y_0^2 + 20*c_6^3*y_0^3 + 15*c_6^2*y_0^4 + 6*c_6*y_0^5 + y_0^6)

they longer , messier, avoid having repeatedly.

at moment i've resorted using program called cvode i'd able in r if possible.


Comments

Post a Comment

Popular posts from this blog

html - How to style widget with post count different than without post count -

i have border below li elements in wordpress widget types appropriate. element inside set block gives them nice easy click full width. the problem, when choose option "show post count" when adding widget such categories or archives displays post count plain text after tag , because tag block pushes post number text line below. whether choose show post count or not show post count still gives exact same element identifyer example "widget_categories". how style post count not on next line still keep links non post count list items full width? here's example on sidebar www.bbmthemes.com/themes/modular/blog-standard/ in order solve problem need remove display:block; for ul.widgets ul a , add in it's place line-height: 33px; what going achieve have same visual effect. what going lose doing links going work on text , not on "box". the other way out of find post count coming , make changes this. example change text count i...
Read more

How to remove text and logo OR add Overflow on Android ActionBar using AppCompat on API 8? -

have researched death , cannot find answer. i have actionbar using appcompat. supporting api v7 (api v8 do). actionbar works api v11. api < v11 has problem of removing icon (and feature) actionbar without supplying overflow. big picture have app logo, app title, , 3 icons squidging same actionbar on low end phone. trying make room! i happy 1 of 2 solutions. either: a method of getting overflow functionality can retrieved pull down. (preferred) method of removing icon , text actionbar below current xml api 11+: <resources> <!-- base application theme api 11+. theme replaces appbasetheme res/values/styles.xml on api 11+ devices. --> <style name="appbasetheme" parent="theme.appcompat.light.darkactionbar"> <item name="android:actionbarstyle">@style/liteactionbarstyle</item> <item name="actionbarstyle">@style/liteactionbarstyle</item> </st...
Read more

IIS->Tomcat Redirect: multiple worker with default -

i trying configure multiple instances of tomcat (7) behind single domain in iis. tomcat instances both run on separate ports. have installed tomcat isapi_rediect.dll , configured iis. everything works redirector, except (you knew coming) being told "all other traffic must mapped 1 of 2 instances (they different applications). more accurate requests except specific pattern need mapped 1 worker , specific pattern (/rpt) needs go other worker. i have 2 workers defined in workers.properties: worker.list=ajp13_1,ajp13_2 worker.ajp13_1.port=8009 worker.ajp13_2.port=8010 and following in uriworkermap.properties: /wxx|/*=ajp13_1 /rpt|/*=ajp13_2 with this, seems fine. except (considerable) number of request not include wxx scope. need map rpt , sub content ajp3_2 , other requests ajp3_1. if add /|*=ajp13_1 then ajp3_2 becomes unreachable , wildchar match maps /rpt/* uri ajp3_1. is there way trying do? appreciate help- have spent time reading same references ...
Read more