testwafer.m

MDPSAS工具箱是马里兰大学开发的

clear all
clear classes

addpath ../../class
addpath ../../class/double
addpath ../../../PhysProp

clear all

N = 30;    % quadrature pts
w = wafer(N);

w = set(w,'solverset',20,'nmax');
w = set(w,'solverset',1,'plot');
w = newton(w);
plot(w)

xlabel('Radial position r (m)')
ylabel('Temperature (K)')
title('Steady state solution')