clc
close all
x=[ ];
xx=x;tt=t;
minx=min(x);mint=min(t);maxx=max(x);maxt=max(t);
x=x';t=t';x=premnmx(x);t=premnmx(t);
x2=x;
[l,n]=size(x);
i=1;
error=10;
thick=2*1.9*0.0001;
m=2;
ss=0.005;s=ss;
q=0.5;
for p=1:50
    i=i+1;
    B=(1+q*s)*diag(ones(n,1))-q*s/2*diag(ones(n-1,1),1)-q*s/2*diag(ones(n-1,1),-1);
    A=(1-q*s)*diag(ones(n,1))+q*s/2*diag(ones(n-1,1),1)+q*s/2*diag(ones(n-1,1),-1);
    for j=1:m
        y=inv(B)*A*x2';
        x2=y';
    end
    y=y';
    error(i)=sum((y-t).^2).^0.5/n;
    ss(i)=ss(i-1)-0.1*ss(i-1)*error(i);
    s=ss(i);
      
end
figure;plot(3:n-3,t(3:n-3),3:n-3,y(3:n-3));  
dx=thick/m;
dt=1;
alfa=s*dx.^2/dt;
ny=inv(B)*A*x';
nt=inv(B)*A*postmnmx(x,minx,mint)';
Y=norm(postmnmx(y,mint,maxt));
alfa=Y/norm(ny)*norm(postmnmx(x,minx,mint))/norm(x)*alfa;
alfa
Errorpercent=error(i)*100