# Revision history [back]

### I've some problems with modelization of differential equations

Hi ! I've a project in mathematics and I'm supposed to modelize 5 functions according to various parameters that vary (hence the use of sliders). But I can't because it is highly developed functions and I had no lessons on this program. I give you the code which gives me this error:

File "<ipython-input-1-943112df05f1>", line 30 T.ode_solve( y_0=[lambda,Integer(0),Integer(1)-lambda,Integer(0)] , t_span=[Integer(0),Integer(150)] , params=[nu,d,B,v,c,p,beta] , num_points=Integer(1000)) ^ SyntaxError: invalid syntax

Code

# t : temps

@interact def interactive_function(nu = slider(0.01, 0.1, 0.005, default=0.03125), d = slider(0, 1, 0.05, default=1), B = slider(500, 30000, 250, default=13333), v = slider(0, 1, 0.05, default=0.2), c = slider(0, 50, 1, default=4), p = slider(0, 1, 0.05, default=0.3), beta = slider(0, 1, 0.002, default=0.014),):

def f_1(t,y,params) :
return [ -params-params*y-y*params*y , y*params*y-(params+params)*y, (1-params)*params*y-params*y, params*params*y-(params+params)*y, (params*y)/(y+y+y) ]

T = ode_solver()
T.function = f_1
T.algorithm="rk8pd"
T.ode_solve( y_0=[lambda,0,1-lambda,0] , t_span=[0,150] , params=[nu,d,B,v,c,p,beta] , num_points=1000 )

f = T.solution
X = [(x,x) for x in f]
Y = [(x,x) for x in f]
Z = [(x,x) for x in f]
A = [(x,x) for x in f]
lambda = [(x,x) for x in f]

P1 = line(X, rgbcolor='green')
P2 = line(Y, rgbcolor='pink')
P3 = line(Z, rgbcolor='red')
P4 = line(A, rgbcolor='brown')
P5 = line(lambda, rgbcolor='green')
show(P1+P2+P3+P4+P5)


Cordially, LordHorus.

PS = Sorry if my english can be wrong, I'm french :p 2 None

### I've some problems with modelization of differential equations

Hi ! I've a project in mathematics and I'm supposed to modelize 5 functions according to various parameters that vary (hence the use of sliders). But I can't because it is highly developed functions and I had no lessons on this program. I give you the code which gives me this error:

 File "<ipython-input-1-943112df05f1>", line 30
T.ode_solve( y_0=[lambda,Integer(0),Integer(1)-lambda,Integer(0)] , t_span=[Integer(0),Integer(150)] , params=[nu,d,B,v,c,p,beta] , num_points=Integer(1000))
^
SyntaxError: invalid syntaxsyntax


Code

