Hola, necesito hacer un programa que calcule el area que encierra la grafica de la funcion f(x)=sen(ax) con "x" de 0 a 2*pi, y "a" de 0 a 10.
El problema es que había diferenciado entre valores de "a" pares e impares y aún así había algunos valores para los que la integral del area me salía cero, por tanto he distinguido otros casos pero es un poco "chapucero".
También pensé hacerlo con el valor absoluto, pero siempre sale el mismo resultado, indepentientemente del valor de la "a", y haciendolo a mano me salen diferentes valores.
Os pongo aquí mi código, por favor a ver si se os ocurre otro mejor:
x=linspace(0,2*pi);
y=sin(a*x);
g=plot(x,y);
syms x
a
if (rem(a,2)==0 & a~=4 & a~=8)
area=4*int(sin(a*x),x,0,pi/2)
elseif(a==4)
area=8*int(sin(a*x),x,0,pi/4)
elseif(a==8)
area=16*int(sin(a*x),x,0,pi/8)
else
area=2*int(sin(a*x),x,0,pi)
end
area=double(area);
title(['El area es ',num2str(area)])
Esto es parte de otro programa, por lo que los valores de la a se especifican en el apartado correspondiente.
Muchísimas gracias, agradeceré enormemente cualquier sugerencia.