[t,y,yprime,delta,ires,rpar,ipar]=sdres1(t,y,yprime,delta,ires,rpar,ipar);

function [t,y,yprime,delta,ires,rpar,ipar]=sdres1(t,y,yprime,delta,ires,rpar,ipar);

ipar_shape=size(ipar);ipar=reshape(ipar,1,[]);
y_shape=size(y);y=reshape(y,1,[]);
yprime_shape=size(yprime);yprime=reshape(yprime,1,[]);
delta_shape=size(delta);delta=reshape(delta,1,[]);
rpar_shape=size(rpar);rpar=reshape(rpar,1,[]);
delta(1) = yprime(1) + 10.0e0.*y(1);
delta(2) = y(2) + y(1) - 1.0e0;
ipar_shape=zeros(ipar_shape);ipar_shape(:)=ipar(1:numel(ipar_shape));ipar=ipar_shape;
y_shape=zeros(y_shape);y_shape(:)=y(1:numel(y_shape));y=y_shape;
yprime_shape=zeros(yprime_shape);yprime_shape(:)=yprime(1:numel(yprime_shape));yprime=yprime_shape;
delta_shape=zeros(delta_shape);delta_shape(:)=delta(1:numel(delta_shape));delta=delta_shape;
rpar_shape=zeros(rpar_shape);rpar_shape(:)=rpar(1:numel(rpar_shape));rpar=rpar_shape;
end %subroutine sdres1