% test of spherical bessel function generator x=[0.1:0.1:10]; jf=zeros(length(x),6); for i=1:length(x) j=sbfj(x(i),10); jf(i,:)=j(1:size(jf,2)); end col=['r' 'g' 'b' 'c' 'm' 'y']; figure(1) hold off a=plot(x,jf(:,1),col(1)) set(a,'LineWidth',2); hold on for i=2:6 a=plot(x,jf(:,i),col(i)) set(a,'LineWidth',2); end axis([0 10 -2 2 ]); grid on title('j_m(x) - Backward Recurrence'); xlabel('x') jf=zeros(length(x),6); for i=1:length(x) j=sbfj_3(x(i),10); jf(i,:)=j(1:size(jf,2)); end col=['r' 'g' 'b' 'c' 'm' 'y']; figure(2) hold off a=plot(x,jf(:,1),col(1)) set(a,'LineWidth',2); hold on for i=2:6 a=plot(x,jf(:,i),col(i)) set(a,'LineWidth',2); end axis([0 10 -2 2 ]); grid on title('j_m(x) - Backward Recurrence - 3 digits'); xlabel('x') % forward recurrence jf=zeros(length(x),6); for i=1:length(x) j=sbfj_f(x(i),10); jf(i,:)=j(1:size(jf,2)); end col=['r' 'g' 'b' 'c' 'm' 'y']; figure(3) hold off a=plot(x,jf(:,1),col(1)) set(a,'LineWidth',2); hold on for i=2:6 a=plot(x,jf(:,i),col(i)) set(a,'LineWidth',2); end axis([0 10 -2 2 ]); grid on title('j_m(x) - Forward Recurrence'); xlabel('x') % forward recurrence 3 sign digits jf=zeros(length(x),6); for i=1:length(x) j=sbfj_f_3(x(i),10); jf(i,:)=j(1:size(jf,2)); end col=['r' 'g' 'b' 'c' 'm' 'y']; figure(4) hold off a=plot(x,jf(:,1),col(1)) set(a,'LineWidth',2); hold on for i=2:6 a=plot(x,jf(:,i),col(i)) set(a,'LineWidth',2); end axis([0 10 -2 2 ]); grid on title('j_m(x) - Forward Recurrence - 3 digits'); xlabel('x') % forward recurrence y jf=zeros(length(x),6); for i=1:length(x) j=sbfy(x(i),10); jf(i,:)=j(1:size(jf,2)); end col=['r' 'g' 'b' 'c' 'm' 'y']; figure(5) hold off a=plot(x,jf(:,1),col(1)) set(a,'LineWidth',2); hold on for i=2:6 a=plot(x,jf(:,i),col(i)) set(a,'LineWidth',2); end axis([0 10 -2 2 ]); grid on title('y_m(x) - Forward Recurrence'); xlabel('x')