**Good Morning,**

**This is the script that I have used to to an integration, as the command ''int'' didn't work because the coputation was really heavy.**

p_old := 0;

g := 0; #index

E_tot_net_Point := 0;

for g to 2 do

h := g+1;

loopstart := SP_time[g];

loopend := SP_time[h];

for tt from loopstart by 0.01 to loopend do

t := tt;

P_tot_net_Fkt := t -> C_motore[g]*omega_motore[g]

if t≠SP_time[g] and t ≠ SP_time[h] then

P_tot_net_Point := P_tot_net_Fkt(t);

P_tot_net_Point_next := P_tot_net_Fkt(t+0.001);

plt_P_tot_net_plt[t] := pointplot([t, P_tot_net_Point]);

E_tot_net_Point := E_tot_net_Point+(P_tot_net_Point*0.01);

plt_E_tot_net_plt[t] := pointplot([t, E_tot_net_Point]);

p_old := P_tot_net_Point;

end if;

end do:

end do:

plots:-display(entries(plt_E_tot_net_plt, 'nolist'), title = "Energy", color = "Red", symbol = point, symbolsize = 20, gridlines = true, size = [1000, 500], labels = ['[s]', '[Watt*s]'], font = [axes, 12])

**I have plotted the function as a series of points over the orizontal axis (t).**

**Now I would like to calculate the last value of the curve (which is itself a point) as a function of the orizontal axes t. How can I do it?**

**Something like E_tot_net_Point = f(t)**

** **

**Thank you for the help**