PS: cambiamo un attimo argomento. Non sono riuscito a far funzionare la funzione qr:
- Codice:
qr =@(x) qsun*exp(-beta*x);
Io ho scritto questo codice:
- Codice:
tspan=[0 L];
y_0=[Tf1,2000,q_sun,0];
fun=@(t,y) [h*(y(2)-y(1))/(rho_u*cp); %bilancio al fluido
y(4); %creazione variabile per la derivata di Ts
-beta_*y(3); %derivata della qr
1/ks*((-beta_*y(3))-h*(y(2)-y(1)))]; %bilancio al solido
come mi hai detto tu, e funziona.
Poi ho provato a scrivere:
- Codice:
y_0=[Tf1,2000,0];
qr =@(t) q_sun*exp(-beta_*t);
fun=@(t,y) [h*(y(2)-y(1))/(rho_u*cp); %bilancio al fluido
y(3); %creazione variabile per la derivata di Ts
1/ks*qr-h*(y(2)-y(1))]; %bilancio al solido
Ma mi da errore.