diff options
| author | priyanka | 2015-06-24 15:03:17 +0530 |
|---|---|---|
| committer | priyanka | 2015-06-24 15:03:17 +0530 |
| commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
| tree | ab291cffc65280e58ac82470ba63fbcca7805165 /800/CH10/EX10.3/10_3.sce | |
| download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip | |
initial commit / add all books
Diffstat (limited to '800/CH10/EX10.3/10_3.sce')
| -rwxr-xr-x | 800/CH10/EX10.3/10_3.sce | 50 |
1 files changed, 50 insertions, 0 deletions
diff --git a/800/CH10/EX10.3/10_3.sce b/800/CH10/EX10.3/10_3.sce new file mode 100755 index 000000000..cb15d3f20 --- /dev/null +++ b/800/CH10/EX10.3/10_3.sce @@ -0,0 +1,50 @@ +//clear//
+clc
+clear
+exec("10.3data.sci");
+w = 0:10:10000;
+
+function W=f(w,x)
+
+ W =zeros(1,1);
+
+pt0=.3*Po;
+y=(1-alpha*w)^.5;
+ph=pt0*(1.5-x)*y;
+pt=pt0*(1-x)*y;
+pb=2*pt0*x*y;
+rt=-k*kt*ph*pt/(1+kb*pb+kt*pt);
+rate=-rt;
+W(1)=-rt/ftO;
+endfunction
+pt0=.3*Po;
+X=ode([0],w0,w,f);
+
+
+for i =1:length(X)
+ y(1,i)=(1-alpha*w(1,i))^.5;
+ph(1,i)=pt0*(1.5-X(1,i))*y(1,i);
+pt(1,i)=pt0*(1-X(1,i))*y(1,i);
+pb(1,i)=2*pt0*X(1,i)*y(1,i)
+end
+
+m1 = X';
+m2=y';
+scf(1)
+plot2d(w',[m1 m2]);
+
+xtitle( 'Figure E10-3.1', 'w', 'x,y' ) ;
+legend(['x';'y']);
+
+scf(2)
+l1=ph'
+l2=pt'
+l3=pb'
+plot2d(w',[l1 l2 l3]);
+
+xtitle( 'Figure E10-3.2', 'w', 'ph,pt,pb' ) ;
+legend(['ph';'pt';'pb']);
+
+
+
+
|