problème parenthèse de mon code scilab pour résoudre un système d’équations non linéaires

**tchiekre michel** · 21/11/2017, 12h20

bonjour, je suis nouveau dans le forum donc je sais pas trop comment l utiliser. j essaie de résoudre un système non linéaire avec scilab . mais je n y arrive pas car il me dit qu il ya une parenthèse qui manque a droite. Je suis débutant en programmation. urgent aidez moi

je met tout mon code. tout marche jusqu a la partie ou je dois resoudre mon systeme.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
/////////////////&/&/&/&/&/&/&/&/&/&/&/&/&/&/&/&/&//&/&/&/ LES CONSTANTES ET PARAMETRES DE MON CODE//&/&/&&//&/&&//&/&/&/&/&/&/&&//&/&&//&/&/&/&/&/&/&/&/&/&/&/&/
////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
 
 
///////////////////////////////////////////////////////////les Longueurs////////////////////////////////////////////////////////////////////////////////////////
 
//longueur de la flamme //Diametre du site combustible // Hauteur couche combustible 
L_fl= 4 ; D= 10^-3; H_c= 0.21 ; 
 
 
//////////////////////////////////////////////////////////Les temperatures////////////////////////////////////////////////////////////////////////////////////
 
//Temperature flamme Temperature ambiante Evolution Temperature experimentale Temperature inflammation temperature braise
T_fl=1083 ; T_inf= 305; dt_T_exp=((T_inf -500)/(103-48))*1.2; Tp_inf=500; T_b=561; 
 
///////////////////////////////////////////////Les vitesses et directions////////////////////////////////////////////////////////////////////////////////////
//Vitesse du vent
Uw=4.6;
//////////////////////////////////////////////DISTANCE Y EN FONCTION DU TEMPS/////////////////////////////////////////////////////
function y=Y(t),y= -0.9484*t+ 103.98, endfunction 
 
//////////////////////////////////////////////////Les coefficients/////////////////////////////////////////////////////////////////////////////////////////
//Emissivite flamme VISCOSITE CINEMATIQUE VISCOg(m²/s) A 500°C Fraction voluique couche combustible Emissivite braise 
 
esp_fl= 1-exp(-0.6*L_fl) ; nu_g=7.77*10^(-5); phi= 0.0026; esp_b= 1 ; 
 
// NBRE DE PRANDLT PR EMISSIVITé DE LA COUCHE COMBUSTIBLE 
 
Pr= 0.71 ; esp_c=0.6 ;
sigma= 5.67*10^(-8); R_eD=((1-phi)*(Uw*D)/(nu_g)) ; 
 
/////////////////////////////////////////////REYNOLDS EN FONCTION DU TEMPS/////////////////////////////////////////////////
 
function y=R_ey(t) , y=(Uw/nu_g)*((-0.9484)*t+ 103.98), endfunction 
 
//////////////////////////////////////////////////Les chaleurs specifiques///////////////////////////////////////
C_pc=1.25*10^3 ; 
 
///////////////////////////////////////////////////Les caracteristiques de l'air/////////////////////////////////////
k_fl= 7.07*10^(-3) ; k_b= 0.0454; 
 
////////////////////////////////////////////////Les proprietes physiques du combustible//////////////////////////////////////////////////////
phi= 0.0026 ; rho_c= 512 ; a_fb= 0.6 ; s= 25.40;
//////////////////////////////////////FACTEUR DE FORME EXPRESSION EN FONCTION DU TEMPS////////////////////////////////////////////////////////////////////////
function y=F(t),y= (10^(-12))*t^6 -4*(10^(-10))*t^5 +4*(10^(-8))*t^4 -2*(10^(-6))*t^3 +5*(10^(-5))*t^2 -0.0005*t +0.0012, endfunction 
 
///////////////////////////////////////////// TEMPERATURE T EN FONCTION DU TEMPS/////////////////////////////////////////////////////////////////////
function y=T(t),y= 6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13, endfunction 
 
 
//////////////////////////////////////////////////////////// LES INTEGRALES A CALCULER//////////////////////////////////////////////////////
 
// /////// I_11///////////////////////////
function y=i11(t),y=((10^(-12))*t^6 -4*(10^(-10))*t^5 +4*(10^(-8))*t^4 -2*(10^(-6))*t^3 +5*(10^(-5))*t^2 -0.0005*t +0.0012)^2, endfunction
I11=intg(53,100,i11)
disp(I11,'I11=')
 
////////// I_12/////////////////////////
function y=i12(t),y= ((10^(-12))*t^6 -4*(10^(-10))*t^5 +4*(10^(-8))*t^4 -2*(10^(-6))*t^3 +5*(10^(-5))*t^2 -0.0005*t +0.0012)*exp(-0.25*s*(-0.9484*t+ 103.98)),endfunction
I12=intg(53,100,i12)
disp(I12,'I12=')
 
//////////////////I_13///////////////////
function y=i13(t),y= ((10^(-12))*t^6 -4*(10^(-10))*t^5 +4*(10^(-8))*t^4 -2*(10^(-6))*t^3 +5*(10^(-5))*t^2 -0.0005*t +0.0012)*((6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13)^4-(T_inf)^4), endfunction
I13=intg (53,100,i13)
disp(I13,'I13=')
 
