clf()
ylim(0,.5)
ylim(0,1)
t=linspace(0,2*pi)
y=sin(t)
h1=axes([.1,.1,.8,.8])
plot(t,y)
h2=axes([0.6,0.6,0.2,0.2])
fill(t,y**2,'r')
h3=axes([0.15,0.2,.3,.3])
polar(0.25*t,0.25*y/t)

Strang-6.4-1 (last edited 2010-12-17 12:45:49 by anvesh)