function [sys,x0,str,ts] = spacemodel(t,x,u,flag)switch flag,case 0,    [sys,x0,str,ts]=mdlInitializeSizes;case 3,    sys=mdlOutputs(t,x,u);case {2,4,9}    sys=[];otherwise    error(['Unhandled flag = ',num2str(flag)]);endfunction [sys,x0,str,ts]=mdlInitializeSizessizes = simsizes;sizes.NumContStates  = 0;sizes.NumDiscStates  = 0;sizes.NumOutputs     = 4;sizes.NumInputs      = 0;sizes.DirFeedthrough = 1;sizes.NumSampleTimes = 1;sys = simsizes(sizes);x0  = [];str = [];ts  = [0 0];function sys=mdlOutputs(t,x,u)x1d=sin(3*t);dx1d=3*cos(3*t);x2d=cos(3*t);dx2d=-3*sin(3*t);sys(1)=x1d;sys(2)=x2d;sys(3)=dx1d;sys(4)=dx2d;