/////////////////I_14////////////////////
function y=i14(t),y=(((10^(-12))*t^6 -4*(10^(-10))*t^5 +4*(10^(-8))*t^4 -2*(10^(-6))*t^3 +5*(10^(-5))*t^2 -0.0005*t +0.0012)*(T_fl-(6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13))*exp(-0.35*(-0.9484*t+ 103.98)/ L_fl)*((Uw/nu_g)*((-0.9484)*t+ 103.98))^(0.5))/(-0.9484*t+ 103.98), endfunction
I14=intg(53,100,i14)
disp(I14,'I14=')
 
////////////////I_15////////////////////
function y=i15(t),y= ((10^(-12))*t^6 -4*(10^(-10))*t^5 +4*(10^(-8))*t^4 -2*(10^(-6))*t^3 +5*(10^(-5))*t^2 -0.0005*t +0.0012)*(T_b-(6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13))*exp(-0.25*s*(-0.9484*t+ 103.98)), endfunction
I15=intg(53,100,i15)
disp(I15,'I15=')
 
////////////////I_16////////////////////
function y=i16(t),y=(10^(-12))*t^6 -4*(10^(-10))*t^5 +4*(10^(-8))*t^4 -2*(10^(-6))*t^3 +5*(10^(-5))*t^2 -0.0005*t +0.0012, endfunction
I16=intg(53,100,i16)
disp(I16,'I16=')
 
////////////////I_21////////////////////
function y=i21(t),y=((10^(-12))*t^6 -4*(10^(-10))*t^5 +4*(10^(-8))*t^4 -2*(10^(-6))*t^3 +5*(10^(-5))*t^2 -0.0005*t +0.0012)^2, endfunction
I21=intg(53,100,i21)
disp(I21,'I21=')
 
////////////////I_22////////////////////
function y=i22(t),y=((10^(-12))*t^6 -4*(10^(-10))*t^5 +4*(10^(-8))*t^4 -2*(10^(-6))*t^3 +5*(10^(-5))*t^2 -0.0005*t +0.0012)*exp(-0.25*s*(-0.9484*t+ 103.98)), endfunction
I22=intg(53,100,i22)
disp(I22,'I22=')
 
////////////////I_23////////////////////
function y=i23(t),y=((10^(-12))*t^6 -4*(10^(-10))*t^5 +4*(10^(-8))*t^4 -2*(10^(-6))*t^3 +5*(10^(-5))*t^2 -0.0005*t +0.0012)*((6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13)^4-(T_inf)^4), endfunction
I23=intg(53,100,i23)
disp(I23,'I23=')
 
////////////////I_24////////////////////
function y=i24(t),y=(((10^(-12))*t^6 -4*(10^(-10))*t^5 +4*(10^(-8))*t^4 -2*(10^(-6))*t^3 +5*(10^(-5))*t^2 -0.0005*t +0.0012)*(T_fl-(6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13))*exp(-0.3*(-0.9484*t+ 103.98)/ L_fl)*((Uw/nu_g)*((-0.9484)*t+ 103.98))^(0.5))/(-0.9484*t+ 103.98), endfunction
I14=intg(53,100,i14)
disp(I14,'I14=')
 
/////////////////I25///////////////////
function y=i25(t),y= ((10^(-12))*t^6 -4*(10^(-10))*t^5 +4*(10^(-8))*t^4 -2*(10^(-6))*t^3 +5*(10^(-5))*t^2 -0.0005*t +0.0012)*(T_b-(6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13))*exp(-0.25*s*(-0.9484*t+ 103.98)), endfunction
I25=intg(53,100,i25)
disp(I25,'I25=')
 
/////////////////I26///////////////////
function y=i26(t),y= (((10^(-12))*t^6 -4*(10^(-10))*t^5 +4*(10^(-8))*t^4 -2*(10^(-6))*t^3 +5*(10^(-5))*t^2 -0.0005*t +0.0012)*((Uw/nu_g)*((-0.9484)*t+ 103.98))^(0.5)*exp(-0.3*(-0.9484*t+ 103.98)/(L_fl)))/(-0.9484*t+ 103.98), endfunction 
I26=intg(53,100,i26)
disp(I26,'I26=')
 
///////////////I27/////////////////////
function y=i27(t),y=((exp(-0.25*s*(-0.9484*t+ 103.98)-0.3*(-0.9484*t+ 103.98)/ L_fl ))*((Uw/nu_g)*((-0.9484)*t+ 103.98))^(0.5))/(-0.9484*t+ 103.98), endfunction
 
I27=intg(53,100,i27)
 
disp(I27,'I27=')
 
///////////////I28/////////////////////
function y=i28(t), y= (((Uw/nu_g)*((-0.9484)*t+ 103.98))^(0.5)*exp(-0.3*(-0.9484*t+ 103.98)/ L_fl)*((6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13)^4-( T_inf)^4))/(-0.9484*t+ 103.98), endfunction 
 
I28=intg(53,100,i28)
 
disp(I28,'I28=')
 
