[x,y,yp]=duivp(x,y,yp);

function [x,y,yp]=duivp(x,y,yp);
%***BEGIN PROLOGUE  DUIVP
%***PURPOSE  Dummy routine for DBVSUP quick check.
%***LIBRARY   SLATEC
%***TYPE      doubleprecision (UIVP-S, DUIVP-D)
%***KEYWORDS  QUICK CHECK
%***AUTHOR  Watts, H. A., (SNLA)
%***DESCRIPTION
%
%   This routine is never called;  it is here to prevent loaders from
%   complaining about undefined externals while testing DBVSUP.
%
%***ROUTINES CALLED  (NONE)
%***REVISION HISTORY  (YYMMDD)
%   750601  DATE WRITTEN
%   890618  REVISION DATE from Version 3.2
%   891214  Prologue converted to Version 4.0 format.  (BAB)
%   920401  Variables declaration and TYPE sections added.  (WRB)
%***end PROLOGUE  DUIVP
%     .. Scalar Arguments ..
%     .. Array Arguments ..

y_shape=size(y);y=reshape(y,1,[]);
yp_shape=size(yp);yp=reshape(yp,1,[]);
%***FIRST EXECUTABLE STATEMENT  DUIVP
error(['stop encountered in original fortran code  ',char(10),';']);
y_shape=zeros(y_shape);y_shape(:)=y(1:numel(y_shape));y=y_shape;
yp_shape=zeros(yp_shape);yp_shape(:)=yp(1:numel(yp_shape));yp=yp_shape;
end
%DECK DULSIA