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, 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) |