Question: I have a situation in FEM..

i am deriving lagrangian shape function symbolically for an 1-D element of order p in FEM....so the eleemnt will have p+1 nodes....

                 

                                   N(k,i) = product((X-Xother)) / product(Xcurrent-Xother)

where: Xcurrent = position of current node for which shape function is derived

            Xother = position of other nodes of the element which i already have in a vector.

             k =number of elements

             i = number of nodes in k th element...

i am having difficulty in generalising the formula for order p.....i.e putting the formula in 'for' loops ...

thanks.

                        

Please Wait...