function y = pb1(n,x)
	y = sin(n*x*%pi);
endfunction

function y = expb1(x,t)
	y = exp(-t) * sin(2*x*%pi);
endfunction

function res = graph(n)
	x = -1:0.05:1;
	plot2d(x,pb1(n,x));
endfunction

function y = graph2()
	x = -1:0.1:1;
	t = 0:0.05:1;
	y = expb1(x,t');
	plot3d(x,t,y);
endfunction


function y = graphf()
	x = -1:0.1:1;
	y = x^2 - %pi^2;
	plot2d(x,y);
endfunction
