напиши код на c++ для решения
(w_m^(s+1)-2*w_m^s+w_m^(s-1))/τ^2 +γ (w_m^(s+1)-w_m^s)/τ+ξ (w_(m+2)^(s+1)-4*w_(m+1)^(s+1)+6*w_m^(s+1)-w_(m-1)^(s+1)+w_(m-2)^(s+1))/h^4 =(w_(m+1)^(s+1)-2*w_m^(s+1)+w_(m-1)^(s+1))/h^2 +α*p_m
(c_m^(s+1)-2*c_m^s+c_m^(s-1))/τ^2 +γ (c_m^(s+1)-c_m^s)/τ=-ξ〖*λ〗_m^2*c_m^(s+1)-λ_m*c_m^(s+1)+α*p_m
c_m^(s+1)=(c_m^s (2+γ*τ)-c_m^(s-1)-α*p_m*τ^2)/(1+γ*τ+τ^2 λ_m+ξ*τ^2*λ_m^2 )
p_m=1/(1-w_m^s )^2 
λ_m=-4*sin^2⁡〖((π*h*m)/2)〗
α=1.0034; ξ=0.002729; γ=0.1
τ=h=0.01
