kod2.m

ok &ouml nemli bir kaybnkao lasdhlas sdohfs naosfnsdn sdondsfgpdsmr ljrg fdbh

clear all
n1 = input ('n1= ');
n2 = input ('n2= ');
a = input ('a = ');
n = input ('n = ');


fidd=fopen('N.txt','w'); 
fprintf(fidd,'%f\n',n);
kapat=fclose(fidd);

delta=sqrt((n1^2-n2^2)/(2*(n1^2)));
w0=sqrt(a/(n1*sqrt(2*delta)));

N=449;
t1=-3*a;
t2=3*a;
aralik=(t2-t1)/449;

i=1;
for j=t1:aralik:t2
   X(i)=j;
   xx=(j*sqrt(2)/w0);
   yy=-(j^2)/(w0^2);
   hg=hermite(n,xx);
   Y(i)=(power(2,0.25))/(power(pi,0.25)*sqrt(power(2,n-1)*faktoryel(n-1)*w0))*hg(n)*exp(yy);
   i=i+1;
end
sonuc2;
plot(X,Y);

     i=1;
       for j=t1:aralik:t2
         X(i)=j;
         xx=(j*sqrt(2)/w0);
   		yy=-(j^2)/(w0^2);
		   hg=hermite(n,xx);
		   Y(i)=(power(2,0.25))/(power(pi,0.25)*sqrt(power(2,n-1)*faktoryel(n-1)*w0))*hg(n)*exp(yy);
         HxHx(i)=Y(i)*Y(i);
         i=i+1;
      end
      
      Int=tabint(X,HxHx,N);

      for j=1:1:n
         Kokler(j)=sqrt((2*pi)^2*n1^2-2*(2j+1)/(w0^2));
      end
      
      for j=1:1:n
	   		fprintf('Kokler(%2i)=  %14.8f\n',j,Kokler(j));
      end
         
  		fprintf('w0=  %14.8f\n',w0);
  		fprintf('Int=  %14.10f\n',Int);