error var(x)???

[address@hidden]

i have an error in the file var.m tpyping "runtesi" using this m-files:

runtesi.m:
# sistema adimensionale della tesi master
# controllare bene i valori numerici
# per quelli mancanti utilizzo i valori del Nylon, come da tesi master
# il file runtesi dà i valori iniziali L,v0,T0,Ta al fine tesi.m
# tutte le unità di misura sono mks
# i valori iniziali sono adimensionali!!!
# dovrei utilizzare una condizione al contorno per la velocità!!!!
# invece sono costretto a usare la condizione per tau0, ovviamente 
improvvisato!!!!
# y(:,1)=v(z)   asterisco
# y(:,2)=T(z)   asterisco
# y(:,3)=tau(z) asterisco
L=1;
T0=520;
Ta=17+273;
v0=3.4;
y0=[1 1 10];
z=linspace(0,L,100)';
y=lsode("tesi",y0,z);
plot(z,y(:,2))
#plot(z,y(:,1))
#plot(z*L,(y(:,2)*(T0-Ta)+Ta)
#plot(z*L,y(:,1)*v0)

tesi.m:
function tesi=f(y,z)
# sistema adimensionale della tesi master
# controllare bene i valori numerici
# per quelli mancanti utilizzo i valori del Nylon, come da tesi master
# il file runtesi dà i valori iniziali L,v0,T0,Ta al fine tesi.m
# tutte le unità di misura sono mks
# y(1)=v(z)   asterisco
# y(2)=T(z)   asterisco
# y(3)=tau(z) asterisco
L=1;
T0=520;
Ta=17+273;
v0=3.4;
R0=5.3e-04;
rho=900;
rhoa=1.189;
Cp=2500;
a=0.37;
b=0.61;
va=0;
vc=1;
eta0=60;
Ea=13500*4.1855;
G=3*1.1e05;
Nua=1.819e-05;
Re=2*(y(1)-va)*rhoa/Nua;
Cd=a*Re^b;
R=sqrt(1/y(1));
Ka=0.2567;
muc=exp((Ea/R)*((1/y(2))-(1/297)));
Etap=eta0*exp((Ea/R)*((1/y(2))-(1/297)));
alfa=(1/R)*0.21*Re^(1/3)*Ka*(1+(8*vc/y(1)-va)^2)^(1/6);
C0=eta0/(rho*L*v0);
C1=rhoa*Cd*L/(rho*R0);
C2=rho*L/v0^2;
C3=eta0*v0/(rho*Cp*L*(T0-Ta));
C4=2*alfa*L/(rho*Cp*R0*v0);
C5=eta0*v0/(G*L);
F1=1-(C0*(1+2*C5*y(1)))/(C5*((y(1)^2)*C0*y(3)));
tesi(1)=((1/F1)*(1+2*C5*y(1)/(C5*y(1)))*(-(y(1)^2/(y(1)^2+C0*y(3)))*C1*R*(y(1)-var)^2+(y(1)^2/(y(1)^2+C0*y(3)))*C2*R^2)-y(3)/(C5*muc*y(1)*F1))*(C0*y(1)/(y(1)^2+C0*y(3)))-(y(1)^2/(y(1)^2+C0*y(3)))*C1*R*(y(1)-var)^2+(y(1)^2/(y(1)^2+C0*y(3)))*C2*R^2;
tesi(2)=C3*y(3)*tesi(1)/y(1)-C4*y(2)/(y(1)*R);
tesi(3)=(1/F1)*(1+2*C5*y(1)/(C5*y(1)))*(-(y(1)^2/(y(1)^2+C0*y(3)))*C1*R*(y(1)-var)^2+(y(1)^2/(y(1)^2+C0*y(3)))*C2*R^2)-y(3)/(C5*muc*y(1)*F1);
endfunction


what's the meaning of that error?
what can i do?

thanks all


________________________________________________________________________
Cerchi un laboratorio fotografico aperto 24 ore su 24?
Stampa le tue foto digitali su Kataweb e le ricevi a domicilio in 48 ore.
http://www.kataweb.it/foto




-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------