///////////////I29//////////////////////
function y=i29(t),y=((((Uw/nu_g)*((-0.9484)*t+ 103.98)*( T_fl-(6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13))*exp(-0.6*((-0.9484)*t+ 103.98)/L_fl)))/((-0.9484)*t+ 103.98)^2), endfunction
 
I29=intg(53,100,i29)
 
disp(I29,'I29=')
 
///////////////I210//////////////////////
function y=i210(t),y=(T_b-(6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13))*exp((-0.25*s*((-0.9484)*t+ 103.98)-(0.3*((-0.9484)*t+ 103.98))/L_fl)*((Uw/nu_g)*((-0.9484)*t+ 103.98))^(0.5))/((-0.9484)*t+ 103.98), endfunction
 
I210=intg(53,100,i210)
 
disp(I210,'I210=')
 
///////////////I211//////////////////////
function y=i211(t),y= (((Uw/nu_g)*((-0.9484)*t+ 103.98))^(0.5)*(exp((-0.3*( -0.9484*t+ 103.98))/L_fl)))/( -0.9484*t+ 103.98), endfunction
 
I211=intg(53,100,i211)
 
disp(I211,'I211=')
 
////////////////I212///////////////////////
function y=i212(t),y=(10^(-12))*t^6 -4*(10^(-10))*t^5 +4*(10^(-8))*t^4 -2*(10^(-6))*t^3 +5*(10^(-5))*t^2 -0.0005*t +0.0012, endfunction
I212=intg(53,100,i212)
disp(I212,'I212=')
 
////////////////I311///////////////////////
function y=i311(t),y=((10^(-12))*t^6 -4*(10^(-10))*t^5 +4*(10^(-8))*t^4 -2*(10^(-6))*t^3 +5*(10^(-5))*t^2 -0.0005*t +0.0012)^2, endfunction
I311=intg(53,100,i311)
disp(I311,'I311=')
 
////////////////I312///////////////////////
function y=i312(t),y=((10^(-12))*t^6 -4*(10^(-10))*t^5 +4*(10^(-8))*t^4 -2*(10^(-6))*t^3 +5*(10^(-5))*t^2 -0.0005*t +0.0012)*exp(-0.25*s*(-0.9484*t+ 103.98)), endfunction
I312=intg(53,100,i312)
disp(I312,'I312=')
 
/////////////I313/////////////////////////
function y=i313(t),y=((10^(-12))*t^6 -4*(10^(-10))*t^5 +4*(10^(-8))*t^4 -2*(10^(-6))*t^3 +5*(10^(-5))*t^2 -0.0005*t +0.0012)*((6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13)^4-(T_inf)^4), endfunction
I313=intg(53,100,i313)
disp(I313,'I313=')
 
/////////////I314/////////////////////////
function y=i314(t),y=(((10^(-12))*t^6 -4*(10^(-10))*t^5 +4*(10^(-8))*t^4 -2*(10^(-6))*t^3 +5*(10^(-5))*t^2 -0.0005*t +0.0012)*(T_fl-(6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13))*exp(-0.35*(-0.9484*t+ 103.98)/ L_fl)*((Uw/nu_g)*((-0.9484)*t+ 103.98))^(0.5))/(-0.9484*t+ 103.98), endfunction
I314=intg(53,100,i314)
disp(I314,'I314=')
 
/////////////I315/////////////////////////
function y=i315(t),y= ((10^(-12))*t^6 -4*(10^(-10))*t^5 +4*(10^(-8))*t^4 -2*(10^(-6))*t^3 +5*(10^(-5))*t^2 -0.0005*t +0.0012)*(T_b-(6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13))*exp(-0.25*s*(-0.9484*t+ 103.98)), endfunction
I315=intg(53,100,i315)
disp(I315,'I315=')
 
///////////////I_316////////////////////
function y=i316(t),y=(10^(-12))*t^6 -4*(10^(-10))*t^5 +4*(10^(-8))*t^4 -2*(10^(-6))*t^3 +5*(10^(-5))*t^2 -0.0005*t +0.0012, endfunction
I316=intg(53,100,i316)
disp(I316,'I316=')
 
////////////////I321///////////////////////
function y=i321(t),y=((10^(-12))*t^6 -4*(10^(-10))*t^5 +4*(10^(-8))*t^4 -2*(10^(-6))*t^3 +5*(10^(-5))*t^2 -0.0005*t +0.0012)*exp(-0.25*s*(-0.9484*t+ 103.98)), endfunction
I321=intg(53,100,i321)
disp(I321,'I321=')
 
////////////////I322 ///////////////////////
function y=i322(t),y=exp(-0.5*s*(-0.9484*t+ 103.98)),endfunction
I322=intg(53,100,i322)
disp(I322,'I322=')
 
///////////////I323 ///////////////////////
function y=i323(t),y=((6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13)^(4)-(T_inf)^(4))*exp(-0.5*s*(-0.9484*t+ 103.98)),
endfunction
I323=intg(53,100,i323)
disp(I323,'I323=')
 
