Untitled

n=7;
N=101;
a=3/n;
x=[];
x1=[];
y=[];
z=[];
z1=[];
dx=[];
ddx=[];
 
 
 
 
for k=2:N;
    x(k)=a*sin(2*pi*n*(k-1)/100);
    z(k)=1.3*a*sin(4*pi*n*(k-1)/100);
 
end
 
for k=2:N-1;
    dx(k)=(x(k+1)-x(k-1))/2;
    ddx(k)=x(k+1)-2*x(k)+x(k-1);
    y(k)=(x(k)*dx(k))/sqrt(0.148*(dx(k))^2-0.148*x(k)*ddx(k));
end
 
for i.html">i=1:100
    x1(i.html">i)=x(i.html">i+1);
    z1(i.html">i)=z(i.html">i+1);
end
 
k=2:N;
 
subplot(221);
plot(k/1.3,x1/1.3);
x1
 
subplot(222);
plot(k/1.3,y/1.3);
y
 
subplot(223);
plot(k/1.3,z1/1.3);
z1