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 --- 50/CH7/EX7.3/ex_7_3.sce | 29 +++++++++++++++++++++++++++++ 1 file changed, 29 insertions(+) create mode 100755 50/CH7/EX7.3/ex_7_3.sce (limited to '50/CH7/EX7.3') diff --git a/50/CH7/EX7.3/ex_7_3.sce b/50/CH7/EX7.3/ex_7_3.sce new file mode 100755 index 000000000..6c8f5c1c5 --- /dev/null +++ b/50/CH7/EX7.3/ex_7_3.sce @@ -0,0 +1,29 @@ +// example 7.3 +// solve by shooting method; + +// u''=2*u*u'; +// u(0)=0.5; u(1)=1; + +// let -> U1(x)=du/dx; +// U2(x)=d2u/dx2; + +// U(x)=[U1(x);U2(x)] + +// hence ; +// dU/dx=f(x,U); + +h=.25; + +ub=[.5,1]; + +up=[0:.1:1]; + +x=0:h:1; + +deff('[w]=f(x,U)','w=[U(2); 2*U(1)*U(2)]') + + + +[U] = shooting(ub,up,x,f); + +// the solution obtained would show the values of u in the first collumn and their corresponding derivatives in the second collumn ; \ No newline at end of file -- cgit