///////////////I324 ///////////////////////
function y=i324(t), y=exp(-0.3*((-0.9484*t+ 103.98)/(L_fl))-0.25*s*(-0.9484*t+ 103.98))*((Uw/nu_g)*((-0.9484)*t+ 103.98))^(0.5)*( T_fl-(6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13))*(1/(-0.9484*t+ 103.98)), endfunction
I324=intg(53,100,i324) 
disp(I324,'I324=')
 
 
 
/////////////////I325 ///////////////////////
function y=i325(t),y= ((10^(-12))*t^6 -4*(10^(-10))*t^5 +4*(10^(-8))*t^4 -2*(10^(-6))*t^3 +5*(10^(-5))*t^2 -0.0005*t +0.0012)*(T_b-(6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13))*exp(-0.25*s*(-0.9484*t+ 103.98)), endfunction
I325=intg(53,100,i325)
disp(I325,'I325=')
 
 
/////////////////I326 ///////////////////////)
function y=i326(t),y=exp(-0.25*s*(-0.9484*t+ 103.98)), endfunction
I326=intg(53,100,i326)
disp(I326,'I326=')
 
/////////////////I331 ///////////////////////)
function y=i331(t),y=((10^(-12))*t^6 -4*(10^(-10))*t^5 +4*(10^(-8))*t^4 -2*(10^(-6))*t^3 +5*(10^(-5))*t^2 -0.0005*t +0.0012)*((6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13)^(4)-(T_inf)^4), endfunction
I331=intg(53,100,i331)
disp(I331,'I331=')
 
 
/////////////////I332 ///////////////////////((6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13)^(4)-(T_inf)^4)
function y=i332(t),y=exp(-0.25*s*(-0.9484*t+ 103.98))*((6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13)^(4)-(T_inf)^4), endfunction
I332=intg(53,100,i332)
disp(I332,'I332=')
 
/////////////////I333 ///////////////////////
function y=i333(t),y=(((6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13)^(4)-(T_inf)^4))^2 ,endfunction
I333=intg(53,100,i333)
disp(I333,'I333=')
 
 
 
/////////////////I334 ///////////////////////
function y=i334(t),y=((Uw/nu_g)*((-0.9484)*t+ 103.98))^(0.5)*(6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13-T_fl)*exp(-0.3*s*(-0.9484*t+ 103.98))*((6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13)^(4)-(T_inf)^4),endfunction
I334=intg(53,100,i334)
disp(I334,'I334=')
 
 
/////////////////I335 ///////////////////////
function y=i335(t),y=(T_b-(6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13))*exp(-0.25*s*(-0.9484*t+ 103.98))*((6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13)^(4)-(T_inf)^4),endfunction
I335=intg(53,100,i335)
disp(I335,'I335=')
 
/////////////////I336 ///////////////////////
function y=i336(t),y=((6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13)^(4)-(T_inf)^4),endfunction
I336=intg(53,100,i336)
disp(I336,'I336=')
 
/////////////////I341 ///////////////////////
function y=i341(t),y=(((10^(-12))*t^6 -4*(10^(-10))*t^5 +4*(10^(-8))*t^4 -2*(10^(-6))*t^3 +5*(10^(-5))*t^2 -0.0005*t +0.0012)*(T_fl-(6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13))*exp(-0.3*(-0.9484*t+ 103.98)/ L_fl)*((Uw/nu_g)*((-0.9484)*t+ 103.98))^(0.5))/(-0.9484*t+ 103.98), endfunction
I341=intg(53,100,i341)
disp(I341,'I341=')
 
 
///////////////////I342 //////////////////////
function y=i342(t),y=exp(-0.25*s*(-0.9484*t+ 103.98)-0.3*(-0.9484*t+ 103.98)/ L_fl)*((Uw/nu_g)*((-0.9484)*t+ 103.98))^(0.5)*(T_fl-(6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13))*(1/(-0.9484*t+ 103.98)), endfunction
 
I342=intg(53,100,i342)
 
disp(I342,'I342=')
//
///////////////////I343 //////////////////////
function y=i343(t),y=((6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13)^(4)-(T_inf)^4)*( T_fl-(6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13))*((Uw/nu_g)*((-0.9484)*t+ 103.98))^(0.5)*exp(-0.3*(-0.9484*t+ 103.98)/ L_fl), endfunction
I343=intg(53,100,i343)
disp(I343,'I343=')
 
 
///////////////////I344 //////////////////////
function y=i344(t),y=((Uw/nu_g)*((-0.9484)*t+ 103.98))*( T_fl-(6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13))^(2)*exp(-0.6*(-0.9484*t+ 103.98)/ L_fl)*(1/( -0.9484*t+ 103.98)^(2)), endfunction
I344=intg(53,100,i344)
disp(I344,'I344=')
 
