analytisk:=proc(u,f)
local v,h,diffekv, solution,z;
if simplify(diff(u(x,y),x,x)+diff(u(x,y),y,y))<>0 then `u ej harmonisk` else
v:=int(-diff(u(x,y),y),x);
v:=v+h(y);
v:=unapply(v,x,y);
diffekv:=simplify(diff(v(x,y),y)-diff(u(x,y),x));
solution:=dsolve(diffekv,h(y));
h:=unapply(rhs(solution),y); #rhs står för right hand side.
print(`Imaginärdelen:`, v(x,y));
f:=unapply(u(x,y)+I*v(x,y),x,y);
f:=unapply(f(z,0),z):
print(`Den analytiska funktionen:`,f(z));

end if;
end proc;