Problems to compute a phase shift using sinusoidal functions

I am having problems trying to compute a phase shift using sinusoidal functions. For example I am trying to achieve this:

I can compute this giving it a vertical transition, but no matter where I add its incorrect when plotted along side the above equation without the phase shift.

This is without the attempted phase shift:

y = 3*sin(x*(2*pi/4))+2;

x = [0:0.1:100];

Answer 1

[Wayne King]

   x = 0:0.1:10;
   y1 = 3*sin(4*x); % phase of zero
   y2 = 3*sin(4*(x-pi/4)); % phase shift of \pi radians
   plot(x,y1); 
   hold on
   plot(x,y2,'r')