Monday, July 15, 2013

(dr)/(d theta) = sin^4(pitheta) Solve the differential equation.

Separate the r and theta variables and integrate both sides:
int dr = int sin(pi*theta)^4 d(theta)
r + c_1=int sin(pi*theta)^4 d(theta)
Let pi*theta=u, and pi*d(theta)=du
r + c_1=int sin(u)^4 d(theta)
r + c_1=(1/pi)int sin(u)^4 du
r+c_1=(1/pi)int sin(u)^2*sin(u)^2 du
Use a trigonometric identity:
r+c_1=(1/pi) int (1/2)(1-cos(2u))(1/2)(1-cos(2u)) du
r+c_1=1/(4pi)int (1-2cos(2u)+cos(2u)^2) du
Let 2u=t, and 2du=dt
r+c_1=1/(4pi) int (1-2cos(t)+cos(t)^2)(dt/2)
r+c_1=1/(8pi) int (1-2cos(t)+(1/2)(1+cos(2t))) dt
r+c_1=1/(8pi)[int 1dt-2int cos(t) dt+(1/2)int dt+(1/2)int cos(2t) dt]
r+c_1=1/(8pi)[t-2sin(t)+1/2t+(1/2)(1/2)sin(2t)+c_2]
Substitute back in t=2u=2(pi*theta)
r+c_1=1/(8pi)[2pi*theta-2sin(2pi*theta)+pi*theta+(1/4)sin(4pi*theta)+c_2]
r+c_1=(2pi*theta)/(8pi)-2sin(2pi*theta)/(8pi)+(pi*theta)/(8pi)+1/(32pi)sin(4pi*theta)+c_2/(8pi)
Simplify terms and combine the constants c_1 and c_2 into a new constant c :
r=theta/(4)-1/(4pi)sin(2pi*theta)+theta/8+1/(32pi)sin(4pi*theta)+c
Your general solution is then:
r(theta)=(3pi)/4-1/(4pi)sin(2pi*theta)+1/(32pi)sin(4pi*theta)+c
Below is a plot from 0 to 2pi of a particular solution when c=0 .

No comments:

Post a Comment

Summarize the major research findings of "Toward an experimental ecology of human development."

Based on findings of prior research, the author, Bronfenbrenner proposes that methods for natural observation research have been applied in ...