Projet enceintes 4 voies

[tonipe]

J'ai créé la fonction Bessel d'ordre 1.
Je ne sais pas encore si mon nombre de lignes est suffisant. Pour des valeur entre 0 et 3, oui. Au delà ce n'est pas sur.
Il y a surement moyen de programmer ça plus subtilement.
La dernière ligne, c'est un +.

Code :

<?php
  function bessel($x)
  {
    $y = $x/2;
    $y = $y-pow($x,3)/pow(2,2)/4;
    $y = $y+pow($x,5)/pow(2,2)/pow(4,2)/6;
    $y = $y-pow($x,7)/pow(2,2)/pow(4,2)/pow(6,2)/8;
    $y = $y+pow($x,9)/pow(2,2)/pow(4,2)/pow(6,2)/pow(8,2)/10;
    $y = $y-pow($x,11)/pow(2,2)/pow(4,2)/pow(6,2)/pow(8,2)/pow(10,2)/12;
    $y = $y+pow($x,13)/pow(2,2)/pow(4,2)/pow(6,2)/pow(8,2)/pow(10,2)/pow(12,2)/14;
    $y = $y-pow($x,15)/pow(2,2)/pow(4,2)/pow(6,2)/pow(8,2)/pow(10,2)/pow(12,2)/pow(14,2)/16;
    $y = $y+pow($x,17)/pow(2,2)/pow(4,2)/pow(6,2)/pow(8,2)/pow(10,2)/pow(12,2)/pow(14,2)/pow(16,2)/18;
    $y = $y-pow($x,19)/pow(2,2)/pow(4,2)/pow(6,2)/pow(8,2)/pow(10,2)/pow(12,2)/pow(14,2)/pow(16,2)/pow(18,2)/20;
    $y = $y-pow($x,21)/pow(2,2)/pow(4,2)/pow(6,2)/pow(8,2)/pow(10,2)/pow(12,2)/pow(14,2)/pow(16,2)/pow(18,2)/pow(20,2)/22;
    return $y;
  }
?>

Forr :
Si Jean-Marc à mis une fonction BESSEL d'ordre 1, je calcule avec.
Je n'ai pas le niveau pour redémontrer la formule, je ferai le maximum pour la calculer le mieux possible, le plus utilement possible.
Juste pour le fun, vous savez calculer Pi à la main ? Et si non, ça ne vous pose pas de problème pour l'utiliser ?

Cordialement, Dominique