Matlab - ayuda area bajo la curva con formula del trapezio

hice un programa para calcular el area bajola curva mediante la formula del trapezio dividiendola en subintervalos para que se acercara al area exacta, pero la funcion que pongo queda determinada por mi, lo que kiero es ke el usuario pueda determinar la funcion ke sea. la ke el usuario quiera. si alguien me puede ayudar le agradeceria mucho :)

pongo el codigo para que seden una idea mas o menos:

% Programa para integrar las funcinoes de tarea
% para acercar la forma del trapezio mediante subintervalos
% al area exacta de la funcion, definiendo elintervalo [a b]
% y la cantidad de subintervalos que queramos (n)
clc
%aki kiero ke el usuario inserte la funcion ke el kiera, y ke el area se kalkule para la funcion ke el introdujo

a=input('Defina el limite inferior [a ---> ')
b=input('Defina el limite superior b]---> ')
n=input('Defina la cantidad de subintervalos --->')


x=a:(b-a)/n:b;
A1=0;
A2=0;
A3=0;
A4=0;
A5=0;
for k=1:n
An1=((x(k)^4+x(k+1)^4)/2)*(x(k+1)-x(k));%area del intervalo
A1=A1+An1; %area acumulada de la funcion x^4

An2=((exp(x(k))+exp(x(k+1)))/2)*(x(k+1)-x(k));
A2=A2+An2; %Area de f(x)=e^x

An3=((sin(x(k))+ sin(x(k+1)))/2)*(x(k+1)-x(k));
A3=A3+An3; %Area de f(x)=sen (x)

An4=((sqrt(1+x(k)^2)+sqrt(1+x(k+1)^2))/2)*(x(k+1)-x(k));
A4=A4+An4; %Area de f(x)=sqrt(1+x^2)

An5=((1/(x(k)+1)+1/(x(k+1)+1))/2)*(x(k+1)-x(k));
A5=A5+An5; %Area de f(x)=1/(x+1)

end
disp('El Area de x^4')
A1
disp('El Area de e^x')
A2
disp('El Area de sen (x)')
A3
disp('El Area de sqrt(1+x^2)')
A4
disp('El Area de 1/(x+1)')
A5