Matlab - Calculo de area de dos funciones

clear all,

t  = 0:pi/180:3*pi;

x = sin(t) ;

y = cos(t);

figure

plot(t,x)

hold on

plot(t,y)

trapz(x,y)

syms x

f1=input('función inferior:  ');

f2=input('función superior: ');

a=input('límite inferior: ');

b=input('límite superior: ');

area=int(f2-f1,x, a, b)

ff1=str2func(['@(x)' char(f1)]);

ff2=str2func(['@(x)' char(f2)]);

x=[];

y=[];

for i=a:0.1:b

x(end+1)=i;

y(end+1)=ff2 (i);

end

for i=b:-0.1:a

x(end+1)=i;

y(end+1)=ff1(i);

end

fill(x,y,'r')