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 --- 29/CH1/EX1.6.9/exa1_6_9.sce | 23 +++++++++++++++++++++++ 1 file changed, 23 insertions(+) create mode 100755 29/CH1/EX1.6.9/exa1_6_9.sce (limited to '29/CH1/EX1.6.9/exa1_6_9.sce') diff --git a/29/CH1/EX1.6.9/exa1_6_9.sce b/29/CH1/EX1.6.9/exa1_6_9.sce new file mode 100755 index 000000000..7895a9897 --- /dev/null +++ b/29/CH1/EX1.6.9/exa1_6_9.sce @@ -0,0 +1,23 @@ +//Caption:initial_values +// example 1.6.7 +//page 13 +//F(s)=(4*s+1)/(s^3+2*s) +s=%s; +syms t; +F=(4*s+1)/(s^3+2*s) +f = ilaplace (F,s,t); +y=limit(f,t,0);//initial value theorem +y=dbl(y); +disp(y,"f(0+)=") +// since F'(s)=sF(s)-f(0+) where L(f'(t))=F'(s)=F1 +F1=(4*s+1)/(s^2+2) +f1= ilaplace(F1,s,t); +y1=limit(f1,t,0);//initial value theorem +y1=dbl(y1); +disp(y1,"f_prime(0+)=") +// since F''(s)=(s^2)*F(s)-s*f(0+)-f'(0+) where L(f''(t))=F''(s)=F2 +F2=(s-8)/(s^2+2) +f2= ilaplace(F2,s,t); +y2=limit(f2,t,0);//initial value theorem +y2=dbl(y2); +disp(y2,"f_doubleprime(0+)=") \ No newline at end of file -- cgit