Programación matlab

Ejercicio Matlab
El área del recinto limitado por las curvas y=f(x), x=c, x=d e y=0, puede calcularse de forma aproximada por una suma de Riemann, es decir, si se divideel intervalo [c,d] en n-1 subintervalos por medio de n puntos x1=c, x2,……xn=d, se tendría:
Area =
Se trata de escribir un programa en MATLAB que calcule el área dela elipse de ecuación x^2/a^2 + y^2/b^2 =1 ; utilizando el método anterior. Además, dibujará en color negro y trazo continuo dicha elipse.
Normas de programación1)Introducción por teclado de los valores de los semiejes a y b, y el número de puntos de la partición n (el primer y último punto de la partición serán el extremo izquierdo yderecho del intervalo respectivamente).
2)Se considerará que todos los subintervalos tienen la misma longitud.
3)Aplicando el método explicado se debe calcular el área de lasemielipse superior (cortada por y=0).
4)Debe aparecer por pantalla el valor del área total de la elipse y se dibujará la elipse como se ha indicado.

Solución:%Lectura de datos por teclado
n = input(‘Introduce el número de puntos de la partición \n’);
a = input(‘Introduce el valor del semieje a \n’);
b = input(‘Introduce el valordel semieje b \n’);

%Intervalo en el eje x
x = linspace(-a,a,n);
subintervalo = x(2) – x(1);

%cálculo de los puntos que definen la semielipse superior y el valordel sumatorio
suma = 0;
for i=1:n-1
y(i) = b*sqrt(1 – (x(i)*x(i))/(a*a));
suma = suma + y(i)*subintervalo;
end

%cálculo del área total y escritura porpantalla
suma = 2*suma;
sprintf(‘El area aproximada de la elipse es: %f’, suma)

%dibujo de la elipse
y(n) = 0;
y_negativa = -y;
plot(x,y,’-k’,x,y_negativa,’-k’)
*