MATLAB HINTS - Code Skeleton
M. Nelson  Spring 2004

Here's a skeleton for your simulation code, you need to flesh it out with the details.
(Save it as hw1p2.m, or something similar)
Note: this code won't execute until you replace the *** SECTIONS *** with valid matlab code.
%----------------------------------------------------------------------
% PHYSL 317 / BIOPH 317 / NEURO 317 / BIOENG 376
% Methods in Computational Neurobiology - Fall 2004
%
% hw1p2.m - Homework #1, Problem #2
% RC circuit simulation
% Simulate step responses for a family of RC values (different values of R).
%
%  ***YOUR NAME***
%  *** DATE***
%----------------------------------------------------------------------

% Initialization
clear;
figure(1);
clf;

% Define constants
DT = 1;
TMAX = 100;
TIME  = DT : DT : TMAX;

% resistance and capacitance values                                   
RVALUE = [ 2, 5, 10, 20, 50];
C = 1; 

% Stimulus parameters
IMAX = 1.0;
TSTART = 10;
TSTOP = 60;

% loop over different R values
for ival = 1:length(RVALUE)

    R = RVALUE(ival) ;

    % set initial condition
    v(1) = 0; 

    % now loop over time points
    for itime = 2 : length(TIME)
        t = TIME(itime);

        *** compute the value of the input current I at this time-step based on TSTART, TSTOP***

        *** compute dv/dt , which depends on v(itime-1), I, R and C***

        % now update the voltage using Euler integration
        v(itime) = v(itime-1) + dvdt * DT;
    end
   
    % plot voltage versus time for this R value
    *** decide what plot symbol to use ***
    *** to make things simple, you can leave off the plot symbol  to get started ***
    plot(TIME, v, ***PLOTSYMBOL***)
    hold on;

    *** add a text label somewhere near the voltage trace ***
    *** you need to compute an X and Y LOCATION to draw the string ***
    myString = ['R=' num2str(R)];
    xloc = ***???***;
    yloc = *** ??? ***;
    text( xloc, yloc, myString);

end

*** set axis limits for the figure***
*** add title, xlabel, ylabel ***

GOOD LUCK!!!