The output has the correct amplitude, since the integral's amplitude is frequency-dependent:
$$
\int \sin(2\pi f t){\,\rm d}t = -\frac{1}{2\pi f} \cos(2\pi f t) + c
$$
Note the integral's amplitude \$1/(2\pi f)\$: that's the scaling you observed.
To get a unity output, the frequency must be \$ f=1/2\pi \$:
$$
\int \sin \left(2\pi \frac{1}{2\pi} t \right){\,\rm d}t = -\cos(t) + c
$$
Also note that the integrator is free to add an arbitrary offset voltage to the output.
We can simulate this using CircuitLab simulator built into this site (!):

simulate this circuit – Schematic created using CircuitLab
We can observe that the resulting \$ c=1 \$:

How can I adjust the size of an integrated signal?
If you want to normalize the output amplitude to be equal to the input amplitude, must multiply the integrator's output by the \$ 2\pi f \$ factor:
$$ H(s) = \frac{2 \pi f}{s}. $$
We may as well remove the DC offset, which is equal to the amplitude of the sine wave:

simulate this circuit

I highly suggest playing with CircuitLab - for such simple exploration, it's entirely adequate and you don't need Simulink.
It is also very flexible, since most of its functions support variables, so e.g. you can set up the frequency and amplitude as variables, and use them in all kinds of sources, Laplace block polynomials, passive component parameters, active component model parameters, etc.
The variables \$f\$ and \$A\$ are interactively adjustable, and the transient response time domain simulation is already set up, so it's very easy to play with those circuits.