Re: GNU plot

[asha g]

Hi Tatsuro,
Thanks for volunteering to help. Well, I have got a
few to run. so maybe it is program specific. I am
using GNUPlot 4.0, RHEL 4.0.
I am sending one scriptfile which is showing the error
messages I described. Will it do just to have the
script file/ Won't you need the function files to run
it? 
thanks for your help
Asha 


--- Tatsuro MATSUOKA <address@hidden> wrote:

> Hello
> 
> Information is not adequate so no one can reply your
> inquiry.
> Please show the your octave script, gnuplot version
> and platform use.
> 
> From 2.9.17?, the handing of gnuplot is changed.
> 
> And now gnuplot version is higher than 4.2.2 is
> required.
> 
> Please give us more detailed infomation if you are
> in replyer stand point.
> 
> Can you answer the question with only these
> information?
> 
> Regards
> 
> Tatsuro 
> 
> --- asha g <address@hidden> wrote:
> 
> > I had mentioned that I just reloaded octave 3.0
> and am
> > trying to run some of my old files written in
> 2.9.12
> > but am getting error messages of the sort:
> > 
> > gnu plot 0.60606060606 1.19184  e-14
> > line 0 invalid command
> > 
> > Please explain what is going on. Does it mean that
> my
> > files are not compatible. All these were plotting
> > earlier.
> > 
> > Thanks
> > Asha 
> > 
> >  
> >    
> >    
> >    
> >   Goals too clearly defined can become blinkers. 
> >   Mary Catherine Bateson 
> > 
> > 
> > 
> >      
>
____________________________________________________________________________________
> > Be a better friend, newshound, and 
> > know-it-all with Yahoo! Mobile.  Try it now. 
> >
>
http://mobile.yahoo.com/;_ylt=Ahu06i62sR8HDtDypao8Wcj9tAcJ
> > _______________________________________________
> > Help-octave mailing list
> > address@hidden
> >
>
https://www.cae.wisc.edu/mailman/listinfo/help-octave
> > 
> 
> 
> --------------------------------------
> GANBARE! NIPPON! Win your ticket to Olympic Games
> 2008.
> http://pr.mail.yahoo.co.jp/ganbare-nippon/
> 

 
   
   
   
  Goals too clearly defined can become blinkers. 
  Mary Catherine Bateson 



      
____________________________________________________________________________________
Be a better friend, newshound, and 
know-it-all with Yahoo! Mobile.  Try it now.  
http://mobile.yahoo.com/;_ylt=Ahu06i62sR8HDtDypao8Wcj9tAcJ

% Script file : cabranpa22sevoltanalyoct.m
% Purpose : To study V/Vo vs X for soma - dendrite branched case, two branches, 
passive, sealed end at both branches
% Current injection at x = 0
% Record of revisions:
% Date                                     Programmer                           
    Description
% March 9, 07                                Asha G                             
     original code
% July 2, 07                                Asha G                              
   modified for OCTAVE and also change VV to V
% Define variables
% N -- discretization by space
% V -- output volts at points x
% l -- length of dendrite in mm
% lambda -- space constant in mm
% L -- length nondimensionalised - l/lambda
% h -- deltax = l/N
% h2 -- h*h
% f -- R.H.S of compact diff scheme
% alpha -- coeff of L.H.S compact diff scheme
% alpha1 -- coeff of R.H.S compact diff scheme
% M -- tridiag matrix
% Q -- double derivative solved by M\f'
% k -- index variable 

% Calling function voltoct_cabranpap

clear all 
% clear global Rm Ri t tau I N 
pack ( ' cabranpa22sevoltanalyoct')
% declare global variables
global N   Rm Ri t tau I  Cm rinp
Rm= 20*10^5; Ri = 20*10^2; t= 1; tau = 10; I = 1* 10^-5; N = 10;Cm = 1*10^-2; 
rinp = 1* 10^6;
lp = 300*10^-3;
dp = 10*10^-3;
lambdap = ((Rm/Ri)*(dp/4))^1/2;
ra = (4*Ri)/(pi*dp^2);
rm = Rm /(pi*dp);
Vo = -I * rinp;

% calling function voltoct_cabranpap
[VVp,Vp,hp,hp2] = voltoct_cabranpap(lp,lambdap,Vo);

% calling function voltoct_cabranpab1
lb1 = 100* 10^-3;
db1 = 5*10^-3;
lambdab1 = ((Rm/Ri)*(db1/4))^1/2;
rab1 = (4*Ri)/(pi*db1^2);
rmb1 = Rm/(pi*db1);

