Re: Values of V and abm

[asha g]

Yes the end closes the for loop. 
I wanted abmm to store all values of the iteration - not just the last one.  So 
I wrote that. Maybe there is a better way of writing it.
But more importantly, what is the bit abt not passing the arguments?
I am sending the entire script file by attachment. But as there are many 
function files to it- it might still be hard to debug.
thanks
asha 

> > I am having some problems in a program I have written.
> It is pretty long. But the section is :
> >
> 
> It's usually not much useful posting code that
> can't be used to
> reproduce the bug.
> 
> > N = 10;
> > niter = 10000;
> > for iter = 1: niter
> >
> > [P,S] = timedepcabactoct(Q,V,h2,Iion,deltaT);
> > for jj = 2:N-1
> > V(jj) = S(jj);
> > end
> > if deltaT*niter < 5
> > I = 0;
> > elseif 5 == deltaT*niter <= 55
> > I = 48.4*10^-12;
> > V(1) = V(2)+(Ri/Ao)*I*hd;
> > V(N) = 4/3*V(N-1)-(1/3)*V(N-2);
> > end
> > disp(V)
> > calling function ratcon2oct
> > [alphan,betan,alpham,betam,alphah,betah] =
> ratcon2oct(V);
> > abm = alpham +betam;
> > m = m+(deltaT.*(alpham -m.*abm));
> > abh = alphah + betah;
> > h = h +(deltaT.*(alphah-h.*abh));
> > abn = alphan + betan;
> > n = n + (deltaT.*(alphan-n.*abn));
> > end
> 
>     ^ If I've followed the code correctly, then this
> end corresponds
> to the for iter=... loop.
> Therefore, iter in the following statements has the value
> niter. Is
> this intentional?
> 
> > abmm(iter,:)= abm;
> > abnn(iter,:)= abn;
> > abhh(iter,:)= abh;
> > abm1= abmm(iter,1);
> > t = (1:niter)*deltaT;
> > plot(t,abm1)
> >
> > warning: implicit conversion from real matrix to real
> scalar
> >
> 
> If this is warning from your oct-file, then you're
> probably passing in
> wrong arguments.
> Again, it is hard to tell unless you give enough code to
> reproduce the bug.
> 
> >
> > I need to get an oscillating plot of abm1, abn1 and
> abh1 vs t. But I don't. Instead I get the warning given
> above which I don't follow. Also, if I try to check V at
> the point indicated I get the vector V repeated (N-2) times
> before the values change. I would have expected it to
> change at every iter count. Please tell me if I am writing
> the loop incorrectly.
> >
> > Thanks
> > Asha
> >
> >
> >
> > _______________________________________________
> > Help-octave mailing list
> > address@hidden
> > https://www.cae.wisc.edu/mailman/listinfo/help-octave
> >
> 
> 
> 
> -- 
> RNDr. Jaroslav Hajek
> computing expert
> Aeronautical Research and Test Institute (VZLU)
> Prague, Czech Republic
> url: www.highegg.matfyz.cz

cabunbact2oct.m
Description: Text Data