?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(TmpEps){ 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)