Scribd
Upload
6 views4 pages

One-Dimensional Wave Equation Solution

The document outlines the solution to the one-dimensional wave equation Utt=c^2(Uxx) for an infinite string problem. It provides the specific solution formula and initial conditions, with examples of plotting the wave function at different time intervals. Additionally, it includes commands for visualizing the wave function in both 2D and 3D plots.

Uploaded by

abhijit3856
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
6 views4 pages

One-Dimensional Wave Equation Solution

The document outlines the solution to the one-dimensional wave equation Utt=c^2(Uxx) for an infinite string problem. It provides the specific solution formula and initial conditions, with examples of plotting the wave function at different time intervals. Additionally, it includes commands for visualizing the wave function in both 2D and 3D plots.

Uploaded by

abhijit3856
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd

Practical-4 wave equation.

wxmx 1 / 4

To find the solution of


one dimensional wave
equation Utt=c^2(Uxx)
Part(a)is an infinite
string problem since x
belongs to R.
Its solution is given by:
u(x,t):=(f(x+c*t)+f(x-c*t))/2+(1/(2*c))*(integrate(g(a),a,x-c*t,x+c*t));

Utt-Uxx=0, U(x,0)=sin(x), Ut(x,0)=x^2, inf< x <


inf ,t>0.

➔ kill(all)$

➔ u(x,t):=(f(x+c·t)+f(x−c·t))/2+(1/(2·c))·(integrate(g(a),a,x−c·t,x+c·t));
Z x+c t
f ( x + c t )+ f ( x − c t ) 1
(%o1) u ( x , t ) := + g (a)
2 2c
x−c t

➔ f(x):=sin(x)$

➔ g(x):=x^2$

➔ c:1$

➔ u(x,t);
sin ( x + t ) + sin ( x − t )
(%o5) +
2
3 2 2 3 3 2 2 3
x +3 t x +3 t x+t x −3 t x +3 t x−t
−
3 3
2
Practical-4 wave [Link] 2 / 4

➔ wxdraw2d(explicit(u(x,0),x,−2,2),title="Wave at t=0");

(%t6)

(%o6)

➔ wxdraw2d(explicit(u(x,1),x,−2,2),title="Wave at t=1");

(%t7)

(%o7)
Practical-4 wave [Link] 3 / 4

➔ wxdraw2d(explicit(u(x,2),x,−2,2),title="Wave at t=2");

(%t8)

(%o8)

➔ wxdraw2d([color=red,explicit(u(x,0),x,−2,2),color=blue,
explicit(u(x,1),x,−2,2),color=green,explicit(u(x,2),x,−2,2)]);

(%t9)

(%o9)
Practical-4 wave [Link] 4 / 4

➔ wxplot3d(u(x,t),[x,−8,8],[t,0,5]);

(%t10)

(%o10)

➔ wxdraw3d(explicit(u(x,t),x,−8,8,t,0,5));

(%t11)

(%o11)

You might also like