Demande de coup de main

Bonsoir j'aimerais que vous m'aidiez à programmer sur Scilab la résolution des équations différentielles:
1){█(∂u/∂t- (∂^2 u)/〖∂x〗^2 =0; 0≤x≤1;0≤t≤0,2@u(x,0)=sin⁡(πx)+x@u(0,t)=0@u(1,t)=1)┤
Cherchons la solution sous la forme uh sous la forme
uh(x, t) = sin⁡(πx)+x+∑_(j=1)^N▒a_j(t) (x^j-x^(j+1) )
2){█(μ^2 (∂^2 u)/(∂x^2 )+(∂^2 u)/(∂y^2 )+1=0 (x,y)ϵ [-1,1]^2@u(-1 ; y)= u(1 ; y)= 0 -1 ≤y ≤1@u(x ; -1)= u(x; 1)= 0 -1≤x≤1)┤
Considérons l’équation différentielle μ^2 (∂^2 u)/(∂x^2 )+(∂^2 u)/(∂y^2 )+1=0 ⟺μ^2 (∂^2 u)/(∂x^2 )+(∂^2 u)/(∂y^2 )=-1 (1)
Cherchons une approximation de u_h de u sous la forme :
u_h (x,y)= ∑_(j=1)^N▒〖a_j (1-x^2 )^j (1-y^2 〗 )^j