function y = vfinitcond()

y = [
 	0.02      % DV1                         1   
	0.02      % DV2                         2 
	0.5     % PPV1                           3
	0.5     % PPV2                           4
	0.0    % Contractile1                  5
	0.0    % Contractile2                  6
	0.0     % d(theta)1                      7
	0.0     % theta2                         8
    0.0     % Renshaw 1                     9
    0.0     % Renshaw 2                     10
	1.5     % alpha-MN1                      11
	1.5     % alpha-MN2                      12
	0.0    % Static gamma-MN 1              13
    0.0   % Static gamma-MN 2               14
    0.0    % Intrafusal static gamma1      15
	0.0   % Intrafusal static gamma2      16
	0.0    % Dynamic gamma-MN 1             17
	0.0    % Dynamic gamma-MN 2             18
	0.0    % Intrafusal dynamic gamma1     19
	0.0    % Intrafusal dynamic gamma2     20
    0.0    % IaIN 1                         21
    0.0     % IaIN 2                        22
    0.02    %  Spindle response 1            23
    0.02     % Spindle response 2            24
    0.0     % IbIN 1                        25
    0.0     % IbIN 2                        26
  ];