function yt1=exm05835_2(flag)
%exm05835_2.m
%flag取 1 時，只進行符號計算
%flag取 2 時，只進行SIMULINK仿真
%flag取 3 時，把SIMULINK仿真結果與解析結果比較

switch flag
case 1			%（1）符號卷積法得解析結果
   syms tao;t=sym('t','positive');		%把t定義為"取正"符號變量
   US1=laplace(exp(-t));				%u(t)的L氏變換
   HS1=laplace(t*exp(-t/2))				%h(t)的L氏變換
   yt1=simple(ilaplace(US1*HS1))		%L氏反變換得卷積的理論解
case 2		   %（2）SIMULINK卷積法
   exm05835_2_2
   yt1=[];
case 3			%（3）計算結果比較
   syms tao;t=sym('t','positive');		%把t定義為"取正"符號變量
   US1=laplace(exp(-t));				%u(t)的L氏變換
   HS1=laplace(t*exp(-t/2))				%h(t)的L氏變換
   yt1=simple(ilaplace(US1*HS1))		%L氏反變換得卷積的理論解
   [t,yt2]=sim('exm05835_2_2');
   yyt1=eval(vectorize(char(yt1)));	%理論解的數值序列
   [dy,kd]=max(abs(yyt1-yt2(:,2)));
   dy12=dy/abs(yyt1(kd))
end