Fonctions

Fonction de Bessel de 1ère espèce modifiée.

%-------------------------------------------------------------------
% utilisation fonction besselj
% 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
% on remplace t par l'imaginaire pur jt
% sol générale
%        y = A I (t) + B K (t)
%               v         v
%   Iv : Bessel 1ère espèce mod  besseli
%   Kv : Bessel 2ème espèce mod  besselk
% avec                           2i
%              v   inf      [t/2]
%    I  = [t/2]   Sigma------------------
%     v            i=0   i! gamma(v+i+1)
%
%            I  (t) - I (t)
%            (-v)      v
%    K  = --------------------
%     v      (2/pi)sin(v pi )
%
%---------------------------------------------------------------
echo on;
%---------------------------------------------------------------
%  fonction Bessel 1ère espèce modifiée
% voir contenu fichier pour explications
%-----------------------------------------------------------------
echo off;
clear all;
t=15:0.5:20; 		% domaine
nu=0:10;     		% ordre nu varie de 0 à 10 pour comparaison
b1m=besseli(nu,t'); 	% erreur corrigée t'=transposée de t
mesh(nu,t,b1m); 	% maillage3D
xlabel('ordre');

ylabel('t');
zlabel('Inu(t)');
title('Bessel 1ère espèce modifiée');