Fonctions

Fonctions de Bessel de 2ème espèce.

%-------------------------------------------------------------------------
% utilisation fonction bessely
% d'aprés 'Apprendre et maitriser Matlab'
% erreur corrigée
%-----------------------------------------------------------------------
% f de bessel d'ordre v
% sol particulière de y"+(1/t)y'+(1-[v/t]²)y=0
% sol générale
%        y = A J (t) + B Y (t)
%               v         v
%   Jv : Bessel 1ère espèce   besselj
%   Yv : Bessel 2ème espèce   bessely
% avec                        i      2i
%              v   inf     (-1) [t/2]
%    J  = [t/2]   Sigma------------------
%     v            i=0   i! gamma(v+i+1)
%
%             cos(v pi) J (t) - J  (t)
%                        v      (-v)
%    Y  = -------------------------------
%     v               sin(v pi )
%
%-------------------------------------------------------------
echo on;
%--------------------------------------------------------------
%  fonction Bessel 2ème espèce
% voir contenu fichier pour explications
%--------------------------------------------------------------
echo off;
clear all;
t=5:0.5:15; 		% domaine
nu=0:5;     		% ordre nu varie de 0 à 5 pour comparaison
b2=bessely(nu,t'); 	% erreur corrigée t'=transposée de t
mesh(nu,t,b2); 		% maillage3D
xlabel('ordre');
ylabel('t');
zlabel('Ynu(t)');
title('Bessel 2ème espèce');