代码搜索结果

function lpc=lpconst(x) if 3*x(1)+x(2)-10=-0.5 lpc=1; else lpc=0; end

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

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

function dy=eq1(x,y) dy=zeros(2,1); dy(1)=y(2); dy(2)=1/5*sqrt(1+y(1)^2)/(1-x);

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*')

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

y=dsolve('D2y+4*Dy+29*y=0','y(0)=0,Dy(0)=15','x')

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

dsolve('Du=1+u^2','t')

function dy=vdp1000(t,y) dy=zeros(2,1); dy(1)=y(2); dy(2)=1000*(1-y(1)^2)*y(2)-y(1);