Scilab

Bonjour;
je suis stagiaire en matériaux composites et j'ai réalisé dans un premier temps un programme sur scilab mais l'encadrant a voulue que je fasses un fichier d'entrés et les resultats doit etre enregistrer aussi sous un fichier excel automatiquement dès l'execution du programme qui est comme suite :
function [ntot,nl,recL,recC,p,nc,N,s,ef,etot,el,ec,alpha]=dimensionnement(r,Rv,Pr,Vf,Cfl,Cfc,Mv,Tb,Nb,Lm)
// INPUT :
// r=81/2 : r est le rayon des ouvertures (mm)
// Rv=242/2 : est le rayon de la virolealpha=asin(r./Rv) (mm)
// Pr=52*10^5 : pression à la rupture (pression interne) (Pa)
// Vf=0.6 : fraction volumique des fibres
// Cfl=2.95*10^9 : contrainte longitudinale
// Cfc=2.95*10^9 : contrainte circonfirenciel
// Tb=4450 : le grammage est la masse en g de 10000m de fil (g/10^4 m soit 10^-7g/mm)
// Nb=2 : nombre de bobines
// Mv=1.81 : Masse volumique (kg/dm^3 soit 10^-3g/mm^3)
// Lm=8 (mm)
// OUTPUT :
//alpha : l'angle entre l'axe et la mêche (rad)
//el : épaisseur la couche longitudinale (mm)
//ec : épaisseur de la couche circulaire (mm)
//etot : épaisseur du composite
//ef : épaisseur du fil sec
//ntot : nombre de couches du composite
//nl : nombre de couches de longitudinale
//nc : nombre de couches de circulaire
//s : section de la mêche (mm^2)
//Ncycles : nombre de cycles longitudinale
//p : pas du bobinage circulaire (mm)
alpha=asin((r+(Lm/2))/Rv)
el=(Pr*Rv)/(2*(cos(alpha))^2*Vf*Cfl)
ec=((Pr*Rv)/(Vf*Cfc))*(1-((tan(alpha))^2)/2)
etot=ec+el
ef=etot*Vf
ntot=ceil(etot/ef)
nl=ceil(el/ef)
nc=ceil(ec/ef)
s=(Tb*Nb)/(10^4*Mv)
N=ceil((el*%pi*Rv*Vf*cos(alpha))/s)
p=ceil((s*nc)/(ec*Vf))
recL=(1-(((2*%pi*Rv*nl)/(N*Lm))*cos(alpha)))*100
recC=((Lm-p)/p)*100
endfunction

// tableau des entrées est comme suite:
T_ent=["r","81/2";"Rv","242/2";"Pr","52*10^5";"Vf","0.6";"Cfl","2.95*10^9";"Cfc","2.95*10^9";"Mv","1.81";"Tb","4450";"Nb","2";"Lm","8"]
// pour modifier les valeurs d'entrer veuillez agir dans l'équation ci dessous tout en respectant l'ordre des parametres
[ntot,nl,recL,recC,p,nc,N,s,ef,etot,el,ec,alpha]=dimensionnement(81/2,242/2,52*10^5,0.6,2.95*10^9,2.95*10^9,1.81,4450,2,8)

// pour executer : dans le fichier text (ctrl+l) puis dans le consol cliquez sur (entrer)