unction varargout=indmot_ode1(t,y,flag)
switch flag
case ''
%no input flag
varargout{1}=FF1(t,y);
case 'mass'
%flag of mass calls mass.m
varargout{1}=MM1(t,y);
otherwise
error(['unknown flag ''' flag '''.']);
end


    function yp = FF1(t, y)
                dy = zeros(4,1);%-RK+FD+FE
                dy(1) = -k2v*z(1) -0*z(2) +k2v*z(3) -0*z(4)  -d2v(z'2v-z'1v)+0;
                dy(2) =  0*z(1) - k2h*z(2) - 0*z(3) +k2h*z(4) -d2h(z'2h-z'1h)+0;
                dy(3) =  +k2v*z(1) -0*z(2) -k2v*z(3) -0*z(4)  -d2v(z'1v-z'2v)-d1v(z'1v-h'v )+k1v(z1v-hv);
                dy(4) =  0*z(1) + k2h*z(2) - 0*z(3) - k2h*z(4) -d2h(z'1h-z'1h)-d1h(z'1v-h'h);
    end

 function n = MM1(t, y)
                dy = RM%zeros(4,1); %RM
                %n1 = (m2v +m2k(lh/lg)^2)

            end