///////////////I345 ///////////////////////
function y=i345(t), y=exp(-0.3*((-0.9484*t+ 103.98)/(L_fl))-0.25*s*(-0.9484*t+ 103.98))*((Uw/nu_g)*((-0.9484)*t+ 103.98))^(0.5)*( T_fl-(6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13))*(1/(-0.9484*t+ 103.98)), endfunction
I345=intg(53,100,i345) 
disp(I345,'I345=')
 
 
///////////////I346 ///////////////////////
function y=i346(t), y=exp(-0.3*((-0.9484*t+ 103.98)/(L_fl)))*((Uw/nu_g)*((-0.9484)*t+ 103.98))^(0.5)*( T_fl-(6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13))*(1/(-0.9484*t+ 103.98)), endfunction
I346=intg(53,100,i346) 
disp(I346,'I346=')
 
 
/////////////////I351 ///////////////////////
function y=i351(t),y= ((10^(-12))*t^6 -4*(10^(-10))*t^5 +4*(10^(-8))*t^4 -2*(10^(-6))*t^3 +5*(10^(-5))*t^2 -0.0005*t +0.0012)*(T_b-(6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13))*exp(-0.25*s*(-0.9484*t+ 103.98)), endfunction
I351=intg(53,100,i351)
disp(I351,'I351=')
 
 
 
/////////////////I352 ///////////////////////
function y=i352(t),y= ((10^(-12))*t^6 -4*(10^(-10))*t^5 +4*(10^(-8))*t^4 -2*(10^(-6))*t^3 +5*(10^(-5))*t^2 -0.0005*t +0.0012)*(T_b-(6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13))*exp(-0.5*s*(-0.9484*t+ 103.98)), endfunction
I352=intg(53,100,i352)
disp(I352,'I352=')
 
/////////////////I353 ///////////////////////
function y=i353(t),y=(T_b-(6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13))*exp(-0.25*s*(-0.9484*t+ 103.98))*((6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13)^(4)-(T_inf)^4),endfunction
I353=intg(53,100,i353)
disp(I353,'I353=')
 
 
///////////////I354 ///////////////////////
function y=i354(t), y=exp(-0.3*((-0.9484*t+ 103.98)/(L_fl))-0.25*s*(-0.9484*t+ 103.98))*((Uw/nu_g)*((-0.9484)*t+ 103.98))^(0.5)*( T_fl-(6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13))*(1/(-0.9484*t+ 103.98))*( T_b-(6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13)), endfunction
I354=intg(53,100,i354) 
disp(I354,'I354=')
 
 
///////////////I355 ///////////////////////
function y=i355(t), y=(T_b-(6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13))^(0.5)*exp(-0.5*s*(-0.9484*t+ 103.98)), endfunction
I355=intg(53,100,i355) 
disp(I355,'I355=')
 
 
 
///////////////I356 ///////////////////////
function y=i356(t), y=(T_b-(6*(10^-9)*t^6 - 2*(10^-6)*t^5 + 0.0002*t^4 - 0.0095*t^3 + 0.2239*t^2 - 2.1302*t + 310.13))*exp(-0.25*s*(-0.9484*t+ 103.98)), endfunction
I356=intg(53,100,i356) 
disp(I356,'I356=')
 
/////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
/////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
/////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
//////////////////////////////////////// SYSTEME NON LINEAIRE DE 3 EQUATIONS A 3 INCONNUES///////////////////////////////////////////////////////////////
////////////////////////////////////////////////////////////(T_fl; a_fb; C_pc)///////////////////////////////////////////////////////////////////////////////
A=esp_fl*sigma*(1/(rho_c*phi*H_c)) ; B= 0.25 * s * esp_b *sigma*(T_b)^4*( 1 /(rho_c * phi )) ; 
 
C= esp_c*sigma * ( 1 /(rho_c * phi * H_c)) ; D= 0.565* Pr^(0.5) * k_fl*( 1 /(rho_c * phi * H_c));
 
E= 0.911*s*(R_eD)^(0.385)*Pr^(1/3)*k_b*( 1 /(rho_c * phi * H_c)); 
 
// EQUATION 1// 
 
// A^2*a_fb*(T_fl)^(8)*(1/(C_pc)^2)*I11 + (T_fl)^(4)*A*B*(1/(C_pc)^2)*I12 - A*(T_fl)^(4)*C*(1/(C_pc)^2)*I13+ (T_fl)^(4)*D*A*(1/(C_pc)^2)*I14 + E*A*(T_fl)^(4)*(1/(C_pc)^2)*I15 + A*(T_fl)^(4)*dt_T_exp*(1/(C_pc)^2) == 0
 
// équation 1 après changement de variable
 
// A^2*(x(1))*(x(2))^(8)*(1/(x(3))^2)*I11 + (x(2))^(4)*A*B*(1/(x(3))^2)*I12 - A*(x(2))^(4)*C*(1/(x(3))^2)*I13+ (x(2))^(4)*D*A*(1/(x(3))^2)*I14 + E*A*(x(2))^(4)*(1/(x(3))^2)*I15 + A*(x(2))^(4)*dt_T_exp*(1/(x(3)^2) == 0
 
 
// EQUATION 2//
 
