From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001 From: priyanka Date: Wed, 24 Jun 2015 15:03:17 +0530 Subject: initial commit / add all books --- 564/CH18/EX18.2/18_2.sce | 37 +++++++++++++++++++++++++++++++++++++ 1 file changed, 37 insertions(+) create mode 100755 564/CH18/EX18.2/18_2.sce (limited to '564/CH18/EX18.2/18_2.sce') diff --git a/564/CH18/EX18.2/18_2.sce b/564/CH18/EX18.2/18_2.sce new file mode 100755 index 000000000..360a186f9 --- /dev/null +++ b/564/CH18/EX18.2/18_2.sce @@ -0,0 +1,37 @@ +pathname=get_absolute_file_path('18_2.sce') +filename=pathname+filesep()+'18_2data.sci' +exec(filename) +clear +del=2*((b/tb)+(a/ta));//δ +A=a*b; +deff("[W41]=f(s1)","W41=(T/(2*A*G))*((s1/tb)+(del*a*s1/(4*A)))"); +deff("[W12]=f1(s2)","W12=(T/(2*A*G))*((s2/ta)+(del*b*s2/(4*A)))"); +deff("[W23]=f2(s1)","W23=-((T/(2*A*G))*((s1/tb)+(del*a*s1/(4*A))))"); +deff("[W34]=f3(s2)","W34=-((T/(2*A*G))*((s2/ta)+(del*b*s2/(4*A))))"); +deff("[Sz]=f4(s1)","Sz=0*s1"); +deff("[Sz1]=f5(s2)","Sz1=0*s2"); +funcprot(0); +s1=[-b/2:0.05:b/2]; +s2=[-a/2:0.05:a/2]; +subplot(2,2,1) +fplot2d(s1,f) +fplot2d(s1,f4) +xgrid(3); +xtitle( 'Wraping in 4-1', ' -y- ', '-w-'); +subplot(2,2,2) +fplot2d(s2,f1) +fplot2d(s2,f5) +xgrid(3); +xtitle( 'Wraping in 1-2', ' -x- ', '-w-'); +subplot(2,2,3) +fplot2d(s1,f2) +fplot2d(s1,f4) +xgrid(3); +xtitle( 'Wraping in 2-3', ' -y- ', '-w-'); +subplot(2,2,4) +fplot2d(s2,f3) +fplot2d(s2,f5) +xgrid(3); +xtitle( 'Wraping in 3-4', ' -x- ', '-w-'); +datatipToggle(); +printf("\nclick on the point on the plot to view its coordinates") -- cgit