Dear Ser/Madame,
To create more points in your function, you could use an interpolation function like pchip or spline to create more points. A example how it works is:
clc; clear all; close all;
a = 0.5;
b = 0.7;
%%PART 1: not enough points
h1 = subplot(2,1,1);
x = linspace(0,1,20);
y = 1/2*sin(x*pi).*exp(x*pi);
xpoints=x(find(x>=a & x<=b));
ypoints=y(find(x>=a & x<=b));
plot(h1,x,y);axes(h1); hold on;
area(xpoints,ypoints); hold on;
line([a a],get(h1,'YLim'),'Color',[1 0 0]);
line([b b],get(h1,'YLim'),'Color',[1 0 1]);
title(h1,'title 1');
%%PART 2: more points
h2 = subplot(2,1,2);
xnew = linspace(min(x),max(x),length(x)*50);
ynew = pchip(x,y,xnew);
xpoints_new=xnew(find(xnew>=a & xnew<=b));
ypoints_new=ynew(find(xnew>=a & xnew<=b));
plot(h2,xnew,ynew);axes(h2); hold on;
area(xpoints_new,ypoints_new); hold on;
line([a a],get(h1,'YLim'),'Color',[1 0 0]);
line([b b],get(h1,'YLim'),'Color',[1 0 1]);
title(h2,'title 2');
If you like you can increase the number 50 to any number you like.
Good luck! Christiaan