//4*A^(2)*(a_fb)^(2)* (T_fl)^(7)*(1/(C_pc)^2)*I21 + 4*A*B*(a_fb)*(T_fl)^(3)*(1/(C_pc)^2)*I22 - 4*A*C*a_fb *(T_fl)^(3)*(1/(C_pc)^2)*I23+ 4*A*D*a_fb *(T_fl)^(3)*(1/(C_pc)^2)*I24 + 4*A*E*a_fb *(T_fl)^(3)*(1/(C_pc)^2)*I25 + A*D*a_fb *(T_fl)^(4)*(1/(C_pc)^2)*I26 + B*D*(1/(C_pc)^2)*I27 - C*D*(1/(C_pc)^2)*I28 + D^2*(1/(C_pc)^2)* I29 + E*D*(1/(C_pc)^2)*I210 -D* dt_T_exp* (1/C_pc)-4*A*a_fb*(T_fl)^(3)*dt_T_exp*(1/C_pc) *I212 == 0
 
// équation 2 après changement de variable
 
//4*A^(2)*(x(1))^(2)* (x(2))^(7)*(1/(x(3))^2)*I21 + 4*A*B*(x(1))*(x(2))^(3)*(1/(x(3))^2)*I22 - 4*A*C*a_fb *(x(2))^(3)*(1/(x(3)^2)*I23+ 4*A*D*a_fb *(x(2))^(3)*(1/(x(3))^2)*I24 + 4*A*E*a_fb *(x(2))^(3)*(1/(x(3))^2)*I25 + A*D*a_fb *(x(2))^(4)*(1/(x(3))^2)*I26 + B*D*(1/(x(3))^2)*I27 - C*D*(1/(x(3))^2)*I28 + D^2*(1/(x(3))^2)* I29 + E*D*(1/(x(3))^2)*I210 -D* dt_T_exp* (1/(x(3)))-4*A*a_fb*(x(2))^(3)*dt_T_exp*(1/(x(3))) *I212 == 0
 
 
 
// EQUATION 3//
 
// (a_fb)^(2)* (T_fl)^(8)*A^(2)*(1/(C_pc)^3)*I311 - a_fb*(T_fl)^(4) *B*(1/(C_pc)^3)*I312 + a_fb*(T_fl)^(4) *C*(1/(C_pc)^3)*I313 - a_fb*(T_fl)^(4) *D*(1/(C_pc)^3)*I314 - a_fb*(T_fl)^(4) *E*(1/(C_pc)^3)*I315 + a_fb*(T_fl)^(4) *dt_T_exp*(1/(C_pc)^2)*I316 - a_fb*(T_fl)^(4) *A*B*(1/(C_pc)^3)*I321 - B^2* (1/(C_pc)^3)*I322 + B*C*(1/(C_pc)^3)*I323 - B*D*(1/(C_pc)^3)*I324 - B*E*(1/(C_pc)^3)*I325 + B*dt_T_exp*(1/(C_pc)^2)* I326 + A*a_fb*(T_fl)^(4)*C*(1/(C_pc)^3)*I331 + B*C*(1/(C_pc)^3)*I332 - C^2*(1/(C_pc)^3)*I333 + D*C*(1/(C_pc)^3)*I334 + E*C*(1/(C_pc)^3)*I335 +C*dt_T_exp*(1/(C_pc)^2)*I336 - A*D*a_fb*(T_fl)^4*(1/(C_pc)^3)*I341 - B*D *(1/(C_pc)^3)*I342 + C*D*(1/(C_pc)^3)*I343 - D^2*(1/(C_pc)^3)*I344 - E*D*(1/(C_pc)^3)*I345 + dt_T_exp*D*(1/(C_pc)^2)*I346 - A*a_fb*(T_fl)^4*E*(1/(C_pc)^3)*I351 - E*B*(1/(C_pc)^3)*I352+ C*E*(1/(C_pc)^3)*I353 - D*E*(1/(C_pc)^3)*I354- E^2*(1/(C_pc)^3)*I355 + dt_T_exp*E*(1/(C_pc)^2)*I356 == 0
 
// équation 3 après changement de variable
// (x(1))^(2)* (x(2))^(8)*A^(2)*(1/(x(3))^3)*I311 - a_fb*(x(2))^(4) *B*(1/(x(3))^3)*I312 + a_fb*(x(2))^(4) *C*(1/(x(3))^3)*I313 - a_fb*(x(2))^(4) *D*(1/(x(3))^3)*I314 - a_fb*(x(2))^(4) *E*(1/(x(3))^3)*I315 + a_fb*(x(2))^(4) *dt_T_exp*(1/(x(3))^2)*I316 - (x(1))*(T_fl)^(4) *A*B*(1/(x(3))^3)*I321 - B^2* (1/(x(3))^3)*I322 + B*C*(1/(x(3))^3)*I323 - B*D*(1/(x(3))^3)*I324 - B*E*(1/(x(3))^3)*I325 + B*dt_T_exp*(1/(x(3))^2)* I326 + A*(x(1))*(x(2))^(4)*C*(1/(x(3))^3)*I331 + B*C*(1/(x(3))^3)*I332 - C^2*(1/(x(3))^3)*I333 + D*C*(1/(x(3))^3)*I334 + E*C*(1/(x(3))^3)*I335 +C*dt_T_exp*(1/(x(3))^2)*I336 - A*D*(x(1))*(x(2))^4*(1/(x(3))^3)*I341 - B*D *(1/(x(3))^3)*I342 + C*D*(1/(x(3))^3)*I343 - D^2*(1/(x(3))^3)*I344 - E*D*(1/(x(3))^3)*I345 + dt_T_exp*D*(1/(x(3))^2)*I346 - A*(x(1))*(x(2))^4*E*(1/(x(3))^3)*I351 - E*B*(1/(x(3))^3)*I352+ C*E*(1/(x(3))^3)*I353 - D*E*(1/(x(3))^3)*I354- E^2*(1/(x(3))^3)*I355 + dt_T_exp*E*(1/(x(3))^2)*I356 == 0
 
