function [code] = impl_ltisys( sys ,name)
%IMPL_LTISYS generate implementation by ltisys library
%   

sys=ss(sys);

if sys.ts ~=0
    error('system should be continuous-time');
end

if nargin<2
    name='sysimpl';
end

nx=size(sys.A,1);
nu=size(sys.B,2);
ny=size(sys.C,1);

code = strcat(...
    sprintf('ltisys<%d,%d,%d> %s(',nx,nu,ny,name),...
    '(const double []) {', sprintf('%e,',sys.A'),'}, ',...
    '(const double []) {', sprintf('%e,',sys.B'),'}, ',...
    '(const double []) {', sprintf('%e,',sys.C'),'}, ',...
    '(const double []) {', sprintf('%e,',sys.D'),'}',...
    ');');
end

