Forum de mathématiques - Bibm@th.net

chahid

Membre

Inscription : 02-02-2021

Messages : 1

[SCILAB] Correction d'un script

qui peut m'aider à trouver les erreurs dans le script ci-dessous et mrc:

clf
//unsteady state column

//List of equations to be solved

//variable speci

//tempfin=var(4);

//Data input
plates=11;
LtoD=3;
Pend=(LtoD)/(LtoD+1);

//list of xd to solve the integral
ydest=[0.75:0.01:0.85];
yd=0.83;
xf=0.11;
Feed=350;

//polinomial
coea2=-16.78846;
coeb2=5.35043771;
coec2=0.07820301;

coea1=0.2774380;
coeb1=0.21006435;
coec1=0.4766639;
//antoine coef
Awa=18.3036;
Bwa=3816.44;
Cwa=227.02;

Aet=18.9119;
Bet=3803.98;
Cet=231.47;

//integration
for i=1:1:11
    yd=ydest(i);
    yeq(i,1)=yd;
   
    for j=1:1:plates
        if yeq(i,j)>0.5044
            then
            xeq(i,j)=(-coeb1+(coeb1^2-4*(coea1*(coec1-yeq(i,j))))^0.5)/(2*coea1);
           
        else
        xeq(i,j)=(-coeb2+(coeb2^2-4*(coea2*(coec2-yeq(i,j))))^0.5)/(2*coea2);
    end
   
        if j<11
            then
            yeq(i,j+1)=Pend*xeq(i,j)+(1-Pend)*yeq(i,1);
    end
end
xw(i)=xeq(i,11);
inv(i)= 1/(ydest(i)-xw(i))
end

// integral computation
for i=1:1:10
difer(i) =inv(i)*(xw(i+1)-xw(i));
end

coefint=polyfit(xw,inv,2);
Integral=(coefint(1)/3)*(xf.^3-xres.^3)+(coefint(2)/2)*(xf.^2- xres.^2)+(coefint(3))*(xf-xres)

function resid=rectdis(var)

    ydest=[0.75:0.01:0.85];
    yd=0.83;
    xf=0.11;
    Feed=350;
    xres=var'(1);
    Residue=var'(2);
    Distillate=var'(3);
   
   
    resid(1)=Residue*exp(Integral)-Feed;
    resid(2)=Feed-Distillate-Residue;
    resid(3)=Feed+xf-Distillate*yd-Residue*xres;
endfunction

a=fsolve([0.05,300,50],rectdis);

yoshi

Modo Ferox

Inscription : 20-11-2005

Messages : 16 992

Re : [SCILAB] Correction d'un script

Bonjour,

Les Règles de Bibmath ne sont pas là, juste pour faire joli...
Les membres sont censés les avoir lues et... les respecter
Extraits :

Présentation
  L'objectif de BibM@th est de créer un lieu d'échange, d'entraide, d'information ouvert à tous. Les utilisateurs sont invités à faire de ce forum un moyen de communication convivial, ouvert. Tout message se doit donc de contenir les formules de "politesse" en usage dans les rapports sociaux : Bonjour, (Bonsoir, Salut), s'il vous plaît, merci...
En cas d'oubli (!), un modérateur (ou l'Administrateur) répondrait en vous incitant à reformuler votre question, fermerait la discussion, et passé un délai de quelques jours, la supprimerait.
(...)
Comment bien poster
(...)
  * Notre but étant de vous aider à résoudre vos difficultés, et non de faire les exercices à votre place, ne postez pas le sujet d'un exercice sans montrer que vous y avez travaillé : il n'y serait probablement pas répondu. A vous d'expliquer ce que vous avez déjà fait, là où vous bloquez, et pourquoi...
(...)

Merci de reposter en respectant notre Charte...

Sujet fermé.

     Yoshi
- Modérateur -

---

Arx Tarpeia Capitoli proxima...

droledeg

Invité

Re : [SCILAB] Correction d'un script

bonjour,
je viens d'installer scilab sur mon mac, j'essaye(première fois) de faire un petit programme(tout petit)
a savoir n=12. i et j réels 1<=i<j<=12
x=1/(i+j)
dis(x).  je suis sur sciottes,je clique il faut enregistrerait faut un nom et apparait:vous n'avez pas lesdits nécessaires????????
sos.  merci

droledeg

Invité

Re : [SCILAB] Correction d'un script

bonjour

  j(ai écrire je suis sur scilab

merci

yoshi

Modo Ferox

Inscription : 20-11-2005

Messages : 16 992

Re : [SCILAB] Correction d'un script

Bonsoir

je suis sur sciottes

sciottes ???
kessékça ?

Au Passage, tu aurais ouvrir TA discussion (la tienne) pas parasiter celle des autres, en cliquant ici : Nouvelle discussion

Sujet fermé

      Yoshi
- Modérateur -

---

Arx Tarpeia Capitoli proxima...