حلش کردم نتایج p اونقد کوچیکه که تو نمودار اصلا مشخص نیست
کد:
clc
x=0.5,k=0.5,p=0;
for i=1:50
p(i+1)=p(i)+k*sin(2*pi*x(i))/(2*pi)
x(i+1)=x(i)+p(i+1)
end
plot(p,x)