[VVb1,Vb1,hb1,hb12]= volt_cabranpab1(VVp,lp,lb1,lambdab1,lambdap);
% calling function voltoct_cabranpab2
lb2 = 200* 10^-3;
db2 = 5* 10^-3;
lambdab2 = ((Rm/Ri)* (db2/4))^1/2;
rab2 = (4*Ri)/(pi*db2^2);
rmb2 = Rm/(pi *db2);

[VVb2,Vb2,hb2,hb22]= volt_cabranpab2(VVp,lp,lb2,lambdab2,lambdap);

for n = 1: 1125
% Calling function compactdiffcaboctp

[fp] = compactdiffcaboctp(VVp,hp2);

% Calling function tridiagdoubleoctp

alpha = 2/11;
alpha1 = 1/10;

[Qp,M] = tridiagdoubleoctp(alpha,alpha1,fp);

% Calling function timedepcabranvoltoctp
gamma = 1;

[Pp,Sp] = timedepcabranvoltoctp(Qp,VVp,hp2,gamma );

for jj = 2: N-1
VVp(jj) = Sp(jj);
end 
VVp(1) = 1;

% section for branch 1
% calling function compactdiffcaboctb1

[fb1] = compactdiffcaboctb1(VVb1,hb12);

% calling function tridiagdoubleoctb1

alpha = 2/11;
alpha1 = 1/10;

[Qb1,M] = tridiagdoubleoctb1(alpha,alpha1,fb1);

% Calling function timedepcabranvoltoctb1

gamma = 1;
ra = (4* Ri)/(pi*dp^2);
rm = Rm / (pi*dp);
rab1 = (4*Ri)/(pi*db1^2);
rab2 = (4*Ri) / (pi * db2^2);

[Pb1,Sb1] = timedepcabranvoltoctb1 ( Qb1,VVb1, hp2,gamma, lambdab1,lambdap);

for jj = 2:N-1

VVb1(jj) = Sb1(jj);

end 

VVb1(N) = (4/3) * VVb1(N-1) - ( 1/3) * VVb1(N-2);


% section for branch 2
% calling function compactdiffcaboctb2

[fb2] = compactdiffcaboctb2( VVb2,hb22);

% calling function tridiagdoubleoctb2

alpha = 2/11;
alpha1 = 1/10;

[ Qb2,M ] = tridiagdoubleoctb2( alpha, alpha1, fb2);

% calling function timedepcabranvoltoctb2 
gamma = 1; 
ra = ( 4* Ri) / ( pi * dp^2);
rm = Rm / (pi* dp);
rab1 = ( 4* Ri)/ (pi*db1^2);
rab2 = (4*Ri)/(pi*db2^2);
[pb2,Sb2] = timedepcabranvoltoctb2(Qb2,VVb2,hp2,gamma,lambdab2,lambdap);

for jj = 2:N-1
VVb2(jj) = Sb2(jj); 
end 

VVb2(N) = ( 4/3)* VVb2(N-1) - (1/3)* VVb2(N-2);

% calling function branchpointanalyoct

[Pbp,Sbp] = 
branchpointanalyoct(VVp,VVb1,VVb2,hp2,hb12,hb22,gamma,dp,db1,db2,rm,ra,rab1,rab2,rmb1,rmb2,lp,lb1,lb2);
VVp(N) = Sbp(N);
VVb1(1) = VVp(N); 
VVb2(1) = VVp(N);

nnnn(n,:) = VVp;
oooo(n,:)= VVb1;
pppp(n,:) = VVb2;

end 

% Calling function file voltoct1_cabranpase

%[VVa,Vo] = voltoct1_cabranpase( lp,lambdap,lb1,ra,dp);

hold on 

x = linspace(0,1,N);
X = x/lambdap;

x1 = linspace(1,1+lb1/lp,N);
X1 = (( x1/lambdap));

x2 = linspace(1,1+lb2/lp,N);
X2 = (( x2/lambdap));
% Xt2 = [X(N) + X2];
%x4 = linspace (0,1+lb1/lp,N);
%X4 = (( x4/lambdap));

Vo = -I*rinp;

y = VVp;
y1 = VVb1;
% yt = [y y1]; 
y2 = VVb2;
%y4 = VVa;

% yt2 = [y y2];
%x3 = linspace(0,t,n);
%y3 = nnnn(1:n,N);
plot(X,y,'rs')
plot(X1,y1,'m*')
plot (X2,y2,'g+')
%plot(X4,y4,'yd')
% plot (x3,y3)

xlabel('X')
ylabel('VV')
legend ("on")
legend ('rs==parent',' m*==branch1','g+ == branch2' )
title('branched, lb1<lb2, db1=db2')
print ('cabranpa22sevoltanalyoct10c.eps','-depsc')
hold off