Solution of General Partial Differential Equations

This script gives the numirical solution of PDE of the form a1*Uxx + a2*Uyy + a3*Uxy + a4*Ux+ a5*Uy + a6*U=f(x,y) where a1,a2,a3,a4,a5 and a6 are functions on x and y or constants. The solustion is carried by iteration method

Number of intervals = Coefficient of Uxx a1(X,Y)= Coefficient of Uyy a2(X,Y)= Coefficient of Uxy a3(X,Y)= Coefficient of Ux a4(X,Y)= Coefficient of Uy a5(X,Y)= Coefficient of U a6(X,Y)= the RHS f(x,y)= u(X0,Y)= u(Xn,Y)= u(X,Y0)= u(X,Yn)= Lower X0=Lower Y0 = Upper Xn=Upper Yn = Output format :