clc
g=inline('x^2')
a=3.0
b=4.9
m=350
val=trapezoidal_rule(g,a,b,m)
syms x
valexact=eval(int(char(g),x,a,b))