Body: <p>This is a computational neuroscience project. We are asked to plot the Morris-Lescar model (question 3 below).</p>
  <p>The model should be plotted with scilab or R.</p>
  <p>however, I don't know where to start? For the moment, the problem is more in writing the code to represent the equations and systems of equations than in actually solving them.</p>
  <p>Any clues as to how to start writing our code would be a huge start!</p>
  <p><a href="https://i.sstatic.net/toQtT.jpg" rel="nofollow noreferrer"><img src="https://i.sstatic.net/toQtT.jpg" alt="Morris-Lecar P1" /></a>
  <a href="https://i.sstatic.net/itO9R.jpg" rel="nofollow noreferrer"><img src="https://i.sstatic.net/itO9R.jpg" alt="Morris-Lecar P2" /></a></p>
  <p>In concrete terms, we don't know where to start with my working group, even for question 3.a..</p>
  <p>We know that nullclines are associated with <span class="math-container">$dV/dt = 0$</span> and <span class="math-container">$dW/dt = 0$</span>, which gives the following entries:</p>
  <p><span class="math-container">$W = \frac{I-g{Ca}M{\infty}(V)(V-V_{Ca})-g_L(V-V_L)}{g_K(V-VK)}$</span></p>
  <p><span class="math-container">$W = W{\infty}(V)$</span></p>
  <p>But here for example, no idea how to write it on scilab, or how to display them graphically by varying I.</p>
  <p>I hope the application is now clearer.</p>
  <p>Appendix: scilab code.</p>
  <pre><code>E_L = -50;
  E_Ca = 100;
  E_K = -70;
  V1 = -1.2;
  V2 = 18;
  V3 = 2;
  V4 = 30;
  g_L = 2;
  g_Ca = 4.4;
  g_K = 8.0;
  C = 20;
  gamma = 0.04;
  
  minf(v)=.5(1+tanh((v-va)/vb))
  
  winf(v)=.5(1+tanh((v-vc)/vd))
  
  tauw(v)=1/cosh((v-vc)/(2*vd))
  </code></pre>
  
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2022-04-11T12:34:57.063Z
Comments.Text: Please clarify your specific problem or provide additional details to highlight exactly what you need. As it's currently written, it's hard to tell exactly what you're asking.
Comments.UserDisplayName: M__
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2022-04-11T17:47:38.260Z
Comments.Text: I tried to make the question more to the point, but I think would be easier to answer if you would be even more specific.
Comments.UserDisplayName: Kamil S Jaron
Tags: <r><homework>
Title: Problem of writing scilab code on a project - Morris-Lecar model
ThreadText: # Problem of writing scilab code on a project - Morris-Lecar model
  
  (Asked by: username_0 on 2022-04-11)
  
  This is a computational neuroscience project. We are asked to plot the Morris-Lescar model (question 3 below).
  
  The model should be plotted with scilab or R.
  
  however, I don't know where to start? For the moment, the problem is more in writing the code to represent the equations and systems of equations than in actually solving them.
  
  Any clues as to how to start writing our code would be a huge start!
  
  [![Morris-Lecar P1](https://i.sstatic.net/toQtT.jpg)](https://i.sstatic.net/toQtT.jpg)
  [![Morris-Lecar P2](https://i.sstatic.net/itO9R.jpg)](https://i.sstatic.net/itO9R.jpg)
  
  In concrete terms, we don't know where to start with my working group, even for question 3.a..
  
  We know that nullclines are associated with $dV/dt = 0$ and $dW/dt = 0$, which gives the following entries:
  
  $W = \frac{I-g{Ca}M{\infty}(V)(V-V\_{Ca})-g\_L(V-V\_L)}{g\_K(V-VK)}$
  
  $W = W{\infty}(V)$
  
  But here for example, no idea how to write it on scilab, or how to display them graphically by varying I.
  
  I hope the application is now clearer.
  
  Appendix: scilab code.
  
  ```
  E_L = -50;
  E_Ca = 100;
  E_K = -70;
  V1 = -1.2;
  V2 = 18;
  V3 = 2;
  V4 = 30;
  g_L = 2;
  g_Ca = 4.4;
  g_K = 8.0;
  C = 20;
  gamma = 0.04;
  
  minf(v)=.5(1+tanh((v-va)/vb))
  
  winf(v)=.5(1+tanh((v-vc)/vd))
  
  tauw(v)=1/cosh((v-vc)/(2*vd))
  
  ```
  
  **Tags:** r, homework
  
  > **Comments:**
  > **username_1**: Please clarify your specific problem or provide additional details to highlight exactly what you need. As it's currently written, it's hard to tell exactly what you're asking.  
  > **username_2**: I tried to make the question more to the point, but I think would be easier to answer if you would be even more specific.  
  
  ---
  
  ## 0 answers
  
  