diff options
Diffstat (limited to 'texmf-dist/doc/support/ketcindy/ketsample/samples/s07slides/fig/.Rhistory')
-rw-r--r-- | texmf-dist/doc/support/ketcindy/ketsample/samples/s07slides/fig/.Rhistory | 238 |
1 files changed, 238 insertions, 0 deletions
diff --git a/texmf-dist/doc/support/ketcindy/ketsample/samples/s07slides/fig/.Rhistory b/texmf-dist/doc/support/ketcindy/ketsample/samples/s07slides/fig/.Rhistory new file mode 100644 index 00000000..8470a278 --- /dev/null +++ b/texmf-dist/doc/support/ketcindy/ketsample/samples/s07slides/fig/.Rhistory @@ -0,0 +1,238 @@ +?dbeta +plot(x,dbeta(2,5,x) +) +?dbeta +plot(x,dbeta(x,2,5)) +x=seq(o,1,by=0.02) +x=seq(0,1,by=0.02) +plot(x,dbeta(x,2,5)) +m=(2)/(2+5) +s2=(2*5)/((2+5)^2*(2+5+1)) +m +s2 +plot(x,dbeta(x,2,5)) +sample(1:20, 10) +tmp=sample(1:20, 10) +sort(tmp) +tmp=sample(1:37,19) +tmp=sample(1:37,19); sorrt(tmp) +tmp=sample(1:37,19); sort(tmp) +search() +library(deSolve) +search() +library() +cat("\\ab\\\\") +samples(40,20) +sample(40,20) +source('~/ketcindy/ketfiles/ketwork/Dotfillrdf1.r', chdir = TRUE) +Tmp="Num=50" +K=strsplit(Tmp,'=') +K +substring(K[1],1) +part(K[1],1) +K[1] +substring(K[[1]][1],1) +substring(K[[1]][1],1,1) +K=K[[1]] +eval(parse(text=K[2])) +Num=eval(parse(text=K[2])) +length(Num) +?seq +seq(1,length(list(2,4,5))) +v=seq(1,length(list(2,4,5))) +v[1:2] +source('~/Dropbox/kettoday/0930bezierR/s0401bezierbasic.r', chdir = TRUE) +windispg(bza,bzb) +Windisp(bza,bzb) +Windisp(bza,bz1) +bza +source('~/Dropbox/kettoday/0930bezierR/s0401bezierbasic.r', chdir = TRUE) +bza +Windisp(bza,bz1) +source('~/Dropbox/kettoday/0930bezierimplicitR/s0212implicitimp1.r', chdir = TRUE) +Makevaltable<- function(MS){
Eps=10^(-6)
Nargs=length(MS)
Fstr=Op(1,MS)
Tmp=Op(2,MS)
Tmp1=strsplit(Tmp,'=',fixed=TRUE)
Xname=Tmp1[1]
Xrange=eval(parse(text=Tmp1[2]))
Tmp=Op(3,MS)
Tmp1=strsplit(Tmp,'=',fixed=TRUE)
Yname=Tmp1[1]
Yrange=evstr(parse(text=Tmp[2]))
Tmp=gsub(Xname,"x",Fstr)
Fstr=gsub(Yname,"y",Tmp)
Tmp=grep(Xname,Fstr)
if(Tmp==0){
Fstr=paste(Fstr,'+0*x',sep="")
}
Mdv=50
Ndv=50
if(Nargs>=4){
Tmp=Op(4,MS)
if(nrow(Tmp)>1){
Mdv=Tmp[1,1]
Ndv=Tmp[1,2]
}else{
Mdv=Tmp
if(Nargs==4){
Ndv=Mdv
}else{
Tmp1=Op(5,MS)
if((is.numeric(Tmp1))&&(length(Tmp1)==1)){
Ndv=Tmp1
}else{
Ndv=Mdv
}
}
}
}
X1=Xrange[1]; X2=Xrange[2]
Y1=Yrange[1]; Y2=Yrange[2]
Dx=(X2-X1)/Mdv
Dy=(Y2-Y1)/Ndv
Xval=c()
Tmp=X1
while(Tmp<X2+Eps){
Xval=c(Xval,Tmp)
Tmp=Tmp+Dx
}
Yval=c()
Tmp=Y1
while(Tmp<Y2+Eps){
Yval=c(Yval,Tmp)
Tmp=Tmp+Dy
}
I=1
for(y in Yv +al){
J=1
for(x in Xval){
Zval[I,J]=eval(parse(text=Fstr))
J=J+1
}
I=I+1
}
return(list(Zval,Xval,Yval))
} +Implicitplot<- function(...){
varargin=list(...)
Eps=10^(-4)
Nargs=length(varargin)
Fstr=varargin[[1]]
if(!is.numeric(Fst)){
MS=Fstr
Tmp=Makevaltable(MS)
Zval=Tmp[[1]]; Xval=Tmp[[2]]; Yval=Tmp[[3]]
}else if(is.character(Fstr)){
MS=list()
for(I in 1:Nargs){
MS=c(MS,varargin[i])
}
Tmp=Makevaltable(MS)
Zval=Tmp[[1]]; Xval=Tmp[[2]]; Yval=Tmp[[3]]
}else{
Zval=varargin[[1]]
Xval=varargin[[2]]
Yval=varargin[[3]]
}
Out=c()
for(J in Looprange(1,length(Yval)-1)){
for(I in Looprange(1,length(Xval)-1)){
a1=Xval[I]; b1=Yval[J]; c1=Zval[J,I]
a2=Xval[I+1]; b2=Yval[J]; c2=Zval[J,I+1]
a3=Xval[I+1]; b3=Yval[J+1]; c3=Zval[J+1,I+1]
a4=Xval(I);b4=Yval(J+1);c4=Zval(J+1,I);
PL=as.matrix(c(a1,b1,a2,b2,a3,b3,a4,b4,a1,b1),ncol=2,byrow=FALSE)
VL=c(c1,c2,c3,c4,c1)
QL=c()
for(K in 1:4){
if(abs(VL[K])<=Eps){
QL=appendrow(QL,PL[K,])
}else if(VL[K]>Eps){
if(VL[K+1]<-Eps){ +Tmp=1/(VL[K]-VL[K+1])*(-VL[K+1]*PL[K,]+VL[K]*PL[K+1,])
QL=appendrow(QL,Tmp)
}
}else{
if(VL[K+1]>Eps){
Tmp=1/(-VL[K]+VL[K+1])*(VL[K+1]*PL[K,]-VL[K]*PL[K+1,])
QL=appendrow(QL,Tmp)
}
}
}
if(nrow(QL)==2){
Out=append(Out,c(Inf,Inf),QL[1,],QL[2,])
}
}
}
Out=Out[2:nrow(Out),]
if(nrow(Out)>0){
Out=Connectseg(Out)
}else{
Out=c()
}
retuirn(Out)
} +Ketinit()##
Setwindow(c((-5),(4)),c((-5),(5)))####
A=c(1,-2);Assignadd('A',A) +imp1=Implicitplot("x^2+x*y+2*y^2-(1)","x=c(-3,3)","y=c(-3,3)") +source('~/Dropbox/kettoday/0930bezierimplicitR/s0212implicitimp1.r', chdir = TRUE) +seq(0,0) +?matrix +dt=matrix(nrow=5,ncol=4) +dt +source('~/Dropbox/kettoday/0930bezierimplicitR/s0212implicitimp1.r', chdir = TRUE) +length(matrix(1:5,ncol=2)) +length(matrix(1:6,ncol=2)) +source('~/Dropbox/kettoday/0930bezierimplicitR/s0212implicitimp1.r', chdir = TRUE) +imp1 +source('~/Dropbox/kettoday/0930bezierimplicitR/s0212implicitimp1.r', chdir = TRUE) +imp1 +Windisp(imp1) +source('~/Dropbox/kettoday/0930bezierimplicitR/s0212implicitimp1.r', chdir = TRUE) +Windisp(imp1) +source('~/Dropbox/kettoday/0930bezierimplicitR/s0212implicitimp1.r', chdir = TRUE) +Windisp(imp1) +source('~/Dropbox/kettoday/0930bezierimplicitR/s0212implicitimp1.r', chdir = TRUE) +source('~/Desktop/Rproject2017-/171001.r', chdir = TRUE) +funP +fun(1,0.5,2) +funP(1,0.5,2) +funN +source('~/Desktop/Rproject2017-/171001.r', chdir = TRUE) +funP +funP(1,0.5,2) +tmp="N=200" +tmp1="N=200" +tmp=strsplit(tmp1,"=") +tmp +tmp=tmp[[1]] +tmp +tmp1="N=200" +tmp=strsplit(tmp1,"=")[[1]] +tmp +tmp=strsplit(tmp1,"=") +tmp +tmp[1] +tmp[[1]] +tmp=list(1,5,6,3) +tmp[1,2] +tmp[1:2] +tmp[[1:2]] +tmp +tmp=list(1,5,6,3) +tmp[3:4 +] +tmp +tmp[3:4] +tmp[[3:4]] +tmp[[3]] +tmp[[3:4]] +source('~/Dropbox/kettoday/0930bezierimplicitR/s0212implicitimp1.r', chdir = TRUE) +dt1=c(3.055,3.012,3.045,3.033,3.070,3.038,3.052,3.044,3.072,3.005) +mean(dt1) +var(dt1) +var(dt1)*0.9 +d_1=5;d_1 +varargin=list('c(yN1,yN2)`=c(yN2,-yN2-4*yN1)','x=c(-5,4.991263)',0,c(1,0),'Num=200') +deq=varargin[[1]]
rng=varargin[[2]]
initt=varargin[[3]]
initf=varargin[[4]]
Num=50 +tmp=strsplit(varargin[[5]],"=")[[1]]
Num=eval(parse(text=tmp[2]))
} +if(length(varargin)>4){
tmp=strsplit(varargin[[5]],"=")[[1]]
Num=eval(parse(text=tmp[2]))
} +Num +tmp=strsplit(deq,"=")[[1]]
tmp1=grep(tmp[1],"`","",fixed=TRUE)
tmp1=substring(tmp1,1,nchar(tmp1)-1) +tmp1 +tmp=strsplit(deq,"=")[[1]] +tmp +tmp1=grep(tmp[1],"`","",fixed=TRUE);tmp1 +tmp1=gsub(tmp[1],"`","",fixed=TRUE) +tmp1 +tmp[1] +tmp1=gsub(tmp[1],"`","") +tmp1 +tmp1=gsub("`","",tmp[1],fixed=TRUE);tmp1 +tmp1=substring(tmp1,2,nchar(tmp1)-1);tmp1 +tmp1=gsub("`","",tmp[1],fixed=TRUE)
tmp1=substring(tmp1,3,nchar(tmp1)-1);tmp1 +Xname=strstplit(tmp1,",") +Xname=strsplit(tmp1,",");Xname +Xname=strsplit(tmp1,",")[[1]];Xname +for(J in 1..(length(Xname))){
func=gsub(Xname[J],"X["+as.character(J)+"]",Xname,fixed=TRUE)
} +for(J in 1:(length(Xname))){
func=gsub(Xname[J],"X["+as.character(J)+"]",Xname,fixed=TRUE)
} +as.character(5) +for(J in 1:(length(Xname))){
tmp=parse("X[",as.character(J),"]",sep="")
func=gsub(Xname[J],tmp,Xname,fixed=TRUE)
} +func=tmp[2]
for(J in 1:(length(Xname))){
tmp=paste("X[",as.character(J),"]",sep="")
func=gsub(Xname[J],tmp,Xname,fixed=TRUE)
} +func +for(J in 1:(length(Xname))){
tmp=paste("X[",as.character(J),"]",sep="")
print(tmp)
func=gsub(Xname[J],tmp,Xname,fixed=TRUE)
} +func=tmp[2]
for(J in 1:(length(Xname))){
tmp=paste("X[",as.character(J),"]",sep="")
print(tmp)
func=gsub(Xname[J],tmp,Xname,fixed=TRUE)
} +func=tmp[2]
for(J in 1:(length(Xname))){
tmp=paste("X[",as.character(J),"]",sep="")
func=gsub(Xname[J],tmp,func,fixed=TRUE)
} +func +tmp=strsplit(deq,"=")[[1]]
tmp1=gsub("`","",tmp[1],fixed=TRUE)
tmp1=substring(tmp1,3,nchar(tmp1)-1)
Xname=strsplit(tmp1,",")[[1]]
func=tmp[2]
for(J in 1:(length(Xname))){
tmp1=paste("X[",as.character(J),"]",sep="")
func=gsub(Xname[J],tmp1,func,fixed=TRUE)
} +func +tmp=strsplit(rng,"=")[[1]]
tname=tmp[1]
tmp=eval(parse(text=tmp[2]))
t1=tmp[1]
t2=tmp[2] +t1 +t2 +tmp=paste("funP=function(",tname,",X){",func,sep="") +tmp +tmp=paste("funP=function(",tname,",X){",func,"}",sep="");tmp +tmp=paste("funP=function(",tname,",X){-",func,"}",sep="") +tmp +?flatten +tt=initt
X0=initf
pdL=c(tt,X0)
;pdL +tt=initt
X0=initf
pdL=c(tt,X0) +pdL +floor(1.5) +flg=0
for(J in 1:(floor((t2-initt)/dt))){
if(flg==0){
kl1=dt*funP(tt,X0)
kl2=dt*funP(tt+dt/2,X0+kl1/2)
kl3=dt*funP(tt+dt/2,X0+kl2/2)
kl4=dt*funP(tt+dt,X0+kl3)
X0=X0+(kl1+2*kl2+2*kl3+kl4)/6
tt=initt+#*dt
pdL=Appendrow(pdL,c(tt,X0))
if(abs(tmp)>inf){flg=1}
}
} +t2 +initt +dt +dt=(t2-t1)/Num
tt=initt
X0=initf
pdL=c(tt,X0)
flg=0 +dt +for(J in 1:(floor((t2-initt)/dt))){
if(flg==0){
kl1=dt*funP(tt,X0)
kl2=dt*funP(tt+dt/2,X0+kl1/2)
kl3=dt*funP(tt+dt/2,X0+kl2/2)
kl4=dt*funP(tt+dt,X0+kl3)
X0=X0+(kl1+2*kl2+2*kl3+kl4)/6
tt=initt+#*dt
pdL=Appendrow(pdL,c(tt,X0))
if(abs(tmp)>inf){flg=1}
}
} +funP +tmp=paste(function(",tname,",X){",func,"}",sep="")
funP=eval(parse(text=tmp)) +funP +tmp=function(",tname,",X){-",func,"}",sep="")
funN=eval(parse(text=tmp)) +tmp=paste(function(",tname,",X){",func,"}",sep="") +tmp=paste("function(",tname,",X){",func,"}",sep="")
funP=eval(parse(text=tmp)) +funP +tmp=paste("function(",tname,",X){-",func,"}",sep="")
funN=eval(parse(text=tmp)) +funN +dt=(t2-t1)/Num
tt=initt
X0=initf
pdL=c(tt,X0)
flg=0
for(J in 1:(floor((t2-initt)/dt))){
if(flg==0){
kl1=dt*funP(tt,X0)
kl2=dt*funP(tt+dt/2,X0+kl1/2)
kl3=dt*funP(tt+dt/2,X0+kl2/2)
kl4=dt*funP(tt+dt,X0+kl3)
X0=X0+(kl1+2*kl2+2*kl3+kl4)/6
tt=initt+#*dt
pdL=Appendrow(pdL,c(tt,X0))
if(abs(tmp)>inf){flg=1}
}
} +source('~/ketcindy/ketlib/ketpic5_2_3.r', chdir = TRUE) +de2=Deqplot('c(yN1,yN2)`=c(yN2,-yN2-4*yN1)','x=c(-5,4.991263)',0,c(1,0),'Num=200') +source('~/ketcindy/ketlib/ketpic5_2_3.r', chdir = TRUE) +de2=Deqplot('c(yN1,yN2)`=c(yN2,-yN2-4*yN1)','x=c(-5,4.991263)',0,c(1,0),'Num=200') +de2 +source('~/ketcindy/ketlib/ketpic5_2_3.r', chdir = TRUE) +de2=Deqplot('c(yN1,yN2)`=c(yN2,-yN2-4*yN1)','x=c(-5,4.991263)',0,c(1,0),'Num=200') +Norm(c(1,5,3)) +Norm(c(1,5,3,4)) +source('~/ketcindy/ketlib/ketpic5_2_3.r', chdir = TRUE) +de2=Deqplot('c(yN1,yN2)`=c(yN2,-yN2-4*yN1)','x=c(-5,4.991263)',0,c(1,0),'Num=200') +source('~/ketcindy/ketlib/ketpic5_2_3.r', chdir = TRUE) +de2=Deqplot('c(yN1,yN2)`=c(yN2,-yN2-4*yN1)','x=c(-5,4.991263)',0,c(1,0),'Num=200') +de2 +tmp="abc'" +tmp1=gsub("c'","",tmp);tmp1 +tmp="abc)'" +tmp1=gsub("c'","",tmp);tmp1 +tmp1 +tmp1=gsub(")'","",tmp);tmp1 +tmp1=gsub(")'","","c(1,2)'");tmp1 +gsub(")'","","c(1,2)'") +gsub(")'","","c(1,2)`") +gsub(")`","","c(1,2)`") +source('~/Dropbox/kettoday/1005deqplot/Deqplot171006.r', chdir = TRUE) +tmp +Windisp(tmp) +Windisp(tmp[,2:3]) +source('~/Dropbox/kettoday/1005deqplot/Deqplot171006.r', chdir = TRUE) +Windisp(tmp[,2:3]) +source('~/Dropbox/kettoday/1005deqplot/Deqplot171006.r', chdir = TRUE) +Windisp(tmp[,2:3]) +source('~/Dropbox/kettoday/1005deqplot/Deqplot171006.r', chdir = TRUE) +Windisp(tmp[,2:3]) +?sprintf +source('~/Desktop/diffofsci2e/ketpic2ereplaced.r', chdir = TRUE) +Beginpicture("1cm") +g1=Plotdata("sin(x)","x") +g1 +Drwline(g1) +Dottedline(g1) +Dashline(g1) +source('~/Desktop/diffofsci2e/ketpic2ereplaced.r', chdir = TRUE) +Dashline(g1) +Setpen(2) +Drwpt(list(c(0,0),c(1,2))) +source('~/ketcindy/ketfiles/ketwork/template.r', chdir = TRUE) +source('~/ketcindy/ketfiles/ketwork/s0505bezier3d.r', chdir = TRUE) +source('~/ketcindy/ketfiles/ketwork/s1501mant.r', chdir = TRUE) +source('~/ketcindy/ketfiles/ketwork/s0505bezier3d.r', chdir = TRUE) +Windisp(bz2d) +bz3d +bz2d +source('~/ketcindy/ketfiles/ketwork/s0505bezier3d.r', chdir = TRUE) +bz3d +bz2d +Setangle() +source('~/ketcindy/ketfiles/ketwork/s0505bezier3d.r', chdir = TRUE) +source('~/ketcindy/ketfiles/ketsample/samples/s07slides/fig/s0706animateall.r', chdir = TRUE) |