[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: GNU plot
From: |
asha g |
Subject: |
Re: GNU plot |
Date: |
Mon, 21 Apr 2008 00:01:36 -0700 (PDT) |
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
- Re: GNU plot, (continued)
- Re: GNU plot, Tatsuro MATSUOKA, 2008/04/21
- Re: GNU plot, asha g, 2008/04/21
- Re: GNU plot, Sergei Steshenko, 2008/04/21
- Re: GNU plot, Tatsuro MATSUOKA, 2008/04/21
- Re: GNU plot, Sergei Steshenko, 2008/04/21
- Re: GNU plot, asha g, 2008/04/22
- Re: GNU plot, John W. Eaton, 2008/04/22
- Re: gnuplot, asha g, 2008/04/25
- Re: gnuplot, Ben Abbott, 2008/04/25
- Re: gnuplot, Sergei Steshenko, 2008/04/25
Re: GNU plot,
asha g <=