for i in xrange(len(xs)): xy = np.array((xs[i], ys[i]), dtype=np.double) uv = timoBeam.displacement(xy[0], xy[1], 0, None) if is_vs: exp_uvs.append(uv[1]) else: exp_uvs.append(uv[0]) if None is axis: axis = xs ax.plot(axis, act_uvs, "r", label="Numerical") ax.plot(axis, exp_uvs, "b", label="Precise") if __name__ == "__main__": WeakformProcessor2D = JClass("net.epsilony.tsmf.process.WeakformProcessor2D") processor = WeakformProcessor2D.genTimoshenkoProjectProcess() processor.process() processor.solve() pp = processor.postProcessor() fig = plt.figure() ax = fig.add_subplot(221) ax.set_title("displacement $v$ along axis $x$ ($y=0$)") ax.set_xlabel("$x$") ax.set_ylabel("$v(x,0)$") y = 0 plot_vs_by_y(pp, ax, y) ax2 = fig.add_subplot(222) ax2.set_title("displacement $u$\n along the left edge ($x=0$)") ax2.set_xlabel("$y$")