How to perform indefinite integration of this function?

Hey guys, I need to perform the following operations as shown in the image. I need to calculate the value of function H for different inputs(x).

I am giving the following command from Symbolic Math Toolbox

%xxxxxxxxxxxxxxxxxxxxxx

syms y t x;

f1=(1-exp(-y))/y;

f2=-t+3*int(f1,[0,t]);

f3=exp(f2);

H=int(f3,[0,x]);

%xxxxxxxxxxxxxxxxxxxxxxx

but the value of 2nd integral i.e. integral in the function H can't be calculated and my output is of the form of

H =

int(exp(3*eulergamma - t - 3*ei(-t) + 3*log(t)), t, 0, x)

If any of you guys know how to evaluate this or have a different idea about this, please share it with me.

Thank you.