//deff('[y]=f(x)','y=[-5*x(1)+2*(x(1))+(x(2))^2;-5*x(2)+4*x(1)+2*(x(2)) ]');
//X=fsolve([1,1],f)
//A^2*x(1)*x(2)^(8)*(1/x(3)^2)*I11 +(x(2)^4)*A*B*(1/x(3)^2)*I12 -A*(x(2)^4)*C*(1/x(3)^2)*I13+ x(2)^(4)*D*A*(1/x(3)^2)*I14 + E*A*(x(2)^4)*(1/x(3)^2)*I15 + A*(x(2)^4)*(dt_T_exp)*(1/(x(3)^2)
//4*(A^2)*(x(1)^2)* (x(2)^7)*(1/x(3)^2)*I21 + 4*A*B*x(1)*(x(2)^3)*(1/x(3)^2)*I22 - 4*A*C*x(1)*(x(2)^3)*(1/(x(3)^2)*I23+ 4*A*D*x(1)*(x(2)^3)*(1/x(3)^2)*I24 + 4*A*E*x(1)*(x(2)^3)*(1/x(3)^2)*I25 + A*D*x(1)*(x(2)^4)*(1/x(3)^2)*I26 + B*D*(1/x(3)^2)*I27 - C*D*(1/x(3)^2)*I28 + D^2*(1/x(3)^2)* I29 + E*D*(1/x(3)^2)*I210 -D* (dt_T_exp)* (1/x(3))-4*A*x(1)*(x(2)^3)*(dt_T_exp)*(1/x(3)) *I212 
 
deff('y=f(x)','y=[A^2*x(1)*x(2)^(8)*(1/x(3)^2)*I11 +(x(2)^4)*A*B*(1/x(3)^2)*I12 -A*(x(2)^4)*C*(1/x(3)^2)*I13+ x(2)^(4)*D*A*(1/x(3)^2)*I14 + E*A*(x(2)^4)*(1/x(3)^2)*I15 + A*(x(2)^4)*(dt_T_exp)*(1/(x(3)^2);4*(A^2)*(x(1)^2)* (x(2)^7)*(1/x(3)^2)*I21 + 4*A*B*x(1)*(x(2)^3)*(1/x(3)^2)*I22 - 4*A*C*x(1)*(x(2)^3)*(1/(x(3)^2)*I23+ 4*A*D*x(1)*(x(2)^3)*(1/x(3)^2)*I24 + 4*A*E*x(1)*(x(2)^3)*(1/x(3)^2)*I25 + A*D*x(1)*(x(2)^4)*(1/x(3)^2)*I26 + B*D*(1/x(3)^2)*I27 - C*D*(1/x(3)^2)*I28 + D^2*(1/x(3)^2)* I29 + E*D*(1/x(3)^2)*I210 -D* (dt_T_exp)* (1/x(3))-4*A*x(1)*(x(2)^3)*(dt_T_exp)*(1/x(3)) *I212 ;x(1)^(2)* x(2)^(8)*A^(2)*(1/x(3)^3)*I311 - x(1)*x(2)^(4) *B*(1/x(3)^3)*I312 + x(1)*(x(2)^4) *C*(1/x(3)^3)*I313 - x(1)*(x(2)^4) *D*(1/x(3)^3)*I314 - x(1)*(x(2)^4) *E*(1/x(3)^3)*I315 +x(1)*(x(2)^4) *(dt_T_exp)*(1/x(3)^2)*I316 - x(1)*(x(2)^4) *A*B*(1/x(3)^3)*I321 - B^2* (1/x(3)^3)*I322 + B*C*(1/x(3)^3)*I323 - B*D*(1/x(3)^3)*I324 - B*E*(1/x(3)^3)*I325 + B*(dt_T_exp)*(1/x(3)^2)* I326 + A*x(1)*(x(2)^4)*C*(1/x(3)^3)*I331 + B*C*(1/x(3)^3)*I332 - C^2*(1/x(3)^3)*I333 + D*C*(1/x(3)^3)*I334 + E*C*(1/x(3)^3)*I335 +C*(dt_T_exp)*(1/x(3)^2)*I336 - A*D*x(1)*x(2)^4*(1/x(3)^3)*I341 - B*D *(1/x(3)^3)*I342 + C*D*(1/x(3)^3)*I343 - D^2*(1/x(3)^3)*I344 - E*D*(1/x(3)^3)*I345 + (dt_T_exp)*D*(1/x(3)^2)*I346 - A*x(1)*x(2)^4*E*(1/x(3)^3)*I351 - E*B*(1/x(3)^3)*I352+ C*E*(1/x(3)^3)*I353 - D*E*(1/x(3)^3)*I354- (E^2)*(1/x(3)^3)*I355 + (dt_T_exp)*E*(1/x(3)^2)*I356 ]');
X=fsolve([0.5,1083,1.25*10^3],f)

