summaryrefslogtreecommitdiff
path: root/texmf-dist/doc/support/ketcindy/ketsample/samples/s07slides/fig/.Rhistory
diff options
context:
space:
mode:
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/.Rhistory238
1 files changed, 0 insertions, 238 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
deleted file mode 100644
index 8470a278..00000000
--- a/texmf-dist/doc/support/ketcindy/ketsample/samples/s07slides/fig/.Rhistory
+++ /dev/null
@@ -1,238 +0,0 @@
-?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)