function lpc=lpconst(x) if 3*x(1)+x(2)-10=-0.5 lpc=1; else lpc=0; end
http://www.codebf.com/read/4/121
function [sol,r1,r2]=randlp(a,b,n) debug=1; a=0; b=10; n=1000; r1=unifrnd(a,b,n,1); r2=unifrnd(a,b,n,1); sol=[r1(1) r2(1)]; z0=inf; for i=1:n x1=r1(i); x2=r2(i); lpc=lpconst([x1
http://www.codebf.com/read/4/123
function dy=shier2(t,y) dy=zeros(2,1); dy(1)=y(1)*(0.9-0.1*y(2)); dy(2)=y(2)*(-0.6+0.02*y(1));
http://www.codebf.com/read/4/124
function dy=eq1(x,y) dy=zeros(2,1); dy(1)=y(2); dy(2)=1/5*sqrt(1+y(1)^2)/(1-x);
http://www.codebf.com/read/4/125
x0=0;xf=0.99999; [x,y]=ode23('eq1',[x0 xf],[0 0]); Y=0:0.01:2; plot(1,Y,'g.') hold on plot(x,y(:,1),'b*')
http://www.codebf.com/read/4/126
function dx=shier1(t,x) dx=zeros(2,1); dx(1)=x(1)*(0.7-0.1*x(2)); dx(2)=x(2)*(-0.8+0.02*x(1));
http://www.codebf.com/read/4/127
y=dsolve('D2y+4*Dy+29*y=0','y(0)=0,Dy(0)=15','x')
http://www.codebf.com/read/4/128
t0=0;tf=10; [t,y]=ode45('eq3',[t0 tf],[0 0]); T=0:0.1:2*pi; X=10+20*cos(T); Y=20+15*sin(T); plot(X,Y,'-') hold on plot(y(:,1),y(:,2),'*')
http://www.codebf.com/read/4/129
dsolve('Du=1+u^2','t')
http://www.codebf.com/read/4/130
function dy=vdp1000(t,y) dy=zeros(2,1); dy(1)=y(2); dy(2)=1000*(1-y(1)^2)*y(2)-y(1);
http://www.codebf.com/read/4/131
虫虫下载站 半导体技术网 电子研发网 源码地带 电源技术网 单片机技术网 医疗电子技术 嵌入式系统与单片机