# # y = X : – B - nuX - (lambda)cX

 (lambda)cX # y = Y : (lambda)cX - (v + nu)Y nu)Y # y = Z : (1 - p)vY - nuZ nuZ # y = A : pvY - (d + nu)A nu)A # y = lambda : ((beta)Y / (X + Y + Z)) Z)) # params = nu : natural death rate == 0.03125 0.03125 # params = d : AIDS death rate == 1 1 # params = B : immigration rate of people likely == 13333.3 13333.3 # params = v : conversion rate HIV-> AIDS == 0.2 0.2 # params = c : number of sexual partners == 4 4 # params = p : infectious HIV-positive proportion == 0.3 0.3 # params = beta : probability of transmission == 0.014 0.014 # t : temps temps @interact def interactive_function(nu = slider(0.01, 0.1, 0.005, default=0.03125), d = slider(0, 1, 0.05, default=1), B = slider(500, 30000, 250, default=13333), v = slider(0, 1, 0.05, default=0.2), c = slider(0, 50, 1, default=4), p = slider(0, 1, 0.05, default=0.3), beta = slider(0, 1, 0.002, default=0.014),): default=0.014),): def f_1(t,y,params) : return [ -params-params*y-y*params*y , y*params*y-(params+params)*y, (1-params)*params*y-params*y, params*params*y-(params+params)*y, (params*y)/(y+y+y) ] T = ode_solver() T.function = f_1 T.algorithm="rk8pd" T.ode_solve( y_0=[lambda,0,1-lambda,0] , t_span=[0,150] , params=[nu,d,B,v,c,p,beta] , num_points=1000 ) f = T.solution X = [(x,x) for x in f] Y = [(x,x) for x in f] Z = [(x,x) for x in f] A = [(x,x) for x in f] lambda = [(x,x) for x in f] P1 = line(X, rgbcolor='green') P2 = line(Y, rgbcolor='pink') P3 = line(Z, rgbcolor='red') P4 = line(A, rgbcolor='brown') P5 = line(lambda, rgbcolor='green') show(P1+P2+P3+P4+P5) Thanks for the help you can bring to me ! Cordially, LordHorus. PS = Sorry if my english can be wrong, I'm french :p Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license. about | faq | help | privacy policy | terms of service Powered by Askbot version 0.7.59 Please note: Askbot requires javascript to work properly, please enable javascript in your browser, here is how //IE fix to hide the red margin var noscript = document.getElementsByTagName('noscript'); noscript.style.padding = '0px'; noscript.style.backgroundColor = 'transparent'; askbot['urls']['mark_read_message'] = '/s/messages/markread/'; askbot['urls']['get_tags_by_wildcard'] = '/s/get-tags-by-wildcard/'; askbot['urls']['get_tag_list'] = '/s/get-tag-list/'; askbot['urls']['follow_user'] = '/followit/follow/user/{{userId}}/'; askbot['urls']['unfollow_user'] = '/followit/unfollow/user/{{userId}}/'; askbot['urls']['user_signin'] = '/account/signin/'; askbot['urls']['getEditor'] = '/s/get-editor/'; askbot['urls']['apiGetQuestions'] = '/s/api/get_questions/'; askbot['urls']['ask'] = '/questions/ask/'; askbot['urls']['questions'] = '/questions/'; askbot['settings']['groupsEnabled'] = false; askbot['settings']['static_url'] = '/m/'; askbot['settings']['minSearchWordLength'] = 4; askbot['settings']['mathjaxEnabled'] = true; askbot['settings']['sharingSuffixText'] = ''; askbot['settings']['errorPlacement'] = 'after-label'; askbot['data']['maxCommentLength'] = 800; askbot['settings']['editorType'] = 'markdown'; askbot['settings']['commentsEditorType'] = 'rich\u002Dtext'; askbot['messages']['askYourQuestion'] = 'Ask Your Question'; askbot['messages']['acceptOwnAnswer'] = 'accept or unaccept your own answer'; askbot['messages']['followQuestions'] = 'follow questions'; askbot['settings']['allowedUploadFileTypes'] = [ "jpg", "jpeg", "gif", "bmp", "png", "tiff" ]; askbot['data']['haveFlashNotifications'] = true; askbot['data']['activeTab'] = 'questions'; askbot['settings']['csrfCookieName'] = 'asksage_csrf'; askbot['data']['searchUrl'] = ''; /*<![CDATA[*/ $('.mceStatusbar').remove();//a hack to remove the tinyMCE status bar$(document).ready(function(){ // focus input on the search bar endcomment var activeTab = askbot['data']['activeTab']; if (inArray(activeTab, ['users', 'questions', 'tags', 'badges'])) { var searchInput = $('#keywords'); } else if (activeTab === 'ask') { var searchInput =$('#id_title'); } else { var searchInput = undefined; animateHashes(); } var wasScrolled = $('#scroll-mem').val(); if (searchInput && !wasScrolled) { searchInput.focus(); putCursorAtEnd(searchInput); } var haveFullTextSearchTab = inArray(activeTab, ['questions', 'badges', 'ask']); var haveUserProfilePage =$('body').hasClass('user-profile-page'); if ((haveUserProfilePage || haveFullTextSearchTab) && searchInput && searchInput.length) { var search = new FullTextSearch(); askbot['controllers'] = askbot['controllers'] || {}; askbot['controllers']['fullTextSearch'] = search; search.setSearchUrl(askbot['data']['searchUrl']); if (activeTab === 'ask') { search.setAskButtonEnabled(false); } search.decorate(searchInput); } else if (activeTab === 'tags') { var search = new TagSearch(); search.decorate(searchInput); } if (askbot['data']['userIsAdminOrMod']) { $('body').addClass('admin'); } if (askbot['settings']['groupsEnabled']) { askbot['urls']['add_group'] = "/s/add-group/"; var group_dropdown = new GroupDropdown();$('.groups-dropdown').append(group_dropdown.getElement()); } var userRep = $('#userToolsNav .reputation'); if (userRep.length) { var showPermsTrigger = new ShowPermsTrigger(); showPermsTrigger.decorate(userRep); } }); if (askbot['data']['haveFlashNotifications']) {$('#validate_email_alert').click(function(){notify.close(true)}) notify.show(); } var langNav = $('.lang-nav'); if (langNav.length) { var nav = new LangNav(); nav.decorate(langNav); } /*]]>*/ if (typeof MathJax != 'undefined') { MathJax.Hub.Config({ extensions: ["tex2jax.js"], jax: ["input/TeX","output/HTML-CSS"], tex2jax: {inlineMath: [["$","$"],["\$","\$"]]} }); } else { console.log('Could not load MathJax'); } //todo - take this out into .js file$(document).ready(function(){ $('div.revision div[id^=rev-header-]').bind('click', function(){ var revId = this.id.substr(11); toggleRev(revId); }); lanai.highlightSyntax(); }); function toggleRev(id) { var arrow =$("#rev-arrow-" + id); var visible = arrow.attr("src").indexOf("hide") > -1; if (visible) { var image_path = '/m/default/media/images/expander-arrow-show.gif?v=19'; } else { var image_path = '/m/default/media/images/expander-arrow-hide.gif?v=19'; } image_path = image_path + "?v=19"; arrow.attr("src", image_path); \$("#rev-body-" + id).slideToggle("fast"); } for (url_name in askbot['urls']){ askbot['urls'][url_name] = cleanUrl(askbot['urls'][url_name]); }