aider moi s il vous plait c est vraiment urgent

**Jerome Briot** · 21/11/2017, 12h46

À toi de simplifier les parenthèses inutiles. Il y en a beaucoup trop dans ton code.

Par exemple nu_g = 7.77*10^(-5) s'écrit sans les parenthèses :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

nu_g = 7.77*10^-5

Voir encore plus simplement :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

nu_g = 7.77E-5

Ou encore y = (Uw/nu_g)*((-0.9484)*t+ 103.98) qui s'écrit plus simplement :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

y = Uw/nu_g*(-0.9484*t+ 103.98)

Et pour finir A^2*x(1)*x(2)^(8)*(1/x(3)^2)*I11 +(x(2)^4)*A*B*(1/x(3)^2)*I12 -A*(x(2)^4)*C*(1/x(3)^2)*I13 qui s'écrit :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

A^2*x(1)*x(2)^8/x(3)^2*I11 + x(2)^4*A*B/x(3)^2*I12 -A*x(2)^4*C/x(3)^2*I13

**tchiekre michel** · 21/11/2017, 12h56

merci beaucoup pour la reponse mais mon probleme se trouve au bas de mon code au niveau du systeme non lineaire. (comme je l ai dit tantot je debute avec scilab j ai donc peur de mal programmer mes equation d ou le nombre de parenthese.)

voici mon système non linéaire

deff('y=f(x)','y=[A^2*x(1)*x(2)^(8)*(1/x(3)^2)*I11 +(x(2)^4)*A*B*(1/x(3)^2)*I12 -A*(x(2)^4)*C*(1/x(3)^2)*I13+ x(2)^(4)*D*A*(1/x(3)^2)*I14 + E*A*(x(2)^4)*(1/x(3)^2)*I15 + A*(x(2)^4)*(dt_T_exp)*(1/(x(3)^2);4*(A^2)*(x(1)^2)* (x(2)^7)*(1/x(3)^2)*I21 + 4*A*B*x(1)*(x(2)^3)*(1/x(3)^2)*I22 - 4*A*C*x(1)*(x(2)^3)*(1/(x(3)^2)*I23+ 4*A*D*x(1)*(x(2)^3)*(1/x(3)^2)*I24 + 4*A*E*x(1)*(x(2)^3)*(1/x(3)^2)*I25 + A*D*x(1)*(x(2)^4)*(1/x(3)^2)*I26 + B*D*(1/x(3)^2)*I27 - C*D*(1/x(3)^2)*I28 + D^2*(1/x(3)^2)* I29 + E*D*(1/x(3)^2)*I210 -D* (dt_T_exp)* (1/x(3))-4*A*x(1)*(x(2)^3)*(dt_T_exp)*(1/x(3)) *I212 ;x(1)^(2)* x(2)^(8)*A^(2)*(1/x(3)^3)*I311 - x(1)*x(2)^(4) *B*(1/x(3)^3)*I312 + x(1)*(x(2)^4) *C*(1/x(3)^3)*I313 - x(1)*(x(2)^4) *D*(1/x(3)^3)*I314 - x(1)*(x(2)^4) *E*(1/x(3)^3)*I315 +x(1)*(x(2)^4) *(dt_T_exp)*(1/x(3)^2)*I316 - x(1)*(x(2)^4) *A*B*(1/x(3)^3)*I321 - B^2* (1/x(3)^3)*I322 + B*C*(1/x(3)^3)*I323 - B*D*(1/x(3)^3)*I324 - B*E*(1/x(3)^3)*I325 + B*(dt_T_exp)*(1/x(3)^2)* I326 + A*x(1)*(x(2)^4)*C*(1/x(3)^3)*I331 + B*C*(1/x(3)^3)*I332 - C^2*(1/x(3)^3)*I333 + D*C*(1/x(3)^3)*I334 + E*C*(1/x(3)^3)*I335 +C*(dt_T_exp)*(1/x(3)^2)*I336 - A*D*x(1)*x(2)^4*(1/x(3)^3)*I341 - B*D *(1/x(3)^3)*I342 + C*D*(1/x(3)^3)*I343 - D^2*(1/x(3)^3)*I344 - E*D*(1/x(3)^3)*I345 + (dt_T_exp)*D*(1/x(3)^2)*I346 - A*x(1)*x(2)^4*E*(1/x(3)^3)*I351 - E*B*(1/x(3)^3)*I352+ C*E*(1/x(3)^3)*I353 - D*E*(1/x(3)^3)*I354- (E^2)*(1/x(3)^3)*I355 + (dt_T_exp)*E*(1/x(3)^2)*I356 ]');
X=fsolve([0,1083,1.25*10^3],f)

problème parenthèse de mon code scilab pour résoudre un système d’équations non linéaires

Scilab

Discussions similaires

Partager

Partager