# Revision history [back]

### graphing ODE's using euler's method

Im trying to plot a graph based on the change in susceptible, exposed and infected people due to guinea worm, included is a plot of the worm population W. This is an attempt to at least replicate some of what is from this mathematical model

(http://mysite.science.uottawa.ca/rsmith43/GuineaWorm.pdf ),

hopefully using the Impulsive differential equations but right now I've been relentlessly trying to graph these equations on a reasonable axis over ten years but for some reason it doesn't. Plus the y axis is completely off from what it should be, does anyone know what is wrong with the code? Thanks for any responses

timedata=[]

Sdata=[]

Edata=[]

Idata=[]

Wdata=[]

t=0

S=2600

E=1

I=1

W=200

bR= 37

k= 8760

m= .0142

r= .9

B= .0255*(1/W)

A= 1

Y=100000

mW=26

dt=.1

timedata.append(t)

Sdata.append(S)

Edata.append(E)

Idata.append(I)

Wdata.append(W)

T=10

Tfinal=(T/dt)

for i in range(0,Tfinal):

t= t+dt

Sprime=bR-B*S*W-m*S+k*I

Eprime=B*S*W-A*E-m*E

Iprime=A*E-k*I-m*I

Wprime=Y*I-mW*W

S= S+(Sprime*dt)

E= E+(Eprime*dt)

I= I+(Iprime*dt)

W= W+(Wprime*dt)

timedata.append(t)

Sdata.append(S)

Edata.append(E)

Idata.append(I)

Wdata.append(W)


Splot=list_plot(zip(timedata,Sdata),color='green',plotjoined=True)

Eplot=list_plot(zip(timedata,Edata),color='orange',plotjoined=True)

Iplot=list_plot(zip(timedata,Idata),color='black',plotjoined=True)

Splot

 2 No.2 Revision kcrisman 12072 ●38 ●126 ●246

### graphing ODE's using euler's method

Im trying to plot a graph based on the change in susceptible, exposed and infected people due to guinea worm, included is a plot of the worm population W. This is an attempt to at least replicate some of what is from this mathematical model

(http://mysite.science.uottawa.ca/rsmith43/GuineaWorm.pdf ),

hopefully using the Impulsive differential equations but right now I've been relentlessly trying to graph these equations on a reasonable axis over ten years but for some reason it doesn't. Plus the y axis is completely off from what it should be, does anyone know what is wrong with the code? Thanks for any responses

timedata=[]

Sdata=[]

Edata=[]

Idata=[]

Wdata=[]

t=0

S=2600

E=1

I=1

W=200

timedata=[]
Sdata=[]
Edata=[]
Idata=[]
Wdata=[]

t=0
S=2600
E=1
I=1
W=200
bR= 37 37
k= 8760   m= .0142   r= .9            B= .0255*(1/W) .0255*(1/W)
A= 1 Y=100000 mW=26 dt=.1 timedata.append(t) Sdata.append(S) Edata.append(E) Idata.append(I) Wdata.append(W) T=10 Tfinal=(T/dt) 1
Y=100000
mW=26
dt=.1

timedata.append(t)
Sdata.append(S)
Edata.append(E)
Idata.append(I)
Wdata.append(W)

T=10
Tfinal=(T/dt)
for i in range(0,Tfinal): range(0,Tfinal):
t= t+dt
Sprime=bR-B*S*W-m*S+k*I
Eprime=B*S*W-A*E-m*E
Iprime=A*E-k*I-m*I
Wprime=Y*I-mW*W
S= S+(Sprime*dt)
E= E+(Eprime*dt)
I= I+(Iprime*dt)
W= W+(Wprime*dt)
timedata.append(t)
Sdata.append(S)
Edata.append(E)
Idata.append(I)
Wdata.append(W)

Splot=list_plot(zip(timedata,Sdata),color='green',plotjoined=True)
Eplot=list_plot(zip(timedata,Edata),color='orange',plotjoined=True)
Iplot=list_plot(zip(timedata,Idata),color='black',plotjoined=True)

Splot
Splot=list_plot(zip(timedata,Sdata),color='green',plotjoined=True) Eplot=list_plot(zip(timedata,Edata),color='orange',plotjoined=True) Iplot=list_plot(zip(timedata,Idata),color='black',plotjoined=True) Splot 


 3 retagged updated 2015-01-13 18:20:17 +0200 FrédéricC 4131 ●3 ●37 ●85 graphing ODE's using euler's method Im trying to plot a graph based on the change in susceptible, exposed and infected people due to guinea worm, included is a plot of the worm population W. This is an attempt to at least replicate some of what is from this mathematical model (http://mysite.science.uottawa.ca/rsmith43/GuineaWorm.pdf ), hopefully using the Impulsive differential equations but right now I've been relentlessly trying to graph these equations on a reasonable axis over ten years but for some reason it doesn't. Plus the y axis is completely off from what it should be, does anyone know what is wrong with the code? Thanks for any responses timedata=[] Sdata=[] Edata=[] Idata=[] Wdata=[] t=0 S=2600 E=1 I=1 W=200 bR= 37 k= 8760 m= .0142 r= .9 B= .0255*(1/W) A= 1 Y=100000 mW=26 dt=.1 timedata.append(t) Sdata.append(S) Edata.append(E) Idata.append(I) Wdata.append(W) T=10 Tfinal=(T/dt) for i in range(0,Tfinal): t= t+dt Sprime=bR-B*S*W-m*S+k*I Eprime=B*S*W-A*E-m*E Iprime=A*E-k*I-m*I Wprime=Y*I-mW*W S= S+(Sprime*dt) E= E+(Eprime*dt) I= I+(Iprime*dt) W= W+(Wprime*dt) timedata.append(t) Sdata.append(S) Edata.append(E) Idata.append(I) Wdata.append(W) Splot=list_plot(zip(timedata,Sdata),color='green',plotjoined=True) Eplot=list_plot(zip(timedata,Edata),color='orange',plotjoined=True) Iplot=list_plot(zip(timedata,Idata),color='black',plotjoined=True) Splot 


 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')[0